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1,0 INTRODUCTION 


It is the intent of this handbook to provide design engineers with a com- 
prehensive guide in a convenient format for design of elastomer dampers 
for application in rotating machinery. Design data has been collected from 
an extensive group of references and, together with some recent information, 
assembled into a single volume. This handbook, which is intended to be a 
self-contained reference, includes some theoretical discussion where 
necessary to insure that the reader has an understanding of how and why 
dampers are used for rotating machines. A step-by-step procedure for the 
design of elastomer dampers is presented along with detailed examples of 
actual elastomer damper applications. In addition, numerous references are 
cited for the interested reader. The authors assume that the reader has a 
basic understanding of rotor dynamics and machinery vibration. 

The design information presented is primarily intended for application in 
rotating machinery. Dampers of this type have been shown to be effective 
for controlling rotor vibrations; however, an historic lack of design data 
and designer experience with elastomer dampers has inhibited the expansion 
of their applications. A primary motivation for this handbook is to pro- 
vide such data in a compact and convenient form and to encourage the use of 
elastomer dampers and elastomeric damping treatments, 

A portion of the design data contained in this handbook is relevant to many 
mechanical design problems and applications of elastomers. For example, 
the elastomer dynamic property data may be utilized in the design of 
structural supports for shock and vibration isolation and insulation. Some 
basic information regarding the effect of structural damping is included in 
this work. However, specific information concerning the design of such 
structural members is beyond the scope of this handbook, but an extensive 
bibliography is provided at the conclusion of each chapter to aid the 
design engineer in finding a solution to the problem. 

The term "elastomer" refers to a large variety of synthesized polymers as 
well as to natural rubber. These materials are capable of relatively 
large strains with essentially full recovery. This property of high elas- 
ticity (from which the term elastomer is derived) is the result of a parti- 
cular type of molecular structure. As obtained from Reference 1.1, the 
principal attributes of this type of molecular structure are: 

« The molecules must be very long and flexible, or to be more 
precise, there must be free rotation about most of the bonds 
joining together neighboring chain atoms; 

$ The molecules must be attached together here and there by 

either permanent chemical bonds or linkages, or by mechanical 
entanglements, either form of attachment being termed a 
cross-link; 

$ Apart from these cross-links, the molecules must be free to 
move readily past one another. This is to say, the inter- 
molecular attractions which exist between all molecules must 
be small. 
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In addition to resiliency, elastomers provide energy dissipation which is 
the key to their ability to convert mechanical energy into heat energy, 
thereby damping* 

Elastomers may be separated into the two general classifications of natural 
rubber and synthetic polymers. Particular synthetic polymers may be identi- 
fied either by their generic type (such as fluorocarbon, chloroprene, etc.) 
or commonly used manufacturer trademarks fl . 2^j*. Specific compositions of 
particular synthetic polymers are further identified by surface hardness. 

This surface hardness, which is generally measured in durometer (or Shore ^ 

A hardness) , is an indication of the static modulus of elasticity (or of 
the static shear modulus) of an elastomer. The specific elastomer types 
for which dynamic properties are presented in this handbook are polybuta- 
diene, 70 and 90 durometer Viton (fluorocarbon elastomer), Buna-N, EPDM, 
and neoprene. Physical properties of these and many other types of elasto- 
mers are also presented in this handbook. These physical properties 
include all material properties (e.g., density and thermal conductivity) 
except for the dynamic moduli and dynamic shape factors which are referred 
to as dynamic properties. An extensive list of references are given from 
which additional material properties for a large variety of elastomer types 
may be obtained. Also, test methods are described in detail by which 
dynamic properties of additional elastomer materials may be determined. 

A number of factors must be taken into consideration in elastomer material 
selection when designing elastomer dampers for rotating machines. The 
dynamic properties of an elastomer are a function of the strain, temperature, 
frequency, and preload to which the elastomer is subjected. In addition, 
many elastomer materials are sensitive to environmental factors, such as 
exposure to oil or other fluids, exposure to extreme temperatures, and aging. 
It is important for the designer to recognize that there are inherent varia- 
tions in elastomer materials of the same type. In particular, there may be 
considerable variation among the same generic elastomer materials produced 
by different manufacturers. Even for a particular manufacturer, the batch- 
to-batch variation in elastomer material may be significant. These varia- 
tions can have a significant effect on the dynamic properties of the 
elastomer . 

In general, support flexibility and damping are provided in vibratory systems 
to reduce the transmission of shock and vibration. Of primary interest in 
this handbook is the reduction of vibrations through the addition of damping. 
The ultimate aim of vibration reduction is to reduce the potential of 
machinery failure and safety problems and to reduce the maintenance require- 
ments of machinery. These goals are applicable to rotor vibration as well 
as structural vibration. 


* Numbers in brackets refer to references in Section 1.1 

® Registered Trademark, Shore Instruments and Manufacturing Co., Inc., 
Jamaica, N.Y. 
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The reduction of the transmission of shock or vibration may be referred to 
as either isolation or insulation. The term "isolation” is used when the 
object is to minimize the transmission of vibration from a machine to its 
environment. The term "insulation" is used when the object is to minimize 
the vibration transmitted from the environment to a machine. When referring 
to vibration isolation or insulation, as distinguished from vibration re- 
duction, the primary concern is with the transmission media. That is, the 
object is not the dissipation of vibration energy but rather the reduction 
of vibration transmission. The purpose of vibration isolation or insulation 
is to reduce noise, general discomfort, machinery fatigue problems, assoc- 
iated safety problems and maintenance, and ultimately, to reduce operating 
costs . 

The purpose of this handbook is to facilitate the use of elastomer dampers 
for applications on rotating machinery. To this end, a brief introductory 
discussion of elastomer dampers intended for those readers who are unfamiliar 
with elastomer dampers is presented. Currently, elastomer dampers enjoy 
broad use for control of unidirectional vibrations (jL . 3 ' . This is parti- 
cularly true in the areas of vibration isolation and insulation. However, 
elastomer dampers have, unfortunately, seen very limited use in rotating 
machinery. As the necessary design information and new elastomers become 
more readily available, the use of elastomer dampers in rotating machinery 
should become much more common. 

Elastomer dampers have a number of advantages over other types of dampers. 
They are easy to design, manufacture, assemble and maintain and have the 
general attribute of simplicity. They are also relatively inexpensive due 
to the requirement for less stringent tolerances and the inherent simplicity. 
Elastomer dampers are very durable and can easily handle shocks, rough 
handling, and an occasional moderate overload. If designed properly, main- 
tenance is very simple and inexpensive. Elastomer dampers can tolerate sub- 
stantial misalignment, which is of particular importance in rotor bearing 
applications. They also provide both support stiffness and damping and thus 
do not require parallel stiffness elements as conventional hydraulic mounts 
(squeeze-film dampers) . 

Elastomer dampers also have certain disadvantages which should be carefully 
considered during the design process. The inherent property of elastomer 
dampers of providing relatively low damping at high frequency is not always 
desirable. This is particularly true in rotating machines which operate at 
high frequencies. However, even when this condition exists, the level of 
damping provided by an elastomer damper may be sufficient to control rotor 
vibrations. Also, elastomer dampers tend to be somewhat sensitive to changes 
in ambient temperature (some materials more than others) and are generally 
useful only in a limited temperature range. Elastomers are also sensitive 
to other environmental factors, such as fluid compatibility. Environmental 
considerations can generally be handled successfully by an appropriate choice 
of elastomer material. One of the most formidable obstacles to designing 
consistent elastomer dampers is the variation that is generally found in the 
material composition of a single generic compound. This variation is not 
only observed between manufacturers, but also between batches from a single 
manufacturer. The designer should be aware that manufacturers may change 
the material composition of an elastomer between batches for cost or avail- 
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ability reasons* Such variations in the material composition can affect 
the dynamic properties « 

It is useful to make several observations concerning the comparison of 
elastomer dampers with more conventional dampers (particularly squeeze-film 
dampers) for rotating machinery. Using proper design data, the predictions 
of stiffness and damping for elastomer dampers are just as good, or better, 
than those for squeeze-film dampers* Elastomer dampers are much less ex- 
pensive and easier to manufacture and maintain than squeeze-film dampers. 
Elastomer dampers also provide much more adjustment flexibility during 
assembly and test procedures than do squeeze-film dampers. In particular, 
an elastomer damper can be designed in a modular fashion so that the elasto- 
mer material or configuration can be easily modified in order to change the 
damper properties. This is especially useful for prototype machines. 

This handbook provides dynamic properties of several types of elastomers 
particularly suitable for use in rotating machinery dampers. (These dynamic 
properties have been determined for a range of operating conditions which 
span most practical rotordynamic applications.) The specific steps required 
for the design of elastomer dampers are reviewed in detail and include the 
procedures and justifications for the selection of a particular elastomer 
material and a particular elastomer geometric configuration, as well as the 
details involved in the design of the actual damper' hardware. In order to 
provide as complete a reference as possible, a full range of physical proper- 
ties, such as density and thermal conductivity, are also provided for a 
large variety of elastomer compositions. A discussion of environmental con- 
siderations such as fluid, ozone, and adhesive compatibility is also pre- 
sented along with specific compatibility data for several elastomer materials. 
For the designer who finds it necessary to determine the dynamic properties 
for a particular sample of elastomer material, or for an elastomer material 
for which the dynamic properties are not available, a thorough discussion of 
the most common dynamic property test methods is presented. This handbook 
includes a detailed description of the procedure for performing the base 
excitation resonant mass test, which, of the tests discussed, is the most 
versatile and least sensitive to instrumentation error. Specific detailed 
examples of existing elastomer dampers are discussed, along with the steps 
encountered in the design process, to provide the reader with a complete 
picture of the design procedure. The design engineer is encouraged to use 
these elastomer damper designs as a basis for his own damper design. 

1.1 References 

1.1 Payne, A, R. , and Scott, J. R. , Engineering Design with Rubber, 
Interscience Publishers, Inc., New York, 1960. 

1.2 Elastomer Material , Desk-Top Data Book, The International Plastics 
Selector, Inc., San Diego, California, 1977. 

1.3 Gobel, E, F., Rubber Springs Design , John Wiley & Sons, New York, 

1974. 
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2,0 FUNDAMENTALS OF DAMPER DESIGN 


An understanding of basic vibration theory is required for the proper design 
and specification of machinery dampers. For reader convenience , a brief 
introduction to vibration theory is presented in this chapter. To fully 
appreciate rotor vibration, an understanding of linear vibration theory is 
fundamental. Thus, the introduction includes discussions of both unidirec- 
tional and rotor-bearing vibration theory. Of course the principal interest 
here is rotor-bearing system vibrations, A theoretical and historical back- 
ground is also presented concerning the use of linear damping models, with 
particular emphasis on rotor-bearing systems. Included is a discussion of 
currently used types of rotor dampers and a comparison with elastomer dampers. 

2,1 Unidirectional Vibration Theory 

This section is intended to provide no more than a brief introduction to 
vibration theory. Many references, such as [ 2 . lj * through [2.3], are avail- 
able which are much more complete in this regard. 

Most real vibratory systems are composed of a continuum with finite mass 
energy-storage and energy-dissipation characteristics. An ideal vibratory 
system may be considered to be composed of a number of individual elements, 
each with the capability to either store potential or kinetic energy or to 
dissipate energy. The vibration phenomenon is a manifestation of an alter- 
nating transfer of kinetic and potential energy. Simple system representa- 
tions use elements which possess only one of these capabilities. These 
representations are referred to as lumped parameter models. The number of 
these elements required to accurately represent a real vibratory system 
depends on several factors, including the geometry system constraints, trac- 
tion loads, and the anticipated excitation frequencies. For the purpose of 
this work, these lumped parameter model elements are assumed to behave 
linearly. This assumption is a reasonable approximation for most small vibra- 
tory motions. 

The three fundamental lumped parameter model elements are: the spring, for 

storing potential energy; the mass, for storing kinetic energy; and the 
damper, for dissipating energy. 

The simplest possible vibratory system representation is a single-degree-of- 
freedom model for free vibration with no damping as illustrated schematically 
in Figure 1 where m represents mass and k represents stiffness (or spring 
constant) . The term "free vibration" indicates that the system is provided 
with a set of initial conditions and then allowed to vibrate freely with no 
external influences. For the model in Figure 1, a mass is attached to a 
rigid support by a linear spring which exerts a force proportional to its 
change in length. This ideal (lumped parameter) spring is assumed massless** 
so that the forces at the ends of the spring are equal in magnitude and 


* Numbers in brackets indicate references found in Section 2.6. 

** See for details on the effect of spring mass on vibratory motion. 
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opposite in direction,, The force exerted by the mass is proprotional to 
acceleration . The motion of the mass is unidirectional with displacement x, 
where x is a function of time. The differential governing equation for this 
system may be obtained through the use of Newton’s secogd law of motion, 
which may be stated as: the rate of change of momentum -r— (mx)* is equal to 
the vector sum of forces exerted on the body and takes place in the direc- 
tion in which the force (vector sum) acts. If the mass is constant, the 
rate of change of momentum is equal to the mass times its acceleration (mx) . 

Since the system being considered is unidirectional, this law reduces to 

mx = £ (forces in the x-direction) - -kx (2-1) 

The negative sign indicates that the spring force is a restoring force 
(resists applied motion) . The form of the restoring force indicates that 
the static equilibrium position is assumed to be at x equals zero. 

Rearrangement of equation (2-1) results in the more familiar form of the 
governing equation given by 

mx + kx = 0 (2-2) 

The solution to this differential equation is 

x = Asinrn t + Bcosm t (2-3) 

n n 

where the constants A and B are determined from the initial conditions of 
the system and where w , known as the undamped natural frequency or resonance, 
is given by: 

a) -^k/m (2-4) 

The units of w are radians per second, 
n 

Equation (2-3) may be rewritten in the form 

X = Csin (a) t + 0) (2-5) 

n 

where C is the amplitude of vibration given by: 

C ="^A 2 + B 2 (2-6) 

and 0 is the phase angle of vibration given by 

0 = arctan (B/A) (2-7) 

The term "phase angle" can have slightly different meanings when used in 
different contexts. This form of vibration is referred to as simple harmonic 
motion due to the harmonic (cyclic) form of equations (2-3) and (2-5). 


* Dots over the x represent derivatives with respect to time, i.e., x = — 

dt 
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The vibratory system of Figure 1 is made slightly more complicated by the 
addition of a viscous damping element (the most commonly used model for an 
energy dissipation element) as shown in Figure 2 where b represents the 
viscous damping coefficient. The force exerted by a viscous damper is pro- 
portional to the relative velocity between the ends of the damper. The 
ideal damper * like the ideal spring, is assumed massless. Again considering 
only free vibration, the differential governing equation can be derived from 
Newton's second law to give 


rax + bx + kx = 0 (2-8) 

The form of the solution to this differential equation depends on the value 
of b relative to the values of m and k (assuming all are positive) . The 
criteria for the selection of the appropriate solution are dependent upon 
the value of the damping coefficient, b, relative to that of the critical 
damping value, b^, which is defined by 

b = 2 Vkin* = 2 mio (2-9) 

c n 

where oi is defined in equation (2-4). The ratio of the damping coefficient 
to the critical damping is referred to as the damping ratio, £, such that 

? = ( 2 - 10 ) 
c 

m For the case in which the damping is less than critical (or 

underdamped) (C < 1) , the solution to equation (2-8) is of the 
form 

x = e ^ )t ^ Tn (Asinw^t 4- Bcosw^t) = Ce s i n ( u t + 9) (2-11) 

where A, B, C and 9 are constants which depend on the initial con- 
ditions, and C and 9 are related to A and B by equations (2-6) and 
(2-7) . is called the damped natural frequency and is given by 

U) d =Vk/m (1 - C 2 ) 1/2 = 0) n (1 - C 2 ) 1/2 (2-12) 

and is given in units of radians per second. Note that for unidi- 
rectional vibration, the damped natural frequency, is less than 

the undamped natural frequency, ui . The vibratory response repre- 
sented by equation (2-11) is referred to as a transient damped 
oscillation. A plot of a typical damped oscillation response is 
presented in Figure 3. 

% For the case in which the damping is exactly critical (or critically 
damped) (C=l) , the solution to equation (2-8) is of the form 

x = (A - Bt)e ^ >t ^ m = (a 4* Bt)e W n t (2-13) 

where A and B are constants which depend on the initial conditions. 
For this case, a typical response plot is presented in Figure 4. 

e For the case in which the damping is greater than critical (or 

overdamped) (t, > 1), the solution to equation (2-8) is of the form 
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Fig. 2 Single-Degree-of-Freedom Free Vibrational Model With Viscous Damping 
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Time (t) 


Fig. 4 Typical Response for Slngle-Degree-of-Freedom Free Vibration System 
With Viscous Damping; p~ 1 (Critically Damped) 
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+ Be 


(2-14) 


Kf . , U VC 2 -1 t 

e bt/2m (Ae n 


-m v? 2 -i t 

n Y 


where A and B, as before, are constants which depend on the initial 
conditions. For this case, a typical response plot is presented in 
Figure 5, 


Most real vibratory systems are underdamped (£ < 1). This is generally true 
even with the addition of external damping supplied by specially configured 
dampers. On occasion, the addition of external damping may convert an un- 
derdamped system into a critically damped or over damped system. However, 
this is a special case. Consequently, it is necessary to have methods for 
empirically evaluating damper effectiveness for underdamped systems, in 
order to determine the need for damping and to compare the effectiveness of 
various damper designs. 

Several measures of damping effectiveness are in common use. Each of these 
measures has its own physical significance and is interrelated. One of 
these measures is damping ratio, £, which has been discussed above. Another 
is logarithmic decrement, 6, which is defined as the natural logarithm of 
the ratio of successive cycles of damped free vibration. That is, 

x 

5 = i ln( r° (2_15) 

n 

where is the amplitude of the first of the cycles, x is the amplitude 
of the Nth, and ^ is the nu mber of consecutive cycles, ^sing equation 
(2-11) , logarithmic decrement is found to be related to damping ratio by 

6 = 2 * C /(1 - C 2 ) 1/2 ( 2 - 16 ) 


For small values of the damping ratio, this relationship reduces to 

6 ^ 2tiX for t« 1 (2-17) 


Positive values of log decrement indicate stability since x 2 <x ^ and the 
amplitude decreases with time. Conversely, negative values of log decrement 
indicate instability. A log decrement value of zero indicates a limit cycle, 
where the amplitude of vibration remains constant with time. For most prac- 
tical dynamic systems, the attainment of a non-negative log decrement is not 
a sufficient stability criterion for two reasons. First, to insure safe and 
efficient operation of most machines, vibrations must not only be damped out, 
but damped out rapidly, as briefly sustained vibrations may be troublesome. 
Second, the log decrement is usually determined analytically for predicting 
system behavior. A safety margin is generally advisable whenever an analy- 
tical model is used to represent a real dynamic system. 

Loss factor or loss coef fJLcient (n) , which represents the ratio of damping 
energy to strain energy Jl2.l) , is also used as a measure of damping effec- 
tiveness. For light damping (£<0.1), the loss factor for a homogeneous 
element (not a system) is approximately equal to twice the damping ratio. 

A classification of damping units i 2.5~1 is presented in Table 1, and the 
interrelationships [^2.1J are presented" in Table 2. When applied to elasto- 
mer elements or damper assemblies, the loss coefficient is defined as 
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'-Time (t) 


Fig. 5 Typical Response for Single-Degree-of-Freedom Free Vibration System 
With Viscous Damping; p> 1 (Over Damped) 
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Response oh Systems under Excitaiion | g I a j Properties op Materials ano Specimens 


Table 1 


DAMPING CLASSIFICATION* 



(!) 

Classification of type 
of damping units 

(2) 

Preferred names 

( 3 ) 

Other names used 

(4) 

Symbols and units 
Material j Members 

< 5 > 

Definition 


Relation between specimen properties 
and system response 





1 

1 

s 

c 

XT, 

c 

z 

< 

a 

< 

1 

s 

0 
2 
P 
s 

1 

(A) Absolute energy UNrrs 

Stored or dissipated per cycle of stress. 

Material properties . 

V and D give unit values. 

Specimen properties: 

U t and 0, give total values for 
specimen S. 

Unil elastic 
Strain energy 

Elastic energy 
strain energy' 

U 

in-lb'in’-cycle 


Area under cr m ,„~F 
curve from 0 to r mi . 


D s k'f 1 1 

f? ’ ‘ ZnU, ’ *; “ Q, “ Ar “ Bm ” B * 

a A 2vt li 

7i ~ Tie) ” 5.35 to 

bet T b/.v y MV, 

77 71 “ 27.3 

Total elastic 
strain energy 



V , = uv t p 

in-lb/cycle 

Area under P mU -X 
curve from 0 to P m „ 

Unit damping 
energy 

Specific damping 
Spec, hysteresis 

D 

in-lb in ’-cycle 


Area within a-t 
hysteresis loop 

Total damping 
energy 



D, = DV t a 
in-lb'cycle 

Area within P-X 
hysteresis loop 

(B) Complex modulus notation 
Appropriate for linear materials. 

Material properties: 

E* - £' + /£" 

Same for G*, K \ A/*. 

Properties of Specimen , 
member, or total part S. 
kf ■= k' t ~ ik” 

Complex 

modulus 


£*, G\ K*. A/* 
lb/in* 

k: 

lb/in 

tn 
II II 

*l*tn 

Absolute 

modulus 


1^*1. M 
l*T- Kl 

i*?i 

|£*| =■ (£' J 4 £" J 1‘ 12 


Note: 

For linear 
materials 
«'/?- 1 
V. 1 V 

Q t -Q 

Storage 

modulus 

Elastic modulus 
Real modulus 

E\ C, K\ M* 
?b/in a 

K 

Tb/in 

E' - j£*| cos A 



Relations among unit 
properties of uniform 
materials 


Loss 

modulus 

Dissipation 

modulus 

£", G ", K", M" 
lb/in 2 

5 ‘ ' 

E" - |£*| sin f> 

'•S' 

(C) Relative energy uni re 
Dimensionless ratios 
of damping energy D 
and strain energy V. 

Loss 

coefficient 

Loss factor 
Damping factor 

r) (eta) 

’ 

« 

D 

2W 

D t 

h% " Z-xU , 

D 

* r ~ ZnV 

E" 
~ ~E r 

Quality 

factor 

Storage coef. 
Q factor 

Q 

Q„ Q l 

_ tiV 
Q -~ 

ZtU, 

* ~dT 

„ ZnU 
Q ~ D 

iO 

H 

^1 *! 

i| 

J |- 


Proper i jes of Composites, Structures. ire. These properties are not covered in this monograph. Nomenclature and units under study by a ASA Committee initiated in 1965. 


(D) T emporal decay of free vibrations jn members 

1 

Measured quantities or prop, 
of systems 

How measured 


Relations between system properties and specimens properties 


Logarithmic 

decrement 

Decrement 

J, (Delta) 

*.(»+ 1 ) = XJn)e-*t 
A -In W 


*--£r 

A -4 

J, ^ jit;. 

4, 

' ( t ) [ X„{t) \ 



* XJfl + 1) 





. (db) 20 , Y.{0) 20/ , t,(0) 


Temporal decay 

v, (Upsilon) 

r 1 In W 


wD, 

o) k'f 

or 


* 1 t f °® Y,(/) n .r„(n) 

Decay 

constant 

(in./in.) per sec 

' 1 x a (t) 


' ' ' 4-Tt/j 

” ,= 2 K 

h ^ 2 *6 

** la 

20 

-y log e v t 8 .68 n. 

constant 

Decibel decay 
rate with time 

T r (UrsrLON) 
(db) per sec 

T, = 8.68 i' t 


Y, = 0.69tu -jj- 

k‘: 

Y, - 4.35 to -~ 
k» 

Ti *• 4.35 tot;, 

4.35 oj 

t '"5T 

(E) Spatial a n enuat ion of waves in slender rods, beams, etc., having distributed parameters 







XJj) - A a (0)f-V 

(l\ X.10) 

X.{y) 

Logarithmic 

attenuation 

Attenuation 


X t ().)-= XJ0)e-*r 
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A/6 C — fraction of critical damping for viscous or equivalents viscous damping, where 
c is the damping constant of the dashpot acting on the mass m in parallel to the 
spring stiffness k and c c is the critical damping constant as given by Eq. (2.12) 

7] = loss factor of spring with complex stiffness given by &(1 + j-q) 

k — decay constant of system in free vibration given in terms of viscous damping 
by k — c/2m 

A — logarithmic decrement of free vibration of system defined in ,terms of sue- 
cessive peak displacements during decay 

D = decay rate of free vibrations expressed in decibels per second during exponential 
decay of amplitude 

A f = bandwidth of steady-state resonance of system as related at the half -power 
points 

Q — quality factor of steady-state resonance of system as- related to the resonant 
amplification factor, Ar 

$ — specific damping capacity of system in terms of ratio between energy dissipated 
per cycle and elastic energy-storage capacity of the system 


* Reference 2.30 

With permission of McGraw Hill Book Co. 
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n = 


(2-18) 


where k and are the real and imaginary parts of the complex stiffness, 
k*. 


k* - k x + ik 2 (2-19) 

This is discussed in more detail in the next chapter. In terms of a viscous 
damping coefficient, the loss coefficient may be represented by 

n = wb/k (2-20) 

Damping in elastomers is generally represented by hysteretic damping which 
is equivalent to the imaginary component of the complex stiffness given in 
equation (2-19) . Hysteretic damping results in a dissipation force pro- 
portional to displacement (as compared to viscous damping which is propor- 
tional to velocity), but ninety degrees out of phase with displacement. 

The review of forced vibration begins with a single-degree~of-freedom, 
undamped system, as illustrated in Figure 6. The forcing function is 
assumed to be sinusoidal and may be realized as a force applied to the mass, 
shown in Figure 6(a). The response of the system is defined by the force or 
motion transmissibility , depending on the form of the forcing function. The 
force transmissibility is defined as the ratio of the force transmitted to 
the foundation to the force applied to the mass. The motion transmissibility 
is defined as the ratio of the motion of the mass to the applied motion of 
the foundation IjZ.lJ. 

For the case in which a sinusoidal forcing function is applied to the mass, 
the governing equation becomes: 

mx + kx = F sinmt (2-21) 

o 

where F sinwt is the forcing function. Equation (2-21) is not homogeneous; 
thus it has both a general and particular solution which together are given 
by 

F /k 

x - Asinrn -f Bcosm ^ ^ ~ sincot (2-22) 

nt nt - l , Z 

1-03 /U) 

n 

and where A and B are constants that depend on the initial conditions. The 
first two terms on the right side of equation (2-22) form the general solu- 
tion to equation (2-21) . This is identical to the solution for the analogous 
free vibration case given by equation (2-3) . In virtually all practical 
systems, there is some damping present which causes the general solution 
terms (also known as the transient solution) to disappear, or to become 
negligible, with time. However, the particular, or steady-state, solution, 
given by the last term in equation (2-22) , does not diminish with time and 
is, therefore, of primary interest here. 

Considering only the steady-state solution, the force transmissibility, T^, 
is given by 
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Fig. 6 Undamped Single-Degree-of-Freedom Forced Vibration Models 
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(2—23) 



03 

n 


where T^, is unbounded when the forcing frequency, to, equals oj , the natural 
frequency. A plot of as a function of the frequency ratio, 03/03 , is 
presented in Figure 7. 


For the case in which the foundation has a prescribed sinusoidal motion, the 
governing equation becomes 


mx 4- kx = kU sino3t (2-24) 

o 

where (U sino3t) is the prescribed motion of the support. Equation (2-24) is 
very similar to equation (2-21) and, thus, has a similar solution given by 

U 

X - AsinoJ t + Bcoso) t H ~ tt sino3t (2-25) 

n n „ 2 . 2 

1-03 / 03 

n 

As for equation (2-22) , the first two terms on the right side of equation 
(2-25) form the transient solution, while the last term represents the 
steady-state solution. 

Again considering only the steady-state solution, the motion transmissibility , 
T m , is given by 

T m = (l-u 2 /^ 2 ) -1 (2-26) 

where T is unbounded for a frequency ratio of unity. Note that, for this 
model, the force and motion transmissibilities are identical. Thus, the plot 
in Figure 7 applies to motion transmissibility as well as to force transmissi- 
bility. 

A viscous damping element is now added to the forced vibration models of 
Figure 6, resulting in a single-degree-of-freedom system with viscous damp- 
ing as illustrated by the models of Figure 8. The forcing function, as 
before, can be represented as either a force applied to the mass (Figure 8a) 
or a prescribed motion of the foundation (Figure 8b) . 

For the case in which a sinusoidal forcing function is applied to the mass, 
the governing equation becomes 

mx 4- bx + kx - F sinwt (2-27) 

o 

where F sinwt is the applied force. For reasons given above, only the steady- 
state solution of equation (2-27) is considered. For positive values of b 
(in general, b is not negative, although it may be zero), a phase angle exists 
between the forcing function, F^sinmt, and the steady-state response x. 

This response is given by 


x = Rsin (o3t-0) 


(2-28) 
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Fig. 7 Transmissibility of a Single-Degree-of-Freedom 
Forced Vibration System With No Damping as a 
Function of the Frequency Ratio 
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where 


R = ^ [(l- W 2 /o) n 2 ) 2 + (25»/« n )^) 1/2 

and 


(2-29) 


9 = arctan [(2cca/aJ n )/(l-*a) 2 /a) n 2 )J (2-30) 

C and w are as defined in equations (2-10) and (2-4), respectively. 

The relationship between the force transmitted to the foundation and the 
applied force is more complex than for the undamped system presented above. 
For the damped system, both the spring and damper transmit force to the 
foundation. The total force transmitted to the foundation, F , is given by 

F t - bx + kx (2-31) 

where the two terms on the right are ninety degrees out of phase with each 
other. The ratio of the transmitted force, F , to the magnitude of the 
applied force, F , becomes* 


where 


= T sin (cot - ¥) 


T =A 


1 + (2£co/co ) 
n 


(l-co 2 /u) 2 ) 2 + (2 £co/w ) : 
n n 


(2-32) 


(2-33) 


¥ = Arctan 


! 2c(w/oo ) 3 

I n 

/ 2 / 2 n 2 2 . 2 n 

i 1 — (a) / a> ) + (44 co /w ) 

— n n 


(2-34) 


and T and ¥ are referred to as the transmissibility and the phase angle, 
respectively. Note that, for positive values of 4, T is bounded for all 
values of co. Thus, the presence of damping is necessary to control vibra- 
tion at a system resonance when comparing equation (2-33) to equation (2-23). 
Plots of transmissibility, T, and phase angle, ¥, are presented as functions 
of the frequency ratio for several values of £ in Figures 9 and 10, respec- 
tively. 


The maximum value of the displacement , x. from equation (2-28), occurs for a 
forcing frequency equal to (1-2£ Z ) 1/2 1. Note that, for positive the 

maximum steady-state response (in terms or displacement) always occurs for 
forcing frequencies less than the undamped natural frequency. This is 
generally true for unidirectional vibratory systems. As will be seen in the 
next section, this is not true for lateral vibration of rotating machines. 


* Equations (2-32) through (2-34) are from Reference 2.1 a 
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Fig. 10 Phase Angle versus Frequency Ratio 
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For the case in which the foundation has a prescribed sinusoidal motion 
(Figure 8b) , the governing equation becomes 

mx + bx 4- kx = bu wcoswt + ku s incut (2-35) 

o o 

where u sincut is the prescribed displacement of the foundation and, thus, 
U Q a)cos(u? is the related prescribed velocity of the foundation (the velocity 
is the time derivative of the displacement) . Again, considering only the 
steady-state solution to equation (2-35), the response is given by 

x = Tu q sin (cut - ¥) (2-36) 

where T and ¥ are defined by equations (2-33) and (2-34), respectively. 
Analogous to equation (2-32) , the motion transmissibility and phase angle 
are defined by 


— = T sin (ut - Y) (2-37) 

O 

from equation (2-36). Thus, the force transmissibility and phase angle are 
identical to the motion transmissibility and phase angle for the correspond- 
ing forms of the forced vibration. Consequently, Figures 9 and 10 apply to 
both damped forced vibration cases. 

In general, practical vibratory systems have more than one degree of freedom. 
Practical systems can usually be represented with reasonable accuracy by a 
finite number of spring, mass, and damper elements, connected together in 
some fashion. This is often referred to as a lumped parameter representation. 
Real systems may exhibit free or forced vibration or a combination of both. 
Damping is often added to systems intentionally to aid in vibration reduction. 
In such cases, one end of the damper is generally attached to a rigid support. 
Reduced support stiffness (by the installation of soft members) is sometimes 
used to isolate (or insulate) vibration. Dynamic absorbers, which them- 
selves have two or more degrees of freedom, are also used occassionally to 
isolate (or insulate) vibration. For a more complete discussion of multiple- 
degree-of-freedom systems, the reader is referred to Reference 2.1 thropgh 
2.4. 


2.2 Vibrations in Rotating Machinery 

Section 2.2 is intended to provide an introduction to the theory of rotor 
bearing vibrations. A few references, [2.6j through [2.9], are available which 
provide a fairly complete treatment of rotor bearing vibration. 

Any machines which have rotating components are classified as rotating ma- 
chines. This classification includes machinery in a broad range of sizes. 

Some typical small rotating machines are gyroscopes, electric motors for 
household appliances, and turbochargers. Medium-size rotating machines 
include aircraft gas turbine engines, small industrial compressors, electric 
motors, and internal combustion engines. Some large rotating machines are 
steam turbine /generator sets, ground-based gas turbine/generator sets, large 
industrial compressors, large diesel engines, and power transmission equip- 
ment. Vibrations of rotating machines can lead to serious problems, parti- 
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ularly in the areas of noise pollution, operator (and consumer) safety, 
and machinery operating expenses. For example, unscheduled downtime and 
repair costs resulting from unexpected machinery failure can have a devas- 
tating impact on operating expenses. Consequently, areas which have seen 
intensive developmental effort in the past several years (with a high level 
of continued effort in the foreseeable future) are vibration monitoring, 
vibration diagnostics, and vibration suppression; the latter includes the 
development of elastomer dampers for rotating machinery. Suppression of 
vibration in rotating machinery is generally accomplished either by re- 
moving the sources of vibration or by the dissipation of vibrational energy 
through damping (or, most frequently, by a combination of the two). 

Rotor vibration is generally classified according to the direction and 
frequency of the vibratory motion. The three directions of motion illus- 
trated in Figure 11 are lateral, axial, and torsional. Lateral motion 
involves a displacement and/or rotation of the rotor centerline. In this 
context, rotation is about an axis intersecting and normal to the rotor 
centerline. Axial motion occurs parallel to the rotor centerline. Tor- 
sional vibration involves a rotation (or twisting) of the rotor about its 
centerline relative to a frame of reference fixed at some point in the 
rotor. Vibrations which occur at the frequency of rotation of the rotor 
are referred to as synchronous vibrations. Those which occur at other 
frequencies are known as nonsynchronous vibrations. In geared systems, 
where two or more rotors are operating at different speeds, a different 
synchronous frequency is associated with each rotor. 

Elastomer dampers have been used to control torsional as well as lateral 
vibration. Further discussion of rotor vibration in this handbook will 
address only lateral vibrations. Both synchronous and nonsynchronous vibra- 
tions can be controlled by elastomer dampers. Synchronous rotor (lateral) 
vibration is usually due to mass unbalance in the rotor, although other 
synchronous forcing functions are occasionally present. Rotor mass un- 
balance results from an eccentricity of the centroidal axis of the rotor 
relative to its geometric center. In practice, this unbalance is 
generally a result of unavoidable imperfections in rotor manufacture and 
assembly. Synchronous vibration is usually stable; although, if large 
motions are sustained, it can often lead to rapid machinery failure. 

Synchronous vibration due to unbalance is most often controlled through a 
combination of rotor balancing and external (or support) damping, particu- 
larly when the rotor operates above one or more critical speeds* (or synchro- 
nous resonances) . Rotor balancing refers to any of a number of procedures 
used for applying intentional unbalance to a rotor in order to compensate 
for existing unbalance and thereby reduce synchronous rotor vibration 
{~2.6, 2.9j. Since a rotor can never be n perf ectly" balanced, some damping 
is required to control vibrations at rotor critical speeds (although 
bearing damping is sometimes sufficient) . Very often, rotor balancing 
by itself is not sufficient to adequately control synchronous rotor 


* which operate at speeds above critical speeds are often 

termed "high speed" rotors 
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vibration* Thus, the use of external damping is often required in addition 
to rotor balancing* Also, the addition of external damping generally per- 
mits the reduction of balancing tolerances, which can result in a substan- 
tial reduction in rotor manufacturing and assembly costs* 

In 1919, H. H. Jeffcott ( 2 .I 0 J presented an analysis of a simple, concen- 
trated-mass, flexible rotor which included the effect of external damping* 
Although the model is very simple, the analysis demonstrated the shifting 
of phase of the center of gravity of the rotor as a critical speed is 
traversed. The analysis has been proved experimentally and has subsequently 
formed the basis of later analyses of synchronous vibration in damped rotors, 
using more realistic rotor models* Jeffcott f s analysis is based on the rotor 
model illustrated in Figure 12. This model consists of a longitudinally and 
laterally symmetric rotor with a single concentrated mass in the center 
supported by a flexible, massless shaft, which is, in turn, mounted on 
simple supports. This analysis also applies to the rotor model presented in 
Figure 13, as long as the values of the model parameters (mass, stiffness, 
damping and eccentricity) are the same for both models. The analysis assumes 
small deflections, linear and symmetric stiffness, and viscous damping with 
no gyroscopic effect. An end view of the rotor model is presented in 
Figure 14, where the origin of the stationary (nonrotating) coordinate sys- 
tem "XY" is located at the undeflected centerline of the rotor. For a 
vertical rotor, the undeflected centerline is coincident with the line of 
bearing centers. For a horizontal rotor, these two lines differ by the static 
gravity sag of the rotor. The point E represents the deflected centerline of 
the rotor, while the point M represents the location of the center of mass of 
the rotor. Since the coordinate system is not rotating, the sketch in 
Figure 14 represents an instant of time when the line EM is at an angle cot 
from the X axis. Rotating of the rotor is counter-clockwise. The coordinates 
of the point E at that instant in time are x and y. The coordinates of the 
point M at the same instant in time are x" and y", related to x and y by 

x" = x + a coscot (2-38) 

y" = y + a sincot (2-39) 

where a is the fixed distance between points E and M. This fixed distance 
is known as the mass eccentricity of the rotor. The differential equations 
of motion for the rotor are 


mx" + bx + kx = 0 (2-40) 

my" + by + ky = 0 (2-41) 

where m is the mass of the rotor and b and k are the external viscous damping 
and stiffness coefficients, respectively. In the X-direction, the inertial 
term, mx", is proportional to the acceleration in the X-direction of the point 
M. However, the damping and stiffness act at the point E and, therefore, 
result in forces proportional to x and x, respectively. The forces in the Y- 
direction act in a similar fashion. Substituting equations (2-38) and (2-39) 
into equations (2-40) and (2-41), eliminating x" and y", and rearranging terms 
gives equations (2-42) and (2-43) , which are valid for a constant value of 
rotational speed, m. 


27 




28 




29 





mx + bx 4* kx = maw coswt 


(2-42) 


my 4- by 4- ky = maw sinwt (2-43) 

Note that the rotor system is assumed to be symmetric with respect to the XY 
coordinate system. That is, the system parameters (damping and stiffness) 
are assumed to have equal values in the X and Y directions. (Further, gyro- 
scopic or bearing cross coupling terms have not been included.) 

Equations (2-42) and (2-43) may be easily solved by substituting an assumed 
solution into the differential equations and solving for the constants, 
since the form of the solution is obvious from Figure 14. Jeffcott solved 
these equations for both the transient (general) and steady-state (particular) 
solutions. Only the steady-state (synchronous) solutions are of interest 
here. From Figure 14, it is clear that 

x = Acos(wt-3) (2-44) 


and 


y = Asin(wt-3) 


(2-45) 


where A is the distance from 0 to E and 3 is known as the phase angle. 


Substituting equations (2-44) and (2-45) and their derivatives into equations 
(2-42) and (2-43) and solving for the constants A and 3 gives: 

2 


A = 


maw __ 

272' T ,2" "2 

V (k-mw ) + b w 


(2-46) 


3 


bw 

arctan ; 

k-mw' 


where 0 < 3 < tt 


(2-47) 


The critical speed of the rotor is defined as the speed at which the value of 
A is a maximum. The equation for the critical speed, w , may be derived by 
differentiating equation (2-46) with respect to w, setting it equal to zero, 
and solving for w, resulting in 


2k 

CR 'V4mk-2b 2 


(2-48) 


It is apparent from equation (2-48) that the effect of increasing the damping 
is to increase the critical speed. 

The phase, 3, increases from zero to 180 degrees (0 to tt radians) as the speed 
is increased and has a value greater than or equal to 90 degrees (0.5 tt 
radians) at the critical speed. By comparison, the phase angle at resonance 
for a unidirectional vibratory system is less than or equal to 90 degrees. 
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The amplitude of rotor precession for any rotational speed, given by equation 
(2-46), is directly proportional to rotor mass eccentricity, a. Thus, it is 
also proportional to the rotor mass unbalance represented by the term (ma) . 
The rotor transmissibility, T , is defined as the ratio of the amplitude, A, 
to the mass eccentricity, a, and is given by 


T 

r 


A _ maj 

a ” 2.2 , ,2 2 

V(k-mw ) + b a) 


(u/o) )' 
n 


~V7l - (co/co ) ^ + 4 1 ^ (to/to ) ^ 

^ n — ' ti 


(2-49) 


where a) and t are given by equations (2-4) and (2-10), respectively. Equa- 
tion (2-49) is clearly analogous to the equation for transmissibility of a 
unidirectional vibratory system, given by equation (2-33) . If the damping of 
a rotordynamic system is reduced to zero, a) becomes equal to and the 
amplitude becomes unbounded at the critical speed for any value of rotor mass 
unbalance. Thus, for operation at a critical speed, rotor amplitude is 
limited only by the existence of damping. 


Nonsynchronous vibrations, unlike synchronous vibrations, are often unstable. 
These unstable vibrations are generally characterized by an onset speed or 
load. Below this onset, this vibration is stable and generally not apparent. 
Above the onset, this vibration may be unstable and is often destructive. 

The addition of external damping generally raises the speed or load at which 
the onset of instability occurs. If the onset is raised beyond the operating 
conditions of the rotor, this vibration is considered to be acceptable. 

There are many possible sources of nonsynchronous vibration. Several of 
these are listed below. References are provided for the interested reader. 


1 . 

Oil whip in fluid-film bearings [j2.ll, 2.12~] 


2. 

Ball-pass excitation in rolling-element bearings 


3. 

Alford’s tip clearance forces [2.120 


4. 

Propeller flutter 02.14] 


5. 

Torque whirl 02 . 15j 


6 . 

Internal damping from material hysteresis, shrink fits, 
spline friction and coupling friction [2.16, 2.1^ 

7. 

Gear mesh excitation 


8. 

External sources of vibration such as the firing 
reciprocating machines, and floor vibrations and 
and pipe vibration from other machines 

of 

fluid 

9. 

Rotor/stator rubs 02.18J 


10. 

Anisotrophy in rotor, e.g. , nonsymmetric shaft stiffness 
(2.19, 2.20j 
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Additional information on general non synchronous rotor vibrations may be 
found in References 2.21, 2.22, and 2.23. 

2.3 External Damping in Rotating Machinery 

The use of external damping to control vibrations in rotating machinery is 
not only desirable but is often necessary to insure safe machinery operation. 
External damping, which i>s generally stabilizing, is defined as damping act- 
ing on motion relative to a stationary frame of reference. External damping 
is most often provided in bearings or bearing supports* Internal damping, 
as distinguished from external damping, is defined as damping which acts on 
motion relative to a rotating reference frame. In general, internal (shaft) 
damping has a destabilizing effect on high-speed rotor operation. 

External damping is often used in conjunction with rotor balancing for con- 
trolling synchronous vibration by reducing vibration amplitude. The use 
of external damping in this context generally permits balancing tolerances to be 
relaxed and increases the tolerance of a rotor to long-term balance degra- 
dation. The existence of external damping also increases rotor system 
tolerance to sudden changes in unbalance which may be due, for example, to 
blade loss in a gas turbine engine. The external damping helps to control 
the transient as well as the steady-state response resulting from such a 
sudden change in unbalance. Occasionally, the application of external damp- 
ing actually eliminates the need for balancing and is a very reliable method 
of controlling nonsynchronous vibrations. Thus, external damping is often 
used for reduction of vibration amplitude as well as for control of poten- 
tial instabilities from nonsynchronous vibration. In many cases the appli- 
cation of external damping is the only means of vibration control, short 
of redesigning the rotor system or changing the operating specifications 
(e.g. , reducing design speed or load). 

External damping may be provided in several ways. In some cases the external 
damping is coincidental in the design of the rotor system and not specified. 

In such cases it is very difficult or, in fact, impossible to either predict 
or control the level of the external damping. The most common sources of 
coincidental external damping are bearings and structure. Fluid-film 
bearings often provide substantial damping to the rotor system. On the 
other hand, rolling-element bearings generally provide very little external 
damping. Consequently, rolling-element bearings must often be designed to 
operate in series with some type of mechanical damper. 

Mechanical dampers can provide a source of external damping. In such cases 
the level of damping is relatively controllable and predictable depending 
upon the type of damper used. Mechanical dampers are most often designed 
to operate in series with bearings as integral parts of the bearing supports. 
Several types of mechanical dampers are in current use. Currently, the most 
commonly used mechanical damper in rotating machinery is the squeeze-film 
damper (hydraulic mount). This type of damper is, essentially, a nonrotating 
fluid-film bearing. The damper is composed of concentric nonrotating journals, 
and the space between the journals is filled with a fluid, usually oil. 
Squeeze-film dampers have a number of advantages of which the most important 
is their acceptance by designers based upon years of experience. If properly 
designed, they are capable of providing high levels of damping. In general, 
the operating temperature of squeeze-film dampers is limited only by the 
entrapped fluid. 
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There are a number of disadvantages to the use of squeeze— film dampers, 
particularly in comparison to elastomer dampers. Squeeze— film dampers 
generally require some sealing or scavenging of fluid, or both, and often 
require a filtered oil supply and/or cooling hardware. This results in an 
increase in the cost and complexity of the hardware. The higher cost of 
squeeze-film dampers is due to the complexity of design, the large number 
of parts required, and the tight tolerances that are generally necessary. 
Squeeze-film dampers may also require a parallel support stiffness component 
which represents yet more hardware and higher cost. For the interested 
reader, additional information on squeeze— film dampers is available from 
References 2.24 through 2.27. 


Another type of mechanical damper which is in relatively common use in 
rotating machinery is the coulomb damper. This type of damper provides 
a dissipative force which is independent of position or velocity and always 
acts in a direction which opposes the relative velocity across the damper. 
Coulomb damping is usually provided by dry friction between two solid sur- 
faces. Coulomb dampers have the advantages of relative simplicity and low 
cost. However, they also have a number of disadvantages which tend to 
limit their application. It is extremely difficult to predict or design 
for the level of damping in a coulomb damper. A threshold force exists 
below which the coulomb damper does not participate in the dynamics of the 
rotor system and does not provide any energy dissipation. Coulomb dampers 
often exhibit substantial inconsistency in operation. Specifically, they 
are generally sensitive to small changes in preload force, operating 
temperature, and mating surface wear, all of which tend to change with time 
during operation. In general, coulomb dampers are not considered to be 
very reliable and, thus, are generally not used for most critical rotor 
applications. They also do not lend themselves to compact designs, a 
requirement for many applications. Additional information on coulomb dampers 
is available from References 2.28 and 2.29. 

Elastomers are beginning to see use as mechanical dampers for rotating 
machinery, and are occasionally used in the form of 0-rings for sealing 
squeeze-film dampers. In this application the elastomers provide an 
additional level of damping. Dampers which employ elastomers as the sole 
source of damping are not currently in common use. This is primarily due 
to the general lack of demonstrated reliability and lack of available design 
data for elastomer dampers. 

For completeness in this discussion of damping, a brief description of visco- 
elastic damping in structural (constrained layer) applications follows, 
although it is not the main theme of this handbook, 

2.4 Viscoelastic Structural Damping 

The application of energy-dissipative material to structural elements for 
achieving a reduction in unwanted vibrations is certainly not new. A method 
often employed in the automotive industry is to spray a panel with a pro- 
tective coating which not only protects the panel from environmental attack 
but also provides a mechanism for the reduction in noise level. This under- 
coating technique is often referred to as unconstrained extensional* damping, 
since the energy loss mechanism is provided by the coating when the base 
material experiences tensile or compressive deformations (Figure 15a). The 
aircraft industry also uses this technique by applying a damping layer in 
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the form of damping tape over the base structure. Shear* damping is another 
popular form of damping treatment which is characterized by the sandwiching 
of an energy dissipative layer between two layers of more rigid materials. 
Deformation of this constrained sandwiched structure (Figure 15b) provides 
shear loading in the damping layer which may reduce vibrations. 

One of the first works in this area was published by Oberst who studied the 
bending vibrational effects of an elastic plate coated on one side with a 
layer of viscoelastic material 2.32 . Oberst found that not only is the 
damping of the composite structure dependent upon the loss factor of the 
viscoelastic material, but that damping is also dependent on the stiffness 
and thickness of the viscoelastic material. 

Another pioneer in this field was Kerwin who is generally acknowledged as 
being first to provide an analysis of the shear damping mechanism of a 
constrained viscoelastic layer (constrained between the thin elastic con- 
straining layer and the basic elastic structure). Kerwin 2.33 establishes 
a complex stiffness relationship, demonstrates that the neutral axis of the 
composite system is frequency dependent, and also demonstrates that the 
normal force in the viscoelastic layer may be neglected, since the magnitude 
of the modulus for the viscoelastic material is generally at least an order 
of magnitude smaller than the modulus of the elastic layers. 

Kerwin T s work, although limited to an infinitely long or simply supported 
beam which is subject to bending, has been extended by Di Taranto and 
Blasingame ^2.34, 2.35j to a more generalized form. Reference 2.34 develops 
a sixth-order, complex, homogeneous partial differential equation which 
describes the axial displacement of a viscoelastic layer which is constrained 
between the two elastic layers of the beam. Mead and Markus 2.36 have 
extended the work of Di Taranto et al by developing a sixth-order partial 
differential equation for the transverse vibrations of the same beam con- 
figuration. Both the work of Di Taranto and Blasingame 2.34, 2.35 and 
Mead and Markus' publication 2.36] are restricted to a free beam which is 
not subjected to axial loading! ^ 

Many other authors have extended the state of analytical and experimental 
investigation to multilayered and other variations of the basic three- 
layered configuration. For example, Derby, Ruzicka, Schubert and Pepi 
2.37, 2.38J have presented comprehensive experimental data relating to 
the viscoelastic shear damped characteristics of composite structures. 
Reference 2.37 provides a review of the design equations including the govern- 
ing d_esign parameters for various shear damped configurations. Derby et al 
-2.3 7 j also provides extensive experimental results for L-shaped, T-shaped, 
circular, tubular, and rectangular composite material cross-section beams. 

The iterative scheme for design implementation presented in (2.37J was 
improved upon by the latter publication [2 . 3Sj to provide a more direct 
utilization of design data. 

Grootenhuis, Agbasiere and Nakra [2.39,^2.40, 2.41*1, along with Frohrib 
2.42 , Nashif and Nicholas !J2.43, 2.44 , Plass 2T. 46 J and Schwarzl 2.47^ 
have published extensive literature on multiple-fayered rectangular beam 
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See Ross, Ungar and Kerwin [~2.3l] for a review of both types of damping 
treatments. 
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treatments as well as on viscoelastic composite structures. Investigations 
such as those published by Parfitt [2-.48] and Plunkett and Lee (2.49^ have 
analyzed the sandwiched viscoelastic-elastic structure for various cut 
lengths of the constraining layer. Reference 2.49 provides a thorough dis- 
cussion of thidv system as well as an optimization of the length of the con- 
straining layers for single and multiple-laywered laminated structures. An 
excellent review of the state of the art in viscoelastic structual damping 
is contained in the work of Nakra "2.50^. Recent applications are provided 
in papers presented at a conference" on viscoelastic damping [2.5lj. 

2 . 5 Summary 

An understanding of basic linear vibration theory is required for effective 
design of machinery dampers. A dynamic system may be represented as a series 
of masses, springs and dampers. Several measures of damping effectiveness, 
as a primary mechanism for vibration control, are in common use, including 
damping ratio, logarithmic decrement and loss coefficient. Dynamic systems 
may be subjected to free or forced vibration. Free vibration exists when a 
dynamic system is given a set of perturbations and then allowed to vibrate 
freely, without an external influence. When a dynamic system vibrates under 
the influence of an external forcing function, usually harmonic, the system 
is said to undergo farced vibration. Although real vibratory systems possess 
an infinite number of degrees of freedom, they can generally be represented 
reasonably well by models with a finite number of degrees of freedom. 

Rotordynamic systems may be subjected to synchronous and nonsynchronous vibra- 
tions. Synchronous vibrations have a frequency equal to the speed of rotation 
of the rotor. Nonsynchronous vibrations occur at frequencies other than 
synchronous. Synchronous vibrations are most often due to rotor mass un- 
balance. Nonsynchronous vibrations, which are frequently unstable, may 
result from any of a number of sources, several of which are listed above. 

Both synchronous and nonsynchronous vibrations, if not controlled, are 
potentially destructive. Damping is one of the most practical means a 
designer has to control synchronous as well as nonsynchronous vibrations. 

The use of external damping to control the vibration of rotating machines has 
been discussed. Reference has been made to the use of dampers in bearing 
supports, and particularly to squeeze-film and coulomb dampers which are 
currently the most popular types of mechanical dampers. Elastomer dampers 
are not currently in common use for rotating machines. However, they have 
an excellent potential in this area. Their primary disadvantage is a lack 
of available design data. This handbook seeks to fill this need by providing 
design data for expanding the use of elastomer dampers. 

The material in this chapter was intended only as introductory background 
information and was not in any way intended to be exhaustive. For more 
detail, numerous references have been cited and are listed below. 
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3.0 REVIEW OF LITERATURE ON ELASTOMERS AND ELASTOMER DAMPERS 


There is substantial literature concerning elastomer materials , most of 
which has appeared since the development of synthetic polymers in the 1940s. 

A survey of this literature is presented in this chapter with emphasis upon the 
dynamic properties and applications of elastomers. Most of this literature 
concerns theoretical and analytical work conducted in the areas of visco- 
elastic theory and elastomer chemistry. A significant quantity of experi- 
mental work with elastomers has also been reported, but much of this is in 
the form of laboratory experiments. Relatively little has been published 
describing practical applications of elastomers to vibration problems. The 
material that does exist is almost entirely in the area of unidirectional 
loading. Few publications are available concerning practical applications 
of elastomer dampers for rotating machines. 

The literature has been divided into four broad categories consisting of : 

® Theory and Elastomer Chemistry 

# Experimental Determination of Elastomer Properties 
« Material Properties 
0 Application 

Obviously, each category may contain information which spans a number of 
these classifications. Accordingly, each is categorized according to its 
main emphasis and the author T s intent. 

3.1 Theory and Elastomer Chemistry 

Three of the more well-known texts jj3.1, 3.2, 3 . jT}* present a detailed develop- 
ment of viscoelastic theory. Bland T s presentation f3, l] is highly analytical 
while those of Fltigge p3.2j and Ferry [~3.3j are oriented more toward physical 
interpretation. Ferry *~s discussion of elastomers is general, including a 
review of elastomer chemistry and the determination of dynamic properties. 
Treatments of viscoelasticity are also presented in other referenced texts. 

For example, Argonne et al, [3.4] » presented a brief treatment of visco- 
elasticity, along with a discussion of elastomer chemistry. An extensive 
bibliography (as of 1966) was also included. Lazan [3.5~] included a section 
on viscoelasticity with particular emphasis on material "clamping. Smalley and 
Tessarzik (j.6^j included a section on thermo-viscoelasticity in which they 
considered the effect of temperature on viscoelastic behavior. Additional 
references on viscoelasticity include Alfrey [*3.7], Eirick ("3.81, Stuart 
[3.9l, and Tobolsky [3 . lCf] . Some of these references include discussions 
on elastomer chemistry as well as viscoelasticity. 


* Numbers in brackets indicate references found in Section 3.6. 
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Several authors have extended the theory of viscoelasticity to consider 
numerical simulation of elastomer mechanical behavior, Herrmann [3.11j 
developed a general numerical approach, based on variational methods, for 
calculating mechanical properties of incompressible and nearly incompres- 
sible materials, including elastomers, Korites and Nelson [3.12] developed 
an analysis to consider the heating effect on the loss factor o£ an elastomer 
element mounted in a constrained layer and have also provided some experi- 
mental correlation. Ting and Tuan 13.13 investigated the temperature 
distribution in a viscoelastic cylinder resulting from cyclic internal 
pressure. 

The molecular composition and chemical structure of polymers is discussed in 
many works. Materials Handbook ;3.14; includes a cursory discussion of the 
chemical structure of a number of generic elastomer materials. Corrosion 
[_3.15j includes a discussion on polymer chains and the effect of various 
chain configurations on the "environmental properties" of elastomers. (The 
term "environmental properties" refers to the capability of an elastomer to 
resist degradation from various environmental sources such as chemical sol- 
vents and surface abrasion.) Moore and Geyer [3.16J presented a theoretical 
evaluation of the adhesion properties of elastomers. (In this context, the 
term "adhesion properties" refers to the surface friction existing between 
an elastomer material and other solid surfaces, such as glass.) Additional 
references on elastomer chemistry include Billmeyer 3.17 , Flory 3.18 , 
Noller 1 3.19}, and Schmidt and Marlies 3.20~j. 

3.2 Experimental Determination of Elastomer Properties 

The properties referred to here include both the general physical and dynamic 
properties of elastomers. In this context, general physical properties refers 
to nonmechanical properties (such as density and thermal conductivity) , as 
well as static mechanical properties. Dynamic properties refer to the 
dynamic mechanical properties, such as shear modulus and Poisson T s ratio. 

The test methods used for determination of these properties and the results of 
these tests are of significance. 

There are numerous references dealing with the determination of the general 
physical properties of elastomers (including the effects of environmental 
conditions). For example, Adkins and Gent [3.21^ investigated the static 
stiffness of natural rubber, cylindrical bushings for various directions of 
loading. Beerbower et al [3,22: investigated elastomer-fluid compatibility 
for hydraulic systems and have developed a method by which this compatibility 
can be predicted for various elastomer-fluid combinations. The American 
Society for Testing and Materials (ASTM) has published standards for many 
general physical property tests for elastomers jj3 . 2 3j . 

Numerous works provide information on elastomer properties. Substantial 
testing has been conducted to. determine the dynamic properties of various 
elastomer materials. Payne [3.24, 3.25, 3.26J described test procedures for 
measuring shape factor and determining the form of stress-strain relationships 
of elastomers. He also provided experimental and theoretical comparisons. 
Although primarily concerned with the static loading of elastomers, he also 
considered the effect of dynamic loading on the mechanical properties. 
Baltrukonis et al [3.27J described a test method for determining the dynamic 
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properties of elastomers using torsional deflection* Rightmire [3.28j des- 
cribed a method for measuring Poisson's ratio of elastomers. Rightmire 
notes that small variations in Poisson's ratio (less than 1%) can result in 
substantial differences in the load capacity of elastomers. Cannon et al 
{j3.29j described a resonant mass test method for elastomer compression 
specimens for which data could be taken at resonance only. Accordingly, 
phase angle measurement was not required for this test procedure, but the 
frequency range of data acquisition was extremely limited. 

In 1973, a symposium on the measurement of dynamic properties of elastomers 
and elastomeric mounts was presented at the Society of Automotive Engineers, 
Automotive Engineering Conference. Fourteen papers were presented at the 
symposium and subsequently published in a proceedings. Although the primary 
interest was with elastomer mounts for automobile suspensions, these papers 
provided a great deal of valuable insight into the general area of dynamic 
property measurement for elastomers. These papers cover a broad range of 
subjects including polymer development, test procedures, and product design. 

Those papers concerned with test procedures were almost unanimous in support 
of the transmitted force testing technique over the resonant beam method. 

Miller [3_. 30i presented a history of dynamic testing of elastomers whereas 
Leonard [3. 31j described Jzhe mechanical fundamentals involved in elastomer 
testing. Nachtigal [~3.32j presented the fundamentals of instrumentation 
and computation for elastomer tests. Brouwer < 3.33 j presented a review of 
the resonant beam test method, with emphasis on its inherent problems. Owens 
[3.341 and Wilson [3.35] described the transmitted force method. Hillberry 
et al^ 0.36j described an instrument used for performing transmitted force 
testing with direct analog readout. Danko and Svarovsky [3.37[J discussed the 
use of minicomputers for the determination of nonlinear damping coefficients. 
Nemec [3.38 1 and Paul and Hillberry [3.39 discussed an error analysis and 
correlation" of several existing resonant "beam test machines. Anderson |3.40j 
described a parameterization of_ elastomers to relate test results to product 
design. Sommer and Meyer ^3.41! discussed the effect of composition, pro- _ 
cessing and test conditions" on the dynamic properties measured. Puydak Ji.4^ 
and Koons [3.431 discussed dynamic testing and quality control in polymer 
development . 

In 1973, a symposium on recent advances and developments in testing rubber 
and related materials was presented by the American Society for Testing and 
Materials. Twelve papers were presented at this symposium and subsequently 
published in a proceedings 3.44 . Most of the papers in this volume were 
primarily concerned with elastomer chemistry or static testing of elastomers. 
However, two of these papers related directly to dynamic testing of elastomers. 
Warner f3.45j presented a detailed analytical background for the dynamic 
testing "of elastomers and Hillberry j*3 • 46~1 presented a review of current (1973) 
test equipment. 

At about the same time, another method for dynamic testing of elastomers, 
referred to as the Base Excitation Resonant Mass (BERM) technique, was being 
developed. In a series of test programs extending over several years^ Chiang 
and Tessarzik 3.47*1, Gupta et al _3.4e[ , and Smalley and Tessarzik [3.49 
developed the BERM "technique and generated a substantial quantity of “test 
data, primarily for polybutadiene. Subsequently, Darlow and Smalley 3.50 
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adapted this technique for use with a rotating load test machine to provide 
a correlation between unidirectional load elastomer test results and rotat- 
ing load applications. A reasonable correlation was found that verified the 
applicability of unidirectional testing to the design of dampers for rotating 
machinery. Smalley et al (^3.5lTj presented a review of several test methods 
for elastomer testing, including the BERM method, as well as a discussion of 
instrumentation, testing philosophy and a brief analytical background for 
elastomer testing. 

Many authors have presented the results of elastomer tests (both static and 
dynamic) for a variety of elastomer materials. Painter presented the 

results of a series of tests with a broad temperature range (BTR) elastomer 
using an analog, nonresonant mass test system. These results showed that 
the effect of temperature for this elastomer material was less than the 
corresponding effect for natural rubber of Buna-N. Gent and Lindley <3.53 
presented experimental results derived from static tests for various shapely 
of elastomers, along with theoretical correlations. Hirsch and Boyce [3.54J 
described the result of a series of tests from which they determined the 
dynamic properties of an ethylene/acrylic elastomer. Maza 3.55J presented 
test results which related the dynamic properties of several elastomer ma- 
terials, for frequencies above a kilohertz. Extensive test results for 
polybutadiene elastomers, using the BERM method, have been reported by Gupta 
et al Smalley and Tessarzik [3.49], Tecza et al [3.56 J, and Darlow 

and Smalley J3.57]. ' 

Results have also been reported which relate the effect of various test 
parameters to the dynamic and static properties of elastomers. The test 
parameters most often of interest are temperature, frequency, strain, preload, 
and elastomer geometry . Ferry et al 3.58]] and Williams et al jJ3 . 59] intro- 
duced the use of a single reduced variable to include the effects or both 
temperature and frequency on the dynamic properties of elastomers. This 
method of parameterization, referred to as the Williams-Landel-Ferry (WLF) 
method, permits accurate predictions to be made for operating conditions 
which are beyond the range of possible test conditions. This is done by 
varying frequency in such a way as to represent a particular variation in 
the temperature, or vice versa. Jones !3.60~J successfully used the WLF 
method for several elastomer materials to determine temperature sensitivity. 
Moghe and Neff [3.6lj presented the results of static tests in which the 
addition of filters was found to increase the stiffness and reduce the 
Poisson T s ratio of elastomers. This had a subsequent effect on the static 
shape factor. Sircar and Lamond [3.62, 3.63] described the strain dependence 
of the dynamic properties of various elastomer materials and different com- 
positional blends. Holownia [3.64~J described an experimental method for 
accurately measuring the static Poisson’s ratio of elastomer materials and 
presented test results relating Poisson’s ratio to the quantity of carbon 
black in the elastomer blend. Gupta et al !j3*48j and others [3.49, 3.50, 

3.56, 3.57, 3.65*1 reported the effect of temperature, heat dissipation and 
preload on the dynamic properties of polybutadiene elastomer. They also 
investigated the effect of strain and found this to be a more significant 
and consistent parameter than heat dissipation on the dynamic properties. 
Smalley et al 3.66, 3.67J reported the results of a series of tests using 
the BERM method" for determining the dynamic properties of elastomer 0-rings. 
The effects of temperature, strain, preload (in terms of squeeze and stretch) 
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and elastomer material on the dynamic properties were investigated. Smalley 

3.68 presented a review of results obtained with polybutadiene and other 
elastomer materials as a function of several test parameters and for a 
variety of test hardware configurations. 

The American Society of Testing and Materials has published testing standards 
for elastomers for a large variety of static and general property tests 3.23 . 
These include ASTM-G-395 for compression set, ASTM-D-695 for compressive 
properties and ASTM-C-177 for thermal conductivity. ASTM has also published 
a standard for dynamic property elastomer tests using a Yerzly mechanical 
oscillograph (ASTM-D-945) . This apparatus is a free vibration device which 
is strictly mechanical and operates on the principle of a horizontal rocking 
lever arm. A sample may be tested in either compression or shear. Norrington 

3.69 noted that this apparatus was only suitable for comparative measure- 
ments"' due to the nonlinearity introduced by the continually changing ampli- 
tude of oscillation. In addition, the measurements can generally be made 
over only a very limited frequency range which is determined by the geometry 
of the apparatus and the stiffness of the elastomer sample. Condensations of 
many of the ASTM test standards for elastomers are given in Elastomeric 
Materials 13.70:, which is a general handbook of elastomer material and 
properties. 

3.3 Material Properties 

Numerous references are available which provide collections of general physical 
properties for elastomers. However, very little data is available for dynamic 
properties of elastomers. Most of the elastomer dynamic properties which are 
available from the open literature are found in reports of test programs, such 
as those cited in the previous section. Data regarding the general physical 
properties of elastomers may also be found in published test results. In 
addition, a number of handbooks have collections of these properties in a 
format which is more compatible for design applications. Several handbooks 
are available which are devoted entirely to elastomers. One of the most com- 
prehensive of these is Elastomeric Materials 3.70 . This handbook provides 
a large amount of data for a great variety of "elastomer materials which are 
classified by both generic type and common tradenames. Payne and Scott 3.71 
provided a modest tabulation of elastomer properties in the Appendix to their^ 
book on the application of elastomers. Other handbooks, which are devoted 
entirely to elastomers and contain useful property information, include those 
by McCann 3.72 , Manufacturing Chemists 1 Association 13.73^, Nielsen 3.74j, 
and Simonds" 3 .~7 5~] . 

w j 

General mechanical handbooks may also be found which provide useful elastomer 
property data. Mechanical Behavior of Materials 3.4; provides general physi- 
cal properties including the static modulus of elasticity for a number of 
common polymers. Mechanical Design and Systems Handbook ; ! 3.76 : provides some 
dynamic properties for several elastomer materials. However, this data is 
provided for particular test frequencies, which are generally above one 
kilohertz and, thus, is somewhat limited for use in rotordynamic applications. 
Mark's Handbook for Mechanical Engineers 3.77 • provides a list of environmen- 
tal properties for several types of elastomers^ Additional handbooks which 
provide information on elastomer properties include Vibration and Acoustics 
Measurement Handbook l3.78~^ and Materials and Design Engineering: Materials 

Selector Issue 3.79~\ 
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3.4 Applications 


The static and unidirectional dynamic applications of elastomers are dis- 
cussed in a substantial quantity of literature* However, limited informa- 
tion is available on rotordynamic applications. Applications-oriented 
literature can be considered to include test results and property reference 
literature, such as that cited above. Application of this data, particu- 
larly the dynamic data, may not be straightforward. 

Payne and Scott 3.71 discussed a number of geometric considerations for 
unidirectional loading, such as a compression specimen shape factor. They 
also provided a general design guide for elastomeric mounts. GSbel 3.80 
provides considerable detail useful for the design of elastomer structures 
for unidirectional and torsional loading (e.g., shaft couplings). This 
discussion includes calculation methods, principles of construction and 
numerous examples. McGuire 3.81 and Weller 3.82 described helicopter 
rotor lead-lag dampers with elastomer elements T Nashif and Halvorsen 
f 3 * 83 j discussed the design of layered viscoelastic damping treatments for 
structural damping. Jones and Henderson 3.84J described the design of 
elastomer elements for guide vane support 'in gas turbine engines. Hertz 
[^3.85 discussed the design of 0-ring seals, taking into consideration the 
modulus of elasticity of the elastomer materials. Foster |3.86j discussed 
vibration isolation using elastomers and described the design of an elas- 
tomer support J: or spinning frames in textile mills. The Shock and Vibration 
Handbook 3.87 provides some general design information concerning the use 
of elastomers for vibration isolation. In a series of articles, Nagel 
discussed various aspects of elastomer component design including some 
specific applications 3.88 through 3.97 . Some design information is also 
available from various 'elastomer manufacturers 1 catalogs, such as that by 
the Minor Rubber Company 3.92^. 

The design information available for elastomer dampers for rotating machinery 
is extremely limited. Most of the information that is available is empirical. 
However, analytical correlations are provided in some of the test result 
literature cited above. Finney and Gupta 3.93 applied a finite element 
analysis for elastomer application. However, in general, these analytical 
design tools are expensive to use. 

There are apparently no complete design guides available for elastomer dampers 
for rotating machines. A few references are available which include detailed 
descriptions of actual elastomer dampers, and a number of patents have been 
obtained for elastomer dampers (including References 3.94 through 3.96) in 
1972. In general, the patent disclosures describe the geometric design of the 
damper only and no elastomer specification is given. Potter 3.97 provided 
a general discussion of elastomer damper applications with specific reference 
to lead-lag rotary wing dampers. He stated that these dampers typically have 
a maintenance free service life of 1500 to 2000 flight hours and an operating 
temperature range of from -54° C (-65° F) to 93° C (200°F) . Anwar and 

Kramberger (3.98 discussed the use of elastomer mounts for fluid-film bear- 

ings and described the results of a series of unidirectional load tests. 
Fischer et al [3.99J described the design of a compliant surface spherical 
bearing for a gyroscope which utilized elastomer elements. In a series of 

references, Tecza and others 3.56, 3.100, 3,101 and 3.102 described in 



detail the damper design and results of a series of tests utilizing an 
elastomer damper for a gas turbine rotor simulator. Smalley [3.68} provided 
a review of this effort in relation to its historical background. Zorzi et 
al [3.103] described a similar elastomer damper design used for a super- 
critical power transmission shaft. Details of these and other elastomer 
damper designs for rotating machines are presented in Section 8.0 of this 
handbook. 

3.5 Summary 

A review of the literature concerning elastomer materials and properties 
and the design of elastomer dampers has been presented in this section. 
Particular areas which have been discussed include viscoelastic theory, 
elastomer chemistry, experimental determination of properties, material 
property references, and applications-oriented literature. While the 
applj cations-oriented literature is the most pertinent to the practical 
design of elastomer dampers, it is also the most scarce. A number of the 
topics touched on in this section are covered in more detail in later 
sections of this handbook. 
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4 * 0 DYNAMIC BEHAVIOR OF VISCOELASTIC MATERIALS 


Viscoelastic materials combine the properties of purely elastic materials 
with those of purely viscous materials. For elastic materials, a one-to-one 
relationship exists between stress and strain, and the stress is independent 
of the rate of strain. When an ideal elastic element is subjected to a 
cyclic deformation, the stress and strain remain in phase. For a linear 
elastic material stress is proportional to strain (i.e. Hooke’s Law). An 
ideal viscous material is characterized by a relationship between stress 
and rate of strain and for a sinusoidal cyclic deformation, the stress 
leads the strain by ninety degrees. For a linear viscous element (repre- 
sented by a viscous damper, described in Section 2.0), stress is directly 
proportional to the rate of strain, with the proportionality constant being 
the viscous damping coefficient. 

Viscoelastic materials are characterized by creep and stress relaxation. 

Creep is the process whereby the deformation of a viscoelastic material 
under a sustained load increases with time to an asymptotic limit. Similarly, 

load relaxation is the process whereby the load required to maintain a 
constant strain decreases with time. Stress in a viscoelastic material is 
related to both strain and the rate of strain and, for cyclic deformation, 
is not in-phase with strain. Viscoelastic materials, like elastic materials, 
may be either linear or nonlinear. Warnaka j4.lj* has shown that several 
typical elastomer materials behave as nearly linear viscoelastic materials 
for frequencies in the range of 100 to 1000 Hertz and for dynamic strain 
levels up to 0.1. The elastic and viscous moduli for elastomers tend to be 
functions of frequency, temperature, strain, geometry and material composi- 
tion. The effects of frequency and temperature are inversely related, as 
discussed below. The basic relationships and parametric influences of 
viscoelasticity are discussed in several references including References 
4.2, 4.3 and 4.4. 

4.1 Basic Relationships 

If a load cycle is applied repeatedly to an elastomer element, the process of 
relaxation causes hysteresis - that is, a load deformation curve which depends 
on the sign of the rate of change of load as well as on the load’s magnitude. 
Figure 16 presents a representative elliptical hysteresis loop. Implicit in 
such a load deflection curve or hysteresis loop is a loss of energy per cycle 
equal to the area of the hysteresis loop. Expressed in another way, elements 
with such a load-deflection curve provide energy dissipation or damping (in 
this case,^hysteretic damping, as discussed in Section 2.0). As shown by 
Lazan Jji-.S ,, the material behavior for an elliptical hysteresis loop is 
linear, an3 the hysteretic damping, 1 ^ 2 ^ may be expressed as 

D 

k 2 = — % (4-1) 


* Numbers in brackets indicate References found in Section 4.7 
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Fig. 16 Load-Deflection Curve for an Elastomer Showing Hysteresis as a Result of the Relaxation Process 




where D s is the area of the hysteresis loop and x 1 is the amplitude. 

It is also possible to define stiffness in terms of this loop as the ratio 
between force at maximum deflection and the value x^ of maximum deflection. 

k l = F l^ x l (4-2) 

where F^ and x^ are defined in Figure 16. 


A linear elastomer element undergoing a repeated deflection cycle defined by 


* io)t 
x = x e 


(4-3) 


gives rise to a cyclic force variation described by 

F = F*e 1Ut 

where: to is the frequency (rad/sec) 

x is the instantaneous displacement (m) 


(4-4) 


F is the instantaneous force (N) 
t is time 

it it 

x ,F are complex numbers defining the amplitude and phase 
of the displacement and force variations (and it is 
understood that only the real parts of equations 
(4-3) and (4-4) are used to establish values for x 
and F) . 


The force and displacement may be related by 

F = k*x (4-5) 

A 

where k is a complex number 

k* = k + ik 2 (4-6) 

and is generally a function of frequency, temperature, dynamic strain, and 

geometry, as well as material. The quantities k^ and k^ are generally 

referred to as the stiffness and damping of the elastomer element and are 

likewise functions of the parameters listed above. The damping ability of 

an elastomer element is sometimes expressed in terms of a loss coefficient 

(nondimens ional) which relates damping to stiffness, or the energy dissipated 

per cycle (D ) to the energy stored at maximum displacement (U ) 
s s 



(4-7) 


where k* = k^(l + in) 
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4.2 The Effects of Geometry 


For a general elastomer element, the geometry can normally be characterized 
by a certain area of the element (A) which is stressed and a certain dimen- 
sion of the element, (h) , which is strained. Figure 17 shows two elastomer 
element geometries which have been tested often for dynamic stiffness and 
damping and which illustrate the effects of geometry. For the shear element, 
the stressed area. A, is its width times its length (b*£), and the strained 
dimension is its thickness, h. For the compression element, the stressed 
area is the top area of the cylinder (ir d /4), and its strained dimension 
is the height, h. 

For the shear element, stiffness and damping can be related to effective 
moduli as follows: 

Stiffness: k = ^ (4-8) 

Damping: k 2 = G~ f ^ (4-9) 

where G^^ are shear moduli, commonly referred to as the effective 

storage and loss moduli, respectively, such as 

G* = G f + iG" 


For the compression element, stiffness and damping can be related to the 
effective moduli as follows: 

Stiffness: k_ =* 3 ~ (4-10) 

1 eff h v ' 

Damping: k 2 = 3 G eff £ (4-11) 

where the factor of 3 arises from the classical elasticity relationship 
E = 2(1 + v)G and the fact that elastomer Poisson’s ratios are typically in 
the range of 0.4945 to 0.4999 (Holowina ^4.6 ). Any errors in applying this 
classical elasticity relationship to dynamic behavior of a viscoelastic 
material are implicit in the statement that G" f , G ^ff are effective moduli. 
G eff 5 G'ff are * generally, functions of the mlteriaf, the frequency, the 
amplitude of loading, the temperature, the initial strain and the geometry. 
They are, therefore, component rather than material properties. 


For static loading the effects of bending can be shown [4.2! to cause the 
following relationship between effective and true modulus for the shear 
element of Figure 17. 


(g' ) = g; 

eff o o 


? 2 
1 + h /3JT 


(4-12) 


where Z is the dimension of the sheared area in the direction parallel to the 
loading, h is the thickness, and the subscript, o, implies static loading. 

In Figure 18, the stiffness correction factor (G^ /G') is plotted as a 

function of h/Z for a rectangular shear element. under static compression 
loading, the interaction of axial, radial and circumferential stresses was 
shown by Gent and Lindley [4.7" , Hattori and Takei j*4.8 , and Payne [4.9^ to 
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Fig. 17 Test Element for Shear and Compression Loading 


60 



STIFFNESS CORRECTION FACTOR 



(h /&) — RATIO OF THICKNESS TO LENGTH PARALLEL TO 
DIRECTION OF LOADING 


Fig. 18 Static Stiffness Correction Factor, Rectangular Shear Specimen. 
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give the following relationship between apparent, or effective, static 
modulus and actual, or true, static modulus: 

(G' ) = G' (1 + BS 2 ) (4-13) 

err o o 

where S is a shape factor defined as the ratio of the loaded area to the 
unloaded area (d/4h for a cylinder) and 6 is a coefficient, shown by Gent 
and Lindley jj4.7j to be equal to 2. 

A more complete expression based on the work of Gent and Lindley j4.7j, which 
attempts to account for bulk compressibility effects, is given by: 

G "* 3 q' 

Yr v°V = (1 + es 2 ) -1 + (4-14) 

where B is the bulk modulus f or_ the_^ material. However, this expression has 
been tested by Moghe and Neff [4.10! using classical elasticity solutions and 
shown to be incomplete for large values of S due to an additional dependence 
on Poisson’s ratio. For values of S less than 10, however, neither the second 
term in equation (4-14) nor the inadequacies due to Poisson’s ratio effects 
are significant. 

In Figure 19, the static stiffness correction factor ((Ggff) 0 /G^) is plotted 
against the ratio of diameter to height, as given by equation (4-13) with 
3=2, for a compression loaded cylinder. The existence of such a square law 
dependence on shape factor has been verified experimentally by Payne ! 4.9j» 

He showed that, for natural rubbers, the shape coefficient, 3, was close to 
2 but, for harder materials the coefficient 3 could fall by as much as 50 
percent. Under low-frequency dynamic loading (10 Hz) Payne also demonstrated 
that the effective storage modulus varied with shape factor according to a 
square law. 

The analytical or experimental quantification of the effects of geometry 
(shape) at higher frequencies appears to be limited. Cannon, Nashif, and 
Jones [4.1lj have demonstrated that shape has a pronounced influence on the 
effective storage modulus of cylindrical specimens, in the frequency range 
of 250 to 500 Hz. Smalley and Tessarzik ^4.12^! have evaluated the effect 
of frequency on shape factor for a particular elastomer material and geometry. 

In practice, more complex elastomer geometries are sometimes used such that 
the elastomer elements are not subjected to pure shear or pure compression, 
but rather a combination of the two. A common example of such a practical 
elastomer configuration is a radially loaded ring cartridge. Figure 20 
illustrates the geometry of typical elastomer ring cartridge. The elastomer 
material separates the inner diameter of one cylinder from the outer diameter 
of a second, smaller cylinder. The loading of interest is radial; that is, 
it tends to move one cylinder radially eccentric to the other. 

Three methods of considering geometric effects are reviewed in this section: 

• Beam-Column Method 

9 Method of Gobel i4.13] 
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STIFFNE 



(d/h) — RATIO OF DIAMETER TO HEIGHT 


Fig. 19 Static Stiffness Correction Factor, Cylindrical Compression Element. 
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Fig. 20 Ring Cartridge Element 
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® Predictions Based Upon Plane Stress Analysis 

Each of these is a f static’ analytical method; that is, it seeks a ratio of 
element stiffness to storage modulus (which is a function of the geometry 
only) and applies that ratio in the determination of dynamic stiffness and 
damping. 

4 .2.1 Beam- Column Method 


In this development, the cartridge is assumed to be made up of a series of 
noninteracting, elemental column/beams, whose variation of thickness due to 
radius effects is small. Figure 21 illustrates the cartridge and the Beam- 
Column element, whose mean loaded area is 

da = 7 £d9 (4-15) 


where £ is the length of the cartridge 


r is the mean radius = (r^ + r^)/2. 


Consider the rigid inner member to be displaced radially by a small radial 
distance e, relative to the rigid outer member. The compression deflection 
of the element is (e cos 9) , where 9 is the angle the radial line to the 
element makes with the direction of displacement. The shear deflection of 
the element is (e sin 9). So, considering the element to act independently 
of adjacent elements, the elemental compression and shear forces, df^, df g 


„ r £d0 

df = e cos 9 E r 

c ( r 2 “ r i> 


(4-16) 


df = e sin 9 G r 

S <r 2 - rp 

and the total force exerted at the mean radius is 


(4-17) 


2tt 


F = 


e r £ 
(r 2 ' r l } 


2 2 
E cos 9 d9 + G sin 9 d9 


e r £ it 
< r 2 ' r l } 


E + G (4-18) 


so that, if we assume E = 3G (Poisson’s ratio v = 0.5 implied), then, 
p ( r o r i ^ 

k » ~ = 2 tt "t* V G£ BEAM-COLUMN RESULT (4-19) 

re ( r 2 “ r 1 ) — ■ — 


A slightly different result is obtained if the tapering of beams and columns 
due to radius effects is accounted for. 


4tt G£ 

r £n(r^/r^) 


ACCOUNTING FOR 
RADIUS EFFECTS 


(4-20) 
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But, as Table 3 shows, the reduction in stiffness due to radius effects is 
less than four percent for radius ratios up to a value of two . For most 
practical applications the radius ratio is much less than two* 



The method of Gobel initially follows a very similar development to the 
Beam-Column method accounting for radius effects, with an intermediate result: 


k = - — ~ r v ( G + E ) 
r in (r 2 /r 1 ) 


(4-21) 


However, Gobel then claims that, for -j ~ r- = 1, the effective value of E 

given by r 2 r l 


E ~ 6 • 5 G 
eff 


(4-22) 


and that, for other length to radius ratios, a correction factor, f^, should 
be applied, so that 

7,5 tt£ Gf 


k r to (r 2 /r 1 ) 


(4-23) 


Gobel presents f .. as a graphical function of 1/ (r 2 - r^) . Payne and Scott 
{ 4 . 2 ] in their book ENGINEERING DESIGN WITH RUBBER present a cubic fit to 
Gbbel's function: 

f = 1 + 0.0097 £ 3 /(r 2 - r x ) 3 (4-24) 

which agrees closely with Gobel ? s graphical presentation. The basis for f^, 
or for the statement E e f£ - 6,5 G is unclear from Gobel 1 s work. These are 
not presented as empirical results, since no test data are given for this 
geometry, Gobel T s method implies that E = 6,5 G is an asymptotic condition 
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for very short length cartridges * One set of static test data for the 
cartridge configuration (Adkins and Gent [4.1l]) has been identified, and 
this data is compared with the prediction of the Beam-Column method and the 
method of Gobel in Figure 22 (after Reference 4.12). It shows good agree- 
ment at higher values of 1 / {*2 - r^) between Gobel and test data, but at 
lower Z/ (r 2 - r^) , the Beam-Column method appears to come closer due to a 
much sharper drop in stiffness of the test element. This comparison suggests 
that perhaps E ; ^ = 6.5 G is not a valid asymptotic condition for short 
cartridges. 

4.2,3 Predictions Based Upon Plane Stress Analysis 

The value of using a plane stress elastic analysis has also been investigated 
[4 . 12j . In this case, only a quarter of the cartridge was modeled and 
appropriate symmetry and antisymmetry boundary conditions were imposed. The 
quarter cartridge was embedded in a rectangular finite difference mesh and a 
finite displacement of the inner member relative to the outer member was 
imposed. A computer program was used to solve for the internal displacement 
field within the cartridge. Since it was not possible to employ a Poisson’s 
ratio of 0.5 in this computer program, a value of 0.4983 was used together 
with a shear modulus value of 1.795 x 10^ N/m^. This combination of Poisson’s 
ratio and shear modulus corresponds to a bulk modulus of 5.28 x 10^ N/m^, which 
is within the range of bulk moduli determined by Rightmire jj4. 15] . 

The stiffness was calculated by determining the component of shear and direct 
stress in the direction opposing the displacement, e, and integrating over the 
active cartridge surface to obtain the net force. The stiffness is, then, the 
ratio of net force to displacement. 

Estimating Dynamic Values of Cartridge Stiffness and Damping. Each of these 
analytical methods described for the ring cartridge can be applied to a ring 
cartridge geometry. This has been done for a particular configuration (and 
two lengths) by Smalley and Tessarzik 0.12 j. Their results are presented in 
Table 4 in nondimens ional form; that is, tHe ratio of stiffness to the product 
of shear modulus and length is presented. 


Table 4 


NOND IMEN S IONAL STIFFNESS FOR CARTRIDGES SELECTED 
FOR TESTING (Inner and Outer Diameters: 1.905 and 2.858 cm) 

(after Smalley and Tessarzik Rj-.12~j) 


Method 

k/G£ For Z = 0.476 cm 

k /GZ For Z = 0.954 cm 

Gobel 

58.7 

62.6 

Plane Stress 

35.1 

35.1 

Beam- Column 

31.4 

! 

31.4 
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REDUCED RADIAL STIFFNESS (K/G JL) 




RATIO OF LENGTH TO RADIAL EXTENT 
-&/(r 2 - r-|) 


Fig. 22 Comparison of Prediction Methods for Cylindrical Cartridge 
(After [4.12]) 
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It is seen that only the Gobel Method reveals any dependence of this non- 
dimensional stiffness on length of the cartridge , and that the level of the 
Gobel predictions is substantially higher than that of the plane stress and 
Beam-Column predictions. 


The assumption necessary to apply the Gobel and plane stress predictions in 
the determination of dynamic stiffness and damping is that: 


* 


k 

* 

G l 



(4-25) 


where k = 4* ik 2 is the complex impedance of the cartridge, G =G' + iG^ 

is the complex modulus and £ is the length. That is, the relative distribu- 
tion of stresses is independent of time under sinusoidal loading and displace- 
ment . 


The simple Beam-Column method, without ^accounting for radius effects, may be 
developed directly using either G or G , as it assumes no interaction between 
elements. In other words, no further assumptions are necessary in applying 
it to dynamic stiffness and damping beyond the significant assumptions in its 
development already presented. 

For a series of tests with two ring cartridge specimens, Smalley and Tessarzik 
[^4.12j found that, using the Beam-Column analysis and the method of Gobel, 
stiffness and damping values were predicted which generally bracketed the test 
results. Neither of these methods appeared to provide consistently better 
correlations than the other, 

4,3 The Effects of Frequency 

As might be expected from a material whose behavior is governed by a time- 
dependent relaxation process, the slope and area of the hysteresis loop, 
together with the stiffness, damping, and loss coefficient, are functions of 
the frequency with which the load cycle is executed. It is generally observed 
that the dynamic stiffness is higher than the static value and that stiffness 
tends to increase with increasing frequency. Shown on a scale which may cover 
many decades of frequency, the general form of elastomer behavior variation 
has three distinct regions [4.16J which are illustrated in Figure 23: 

• Rubbery 

• Transition (leathery) 
m Glassy 

In the rubbery region, the behavior is, as implied, "rubbery". Elastomer 
elements are resilient and the damping is usually of the order of 10 to 30 
percent of the stiffness. The hysteresis loop is a rather narrow ellipse. 

As the frequency increases, the transition region is encountered. Both 
stiffness and damping increase more rapidly in this region but the damping 
moves ahead of the stiffness and the hysteresis ellipse becomes broader. 

For a limited frequency range, damping (h^) becomes higher than stiffness 
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Fig. 23 Regimes of Elastomer Dynamic Behavior 
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(k-^) and the loss coefficient may become greater than 1,0, For further 
increases in frequency, the stiffness increases more than the damping, and 
in the glassy region both quantities level off. The rubbery and transition 
regions are those normally encountered in elastomers and offer the best 
opportunity for use in vibration control. 

4.3.1 Representation of Frequency Dependence 

Observed variation of stiffness and damping with frequency gives rise to 
discrete data points. For mathematical purposes it is convenient to seek, 
for an element, a continuous relationship between force and displacement 
(or stress and strain) which closely matches observed behavior. There are 
several possible approaches to this problem, with different levels of sophis- 
tication and convenience. 


Generalized Viscoelastic Models . One approach in representing the frequency 
variation of stiffness and damping of an elastomer element is to replace the 
elastomer with a mechanical analog, consisting of a series of springs and 
viscous dashpots. The simplest models are the Maxwell and Voigt elements 
(Figure 24). Using these basic elements as building blocks, more complex 
systems of springs and dashpots may be developed, such as those shown in 
Figure 25. 

The application of a triparameter model, shown in Figure 26, to elastomer 
specimen test data was demonstrated by Chiang, Tessarzik and Badgley J4.ll]. 
For a selected general viscoelastic model, the methodology to determine the 
element values to best fit a particular variation of stiffness and damping 
with frequency was developed and demonstrated by Gupta, Tessarzik and 
Cziglenyi !4.18j. 

Generalized Force - Displacement Relationships . An alternative approach for 
representing the variation with frequency is to assume a general linear 
differential relationship between force and displacement: 


a F 
o 


+ a n 


dF 

dt 


+ a n 


d 2 F 


dt 


2 + a 3 


d 3 F 

dt 3 


, dx 
c x + c- -r - 
o 1 dt 


UA , 

TT + C 


A 

2 dt 2 


. (4-26) 


and for sinusoidal loading and displacement: 

F = F*e it0t 

* io)t 
x = x e 


So that equation (4-26) becomes: 


(4-27) 

(4-28) 


2 3 2 

(a + iwa. — to a n — ico a 0 + . ..)F* = (c + icuc- - w c 0 - ...)x* (4-29) 

o 1 2 3 o 1 Z 


and 


F* 

x* 


2 . 3 

c + iajc- - a) c 0 - ia) c n + . . , 
o 1 2 3 


a + ioua- 
o 1 


2 

a) a 0 


(4-30) 


ia) a^ + ... 
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A close connection exists between this generalized force displacement 
relationship and the generalized viscoelastic models since, for a given 
order of expansion on left- and right-hand sides of equation (4-26), there 
is a corresponding combination of springs and dashpots (4.16j which will 
give the same force-displacement relationship. 


For example, the differential equation for the triparameter model (Figure 26) 
is 

F + — — = k y + — ( k =k)^ 

k T1 dt T2 y k T1 U T2 TV dt (4-31) 

so that the nonzero coefficients of equation (4-26) for the triparameter 
model are 


a - 1 
o 

a l = b TF k Tl 
C o = k T2 

C 1 = b Tl (k T2 + k Tl )/k Tl 


(4-32) 


Direct Fit of Mathematical Expression . A direct approach to this problem is 
to impose upon the data some correlation function and to seek those coeffi- 
cients which give the best match between the function and the data. Polynomial, 
trigonometrical, or power-law relationships are some of the options available 
for this approach. It was found by Gupta et al [4.10 that the use of a 
power-law relationship was quite adequate to represent the behavior of elements 
of polybutadiene in the frequency range 100 to 1000 Hz. For engineering pur- 
poses, a known relationship between stiffness, damping, and the parameters 
discussed in this chapter for elastomer elements of different geometries is 
all that is required to use these elements in design for control of machinery 
vibrations. In essence, this approach may be considered as the use of a 
generalized Voigt element in which the storage and dissipative elements are 
a continuous function of frequency, temperature and strain. 

Other methods which deal with the frequency dependence of viscoelastic materials 
have been developed. The reader is referred to recent publications such as 
! 4.281 which deal with the constitutive behavior of viscoelastic media by 
application of generalized derivatives. 

4.4 Influence of Temperature 

In general, elastomer behavior is strongly influenced by temperature. Thus 
complete characterization of elastomer elements must provide their dynamic 
behavior over the frequency and temperature ranges of interest. In consi- 
dering the effects of temperature, the following observations are often made: 

® Temperature tends to be most influential where frequency 
is also strongly influential. 
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® Increasing temperature tends to produce the same results 
as decreasing frequency and vice versa. 

Starting from these observations, a method of reducing the dynamic behavior 
of elastomers (and other polymers) covering a range of frequencies and 
temperatures to a single curve relating appropriate reduced properties was 
developed by Ferry, Fitzgerald, Grandine and Williams [4.19^. The method, 
sometimes referred to as the WLF law, is derived from a paper by Williams, 
Landell and Ferry [4.20 in which the mechanics of its implementation are 
described and demonstrated. The mathematical statement of this method is 
as follows: 


(G r } (o,T " (G r ) a Q w,T 
0 c 


(4-33) 


Log in a T = 


-8.86 (T - T ) 
c 


=10 T 101.6 + (T - T ) 

c 


T p 

(G ) T - (G) a 7n~ 
r to, I to , 0 Ip 

c 


where T^ is a characteristic temperature given by 

T = T + (50 + 5)°K 
c t — 


(4-34) 


(4-35) 


(4-36) 


and 


T is the transition temperature between the rubbery 
and glassy states of the material. 


G is the observed modulus at temperature 0 and 
frequency 

G is the reduced modulus, 
r 

is the coefficient by which test frequency is 
multiplied to obtain the corresponding reduced 
frequency at temperature 0^ . 

p,p are the densities of the elastomer material at the 
c 

specified temperature and the characteristic 
temperature, respectively. 


By these relationships, knowledge of a modulus as a function of frequency at 
several temperatures may be interpreted as knowledge of reduced modulus over 
a wide frequency range at the characteristic temperature. This knowledge may 
be further interpreted to give the variation of modulus over a wide frequency 
range at any temperature between T^ and the maximum temperature tested. 

The method has been demonstrated to be applicable to material properties of 
storage and loss moduli or to element properties of stiffness and damping. 

The additional information needed to implement and use the methods beyond 
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the elastomer dynamic data is the value of T c . The transition temperature, 

T t , may be measured independently or T c may be inferred from a sufficiently 
large volume of dynamic data. With an assumption of constant density (p - p ) , 
equation (4-35) is often simplified to give 


(G ) 


r'to.T 


^w,T T 


(4-37) 


In cases where high filler loadings are present 4 . 21) , equation (4-34) may 
not be an adequate representation of the relationship between temperature 
effect and frequency effect. In this case, it may be necessary to use the 
more general representation 


C_ (T - T ) 

T _ I c 

Og 10 a T C 0 + (T - T ) 
2 c 


(4-38) 


where and C 2 (and if necessary T c ) are evaluated for best fit from test 
data by an iterative numerical procedure. Alternatively, T c may be specified 
by equation (4-36) or chosen as some convenient temperature. The determina- 
tion and use of Log^a T is discussed in more detail in the chapter on elastomer 
properties. 


4.5 The Effect of Static and Dynamic Strain 


Large values of static strain (preload) can have a significant effect on the 
dynamic properties of elastomer elements. However, for simple geometries 
(uniform strain and compression elements) the effect of modest preload levels 
(five percent or less) is generally not large. On the other hand, the effect 
of dynamic strain is generally quite pronounced f4.12, 4.22, 4.23, 4.24 and 
4.25!. This effect has been shown to be mainly dependent on two factors: 
chemical and thermal 4.24, 4.25"'. 

The chemical factor involves the way the bonding matrix is disrupted when 
strained. Two mechanisms are involved in this process. One, often referred 
to generally as strain softening, is substantially reversible in short times 
and involves the physical configurations of the particles in the moving matrix. 
Strain softening occurs in any elastomer sample but is most pronounced in 
materials with high carbon black loading. The other mechanism, referred to 
as structure, is substantially irreversible over periods of an hour or less 
and involves the weak chemical and physical interactions between the particles 
at rest. The division between these two mechanisms is not well defined and 
their effects have not been generally quantified. 


The thermal factor, involving self-heating of the elastomer element, is most 
significant at high strains and relatively unimportant at low strains. This 
self-heating, due to thermal dissipation, has been analyzed and quantified 
[4.12 but has been shown not to account entirely for the effect of dynamic 
strain (which is reasonable, since strain softening was not considered) 4.26 . 
The effect of self-heating can be demonstrated with reasonable simplicity by 
an analysis of the thermal dissipation in an elastomer element under uniaxial 
stress [4.12 j. 

A model intended to tie the effects of high dissipation to the change in 



properties has been developed >4.12j. This model simply needs to have 
available values of and G" as functions of frequency and temperature- 
The hypothesis of this model is that large amplitude vibrations, when 
imposed upon an elastomer element, give rise to dissipation of energy 
throughout the element. As a result, there is a tendency for the elastomer 
to heat up. The temperature distributions in the elastomer are governed by 
the equations of thermal conduction and appropriate boundary conditions. 

As a result of the distribution of temperature in the elastomer and the 
dependence of elastomer moduli on temperature, there is a distribution of 
material properties through the elastomer. Load equilibrium within the 
elastomer, accounting for this distribution of properties, governs the 
forces generated by the elastomer and its apparent stiffness and damping* 

The implementation of such a model in two or three dimensions could be very 
complicated. However, to typify the results given by such a model and to 
identify the major items of information needed while minimizing the com- 
plexity, a one-dimensional development has been undertaken 4.12 . An 
important assumption is that the time scale associated with the vibrational 
phenomena (less than 0.01 second) is short compared to the timg scale assoc- 
iated with changes in the thermal condition of the material (h p C p /k, which 
is of the order of 40 seconds and higher) . 

4.5.1 One-Dimensional Model for Large Dissipation Effects 4.12 

This model is directed at either a shear or a compression element in which a 
uniaxial state of stress prevails. The temperature distribution is one-dimen- 
sional and the flow of heat is normal to the stressed area. 

As a result of imposed conditions, an elastomer element in the form of a 
right circular cylinder is undergoing time varying changes in temperature, 
internal stress and displacement. The cylinder is as shown in Figure 27. 

It is assumed that stresses are uniaxial. 

At any point 

a = a(Z,t)e 1Ut (4-39) 

U z = U(Z,t)e ia)t (4-40) 

T z = T (Z,t) (4-41) 

where a is the normal axial stress 
z 

is the axial displacement 
T^ is temperature; 

and it is understood that only the real parts of equations (4-39) and (4-40) 
apply. a and U are complex numbers containing both amplitude and phase infor- 
mation. It is assumed that the amplitude functions a,U are undergoing much 
slower changes with time than a z and U z . In Reference 4.12 , the governing 
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Fig. 27 Compression Loaded Cylinder 



equations for this thermo-viscoelastic model for compression (Figure 27) 
and shear loading (Figure 28) are derived to give: 

Energy equation : 
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Equilibrium equation : 
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where y - 1 for a shear element, and for a compression element, y = 3. 

The additional equations to be satisfied are: 

Boundary conditions : 

T = T a at Z = 0,h (4-44) 


U = 0 at Z = 0 


(4-45) 


For imposed displacement at top of element : 

V = x* at Z = h (4-46) 

For base excitation, equation of motion : 

G* Air - Mg) 2 U = m Z x* (4-47) 

dZ O 

where ^ 

x is the amplitude of displacement of top relative to base 
* 

x q is base excitation amplitude (absolute) 

M is the supported mass 
a) is the base excitation frequency 
A is the stressed area 
T q is the temperature at base and top 
cf is defined above. 

4.5.2 Solution of the Equations 

As a result of the dependence of G* on temperature, the system of equations is 
nonlinear. While steady-state results are of most interest, solution in the 
time domain has been selected as more convenient to implement than an itera- 
tive solution. The solution method is a modification of the second-order 
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Crank-Nicholson Method [4.25], in which nonlinear terms are evaluated by 
extrapolation forward one-half time-step. 


The energy equation is written to apply at the time t + At/2. 
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(4-48) 


The left-hand side has a linear operator and consists of known quantities at 
time, t, and unknown quantities at time, t + At. The right-hand side is 
unknown since it is a function of temperature at time t + At/2, but can be 
approximated by the linear extrapolation from known quantities: 


W ■ (3/2)R t - <1/2) Vit 


(4-49) 


The equation for U is solved at time t + At, once and the associated 

values of G* are known. 

The time integration solution is continued until results for successive time 
steps indicate that thermal equilibrium has been reached. The equilibrium 
distribution of temperature, material properties, and deflection will deter- 
mine the steady-state behavior which will be observed once the initial thermal 
transients have died out. 


An interesting feature of thermal transient behavior occurs when the base 
excitation boundary conditions apply. The frequency of excitation and sup- 
ported mass can be such that the element causes subresonant vibrations when 
the element is cold but vibrates above resonance when the element gets hot 
due to self-heating. In this case, there is a wide swing in phase angle 
during the thermal transient and an inability to hold subresonant conditions 
for levels of dissipation higher than a certain limit. Observed stiffness 
and damping of the element at any time during the time integration, including 
equilibrium, are calculated as: 
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It is important to note that a number of material properties are needed to 
implement this analysis: in particular, density, specific heat and conduc- 

tivity. In addition, if this analysis is to be applied with significance 


81 



to compression elements in which the stiffness correction factor (for 
geometry) is of order 1.5 and above, some account must be taken of the stiff- 
ness correction factor. An assumption may be made that the stiffness correc- 
tion factors follow the functional dependence on frequency and geometry 
determined empirically for cylindrical elements and that the increase in 
effective modulus should take place on a local basis, and be reflected in 
the equilibrium equation: 
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(4-52) 


In any case, to be truly representative, the results of this analysis must 
be modified to include the effects of strain softening and structure. 


4 . 6 Summary 

Elastomer materials are generally represented as linear viscoelastic materials 
with storage and loss moduli which are dependent upon frequency, temperature, 
strain, geometry and material composition. The dynamic properties of visco- 
elastic elements are characterized as complex stiffness, where the real part 
is the typical stiffness and the imaginary part is the hysteretic damping 
coefficient. The loss coefficient is defined as the ratio of the imaginary 
part of the stiffness to the real part. Similarly, the basic dynamic property 
of viscoelastic materials is the complex shear modulus. The real and imaginary 
parts of the shear modulus are the storage and loss moduli, respectively. 

Since elastomer materials are essentially incompressible, the Poisson 1 s ratio 
is very nearly equal to 0.5, and the complex compressive modulus is very nearly 
equal to three times the complex shear modulus. The effects of geometry, 
frequency, temperature and strain on the dynamic properties of elastomers are 
very significant. 

The effects of geometry have been considered for three basic loading configu- 
rations: pure shear, pure compression, and combined shear and compression in 

a ring cartridge. The effect of bending on a shear specimen with a length- to- 
thickness ratio of at least four has been shown to be negligible. The effect 
of end conditions on a compression specimen has been shown to be significant 
for diameter- to- thickness ratios greater than 0.5. Since it is generally not 
advisable to use compression specimens with smaller ratio’s than this which 
may have the tendency towards column buckling, the effect of end conditions 
becomes important. Several analytical methods have been used to evaluate the 
effect of combined shear and bending in a ring cartridge. Some of these 
methods have been shown to give results which differ by as much as a factor 
of two. Experimental correlations have tended to fall in the range between 
the analytical prediction extremes, with none of the analytical methods 
giving consistently superior correlations. 

Over many decades of frequency, elastomer behavior is characterized by three 
regions: rubbery, transition and glassy. Various modelling approaches have 

been used to represent the frequency dependence of elastomer dynamic proper- 
ties including generalized viscoelastic models, generalized force-displacement 
relationships, and direct fit of mathematical expressions (particularly using 
power- law relationships) . 

The effect of temperature on the dynamic properties of elastomers has been 
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found to be opposite that of frequency. That is, increasing temperature 
tends to produce the same results as decreasing frequency and vice versa. 

A method of reduced variables, sometimes referred to as the WLF law, has 
been devised to combine the effects of temperature and frequency using a 
single reduced variable. 

The effect of static strain (preload) on the dynamic properties of elastomers 
has been shown to be relatively insignificant for moderate values of preload. 
The effect of dynamic strain, which is extremely significant, can be attri- 
buted to two factors: a chemical factor (strain softening and structure) 

which is most significant at low to moderate strain values, and a thermal 
factor (self-heating) which is most significant at high strain values. While 
strain softening and structure are difficult to quantify, an analysis of 
self-heating is straightforward to a simple geometry and a uniaxial load. 
Using such an analysis, self-heating, as expected, has been shown to account 
for some, but not all, of the effects of dynamic strain. 
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5 . 0 MEASUREMENT OF DYNAMIC PROPERTIES OF ELASTOMERS 


Various test methods have been used for determining the dynamic properties of 
elastomers. Experimental determination of these properties is necessary to 
correlate with, and supplement, the results of analytical procedures such as 
those discussed in Section 4,0. While providing important insight into the 
behavior of elastomers, analytical studies generally have been limited in 
their characterization of the effects of many significant parameters. Con- 
sequently, the emphasis has been on empirical determination of elastomer 
dynamic properties; this trend will probably continue into the near future. 

In general, empirical relationships for elastomer dynamic properties involve 
the use of a Voight model (Section 4.0) where the damping is assumed to be 
hysteretic (Section 2.0) and is represented as the imaginary part of the 
complex stiffness, k*. Since it is not possible to measure k* directly, it 
is necessary to measure certain physical quantities (such as force and dis- 
placement) and compute k*. If general material properties (rather than test 
element properties) are desired, storage and loss moduli may be calculated 
from k* by considering the appropriate geometric relationships, as discussed 
in the previous section. The values of the storage and loss moduli for a 
given test element will be functions of frequency, temperature and strain. 
Empirical relations may be formulated to represent the effects of these 
parameters. 

When designing an elastomer damper, it is desirable to use existing dynamic 
property data for elastomers (such as that provided in Section 6.0). However, 
this is not always possible. For example, if the design requires an elastomer 
material (or operating or environmental condition) for which the appropriate 
data is not available, it may be necessary to generate this data through a 
test program. This may also be desirable where special elastomer configura- 
tions are to be used or when critical applications require extremely accurate 
data and the elimination of variations due to minor compositional changes, 
which can influence the dynamic properties Q> . 1J ^ . It has been shown that 
translational load test results are adequate for predictions of properties 
for rotating load dampers [5.2~j. However, consideration of the effects of 
both compressive and shear deformation for appropriate geometries, such as 
radial buttons (Sections 6,0 and 8.0) is required. 

Dynamic test systems may be considered to consist of two basic subsystems s 
1) the mechanical hardware, with associated control electronics, that 
exercises the test specimen, and 2) the sensors and measurement electronics 
(generally referred to as instrumentation) which acquire the basic vibration 
data from which specimen properties are determined. Advances in electronics 
and, in particular, the improved capability of measuring phase angles between 
two signals have influenced a gradual shift towards simpler mechanical systems 
These systems have fewer components of uncertain dynamic characteristics with 
more sophisticated electronics. Recently, the instrumentation subsystem has 
been expanded to include minicomputers for; test control, data acquisition 
and calculation of dynamic properties. Additional discussions on test 
instrumentation may be found in [5*3, 5.4 and 5.5J. 


t Numbers in brackets indicate References found in Section 5.7. 
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In this chapter, cursory descriptions are presented of a number of commonly 
used test methods. These test methods may be separated into three categories: 
forced resonant tests, free resonant tests and forced nonresonant tests. 

Details of those and other test methods may be found in References 5.6 
through 5.15. 

Quick estimates of dynamic stiffness of elastomer specimens can be made from 
static stiffness values. Static stiffness can either be measured directly 
from static load tests or Shore ® (durometer) hardness, which can be_measured 
very easily using a durometer tester. The relationship between Shore ® A 
hardness and static shear modulus is illustrated in Figure 29 5.16 . Dynamic 

stiffness may then be estimated using frequency multipliers which are generally 
in the range of 1.5 to 5 for low frequencies (larger for higher frequencies) 
and low dynamic strain (smaller for high strain) . These multipliers may be 
deduced from static and dynamic properties of generic materials, although 
they are generally functions of filler loading, as well as of elastomer 
material. For example, a frequency multiplier of about four has been used 
for a particular flourocarbon elastomer mix. However, for the engineer who 
requires elastomers for their damping characteristics, a static measurement 
offers no information on energy dissipation. 

5.1 Forced Vibration Resonant Test Methods 


These test methods employ resonant amplification of test specimen deformation. 
Most of these test machines are essentially single-degree-of-f reedom vibratory 
systems, where the elastomer specimen provides most of the flexibility and 
damping and the system inertia may be fixed or variable. For fixed-inertia 
machines, the resonant frequency can be changed only by changing the elastomer 
specimen material or geometry. Most of the early test machines fall into this 
category due to the lack of availability of precision instrumentation (par- 
ticularly for measuring phase angles) . Most of the forced resonant methods 
permit data to be taken only at resonance. The elastomer specimen stiffness 
and damping are calculated from the resonant frequency and bandwidth and the 
characteristics of the test machine, with no need for phase angle measurement. 
With advances in instrumentation, newer methods evolved, such as the Base 
Excitation Resonant Mass (BERM) method. The BERM method uses accurate phase 
angle measurement to broaden the frequency range at which data can be taken 
and reliable data may be acquired at off-resonance conditions. 

Two forced resonant methods which use fixed inertias are the Tangorra Hystere- 
simeter [*5.17 and the vibrating reed 15.18, 5.19!. The Tangorra Hysteresimeter 
is essentially a dynamic hardness tester (Figure ~^0) , It consists of a vertical 
rod made of nonmagnetic material and fitted with a Shore A hardness indentor 
at the lower end. Attached to the rod are two permanent magnets, each of which 
is surrounded by a coil. One magnet/coil combination is used to excite an 
oscillatory vertical motion in the rod, and the other is used to measure the 
motion. The frequency is varied around the system resonance, and the stiff- 
ness and damping of the specimen are calculated at the resonant frequency. This 
instrument has the advantages of low cost, fast operation, and insensitivity 
to sample size or configuration; however, it is restricted to low strains and 
a single frequency for a particular specimen. 




Registered Trademark, Shore Instrument and Manufacturing Co., Inc. 
Jamaica, New York. 
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Fig. 29 Dependence of Static Shear Modulus 
on Hardness (After Reference [5.16]) 
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A vibrating reed* shown schematically in Figure 31, is generally used for 
plastics or other "stiff" materials . However, by bonding an elastomer 
specimen to a metal strip or by adhering a small piece of metal at the 
bottom of an elastomer strip to respond to the exciting transformer, the 
vibrating reed may be used for elastomers. The storage and loss elastic 
moduli are calculated from the resonance frequency location, bandwidth, 
specimen dimensions and related machine constants. The vibrating reed has 
the advantages of being easily enclosed in an environmental chamber and 
suitable for high damping materials. However, it is limited to low strains 
and a single frequency for a particular specimen. Considerable operator 
skill is required for acquisition of accurate data. 

Several forced resonant test methods have been used which provide for variable 
system inertia so that more than one frequency can be considered for a single 

test specimen. A machine described by Moyal and Fletcher Q5.2(Tj and shown 

in Figure 32 is a typical example of a forced resonance machine employing 
electrical vibrators for excitation. As seen in Figure 32, loudspeaker 
magnets (C) are used in conjunction with alternating current coils (B) on a 
steel bar (A) supported by springs (D) to produce a cyclic force to vibrate 
the elastomer sample. The frequency of the alternating current, and thus of 

the vibration, is varied to find the mechanical resonance of the system. The 

stiffness and damping of the test specimen are calculated from the amplitude 
and frequency of resonance and the driving current and mass of the system. 
Alternatively, the frequency can be left constant and the resonance found by 
varying the suspended mass (M) . 

Forced resonance test machines which take advantage of the properties of 
electro- and magneto-striction in metals for producing a cyclic deformation 
force, have been developed. These properties characterize the change in 
length of metals which occurs when it is subjected to an electric or magnetic 
field. In essence, resonance frequency and bandwidth of a metal rod under- 
going electro- or magneto-striction are modified slightly by pressing an 
elastomer test piece against the end of the rod. The storage and loss moduli 
of the elastomer test piece are calculated from the changes in resonant fre- 
quency and bandwidth. There are restrictions on test piece size and fre- 
quency is not continuously variable^ An example of a magneto-strictive 
method was presented by Nolle jj> . 21 examples of electro-strictive methods 
were presented by Dietz et al > 5,22 and Rhinehart [5.23,. 

In the middle 1950s, engineers at Chrysler Corporation developed the resonant^ 
beam tester for measurement of dynamic properties of elastomeric springs ;5.24j. 
A simplified line drawing of the resonant beam tester is shown in Figure 33. 

The concept of the resonant beam tester took advantage of two facts: 1) through 

mechanical amplification by the beam (lever) ratio, the force required to 
exercise the specimen could be reduced and the velocity sensor sensitivity 
increased; and 2) by operating the system at resonance (with the test specimen 
acting as the system suspension and the beam acting as the mass) , the magni- 
tude of the excitation force could be greatly reduced. Best of all, the 
analytical relationships at resonance between test parameters and test speci- 
men stiffness and damping were reduced to very_ simple forms which excluded 
dif f icult-to-measure phase relationships (j>.14 f . The three measured quantities 
were forced and displacement amplitudes and angular frequency. However, 
experience proved that repeatable results were difficult to obtain in spite 
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Fig. 31 Vibrating Reed (After Reference [5.8]) 
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Fig. 32 Forced-Vibration Resonance Machine for Dynamic Tests (After Reference [5.20]) 
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Fig. 33 Resonant Beam Testing Machine (After Reference [5.9]) 
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of considerable refinement to the mechanical system. Inclusion of test 
specimen support stiffness and beam pivot stiffness and damping was shown to 
be necessary to eliminate significant errors (Figure 33) [3.14*1. Correct 

determination of test specimen stiffness and damping with all system con- 
stants included requires an iterative solution, manageable only by a digital 
computer. Furthermore, operator error has shown to be significant in the use 
of an oscilloscope to determine the resonance condition [3 . 11~| . 

Perhaps one of the earliest dynamic test machines , predating the resonant beam 
tester, was a test rig that incorporated many desirable features. This rig 
was, in many ways, the forerunner of today’s advanced test machines. The 
schematic in Figure 34 shows a supporting frame that holds the test specimen 
in an in-line position with the axis of a shake table. The test specimen 
was loaded by a sizeable mass, and the stiffness and damping could be readily 
determined at the system resonance frequency. 

By introducing the concept of variable mass so that the test frequency range 
may be covered by a number of resonant frequency points by replacing the 
velocity transducer by accelerometers mounted directly on either side of the 
elastomer, and by providing a displacement sensor for amplitude measurement 
across the elastomer, one can envision the evolution of advanced testing 
facilities , 

There were at least two reported applications of the base excitation resonant 
techniques to elastomer material testing in the period between 1958 and 1968 
[5.25, 5.26^). In both cases, electromagnetic shakers were used to excite an 
elastomer-supported mass system. Measured input-output accelerations were 
used for the calculation of material properties. In the earlier version, the 
mass was not varied in the vibrating system and a fairly complex, custom-made 
analog instrument was used for data reduction. In the second application, the 
mass was variable and test data was acquired only for the system resonance 
condition, thus making phase angle measurements unnecessary and simplifying 
the data reduction process. An interesting variant of the base excitation 
method was reported in Reference 5.27; a right circular cylinder was excited 
in torsion, and input-output measurements were made on both sides of the sys- 
tem resonance frequency. 

Progress in the development of accurate digital phase meters and in the avail- 
ability of phase-matched dual-channel tracking filters prepared the way for 
the development of the Base Excitation Resonant Mass (BERM) test method j3.28j. 
A schematic of the test rig, Figure 35, shows a unidirectional vibrating mass- 
spring system driven by an electromagnetic shaker. For either large test 
specimen preloads or very large masses (for low frequency testing) , two 
radial guide bearings and two air cylinders are provided. These elements 
have small but finite stiffnesses and parasitic damping which must be indepen- 
dently measured and accounted for in the data reduction process. The full 
dynamic model of the test rig and the elastomer material test samples is 
shown in Figure 36. A significant feature of this test method is that data 
can be taken throughout a broad range of frequencies about the system resonance 
rather than only at resonant frequency. A detailed discussion of this test 
method is presented in a later subsection. 
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Fig. 34 Bootstrap Dynamic Rate Tester (After Reference [5.24]) 
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Fig. 35 Schematic of Base Excitation, Resonant Mass 
Test Rig (After Reference [5.28]) 
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Fig. 36 Dynamic Model of Base Excitation, Resonant Mass 
Test Rig (After Reference [5.28]) 
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5 .2 Free Vibration Resonant Test Methods 


These methods employ a rigid mass supported by a flexible elastomer specimen. 
The elastomer specimen is initially deformed and then released so that the 
mass is permitted to oscillate. These transient oscillations decay with 
time due to the damping in the elastomer. The elastomer properties are deter- 
mined from the form of the damped oscillation. The frequency of the oscilla- 
tion is used to calculate the elastomer stiffness, and the logarithmic decre- 
ment (Section 2.0) of the decaying waveform is used. to calculate the damping. 
In general, these test methods have the advantages of simplicity, ease of 
application, and suitability for use in environmental chambers. The princi- 
pal disadvantage is that, since the amplitude of oscillation (dynamic strain) 
is constantly changing, nonlinear elastomer behavior may not make it possible 
to account for effects of dynamic strain. These methods are also generally 
restricted to relatively low frequencies (less than 100 Hz) . 

The most commonly used test machine of this type is the Yerzley Oscillograph, 
shown in Figure 37, This device has been standardized by the American 
Society for Testing and Materials (D945-72) 0.29^. The Yerzley Oscillograph 

is strictly mechanical and operates on the principle of a horizontal rocking 
level arm. An elastomer specimen can be tested in either compression or 
shear. A pen attached to the end of the lever arm traces a damped sinusoidal 
curve on a rotating drum chart. This instrument can also be used for measur- 
ing the static stiffness of elastomer specimens. 

A device even simpler than the Yerzley Oscillograph is the torsion pendulum 
jj5.18j, which is illustrated schematically in Figure 38. It consists of a 
suspended elastomer specimen which is clamped at the bottom to a moment of 
inertia bar. The pendulum is given an initial angular displacement and is 
released. A damped sinusoid is produced by measuring torque, linear dis- 
placement, angular displacement, or acceleration, and the dynamic properties 
of the elastomer are calculated from the natural frequency and logarithmic 
decrement. A more refined version of the torsion pendulum has been developed 
by the Research Association of British Rubber Manufacturers (RABRM) p>,15j. 

The operation of this device is essentially as described above, but it has 
been improved to provide increased horizontal stability and easily modified 
inertia, 

5.3 Nonresonant Forced Vibration Test Methods 


These methods involve the application of a sinusoidally varying force to an 
elastomer specimen under nonresonant conditions. In general, both stress 
and strain (force and deflection) are measured and used to calculate the 
dynamic properties of the elastomer specimen. These machines are generally 
restricted to low frequencies (less than 100 Hz). In general, within the 
limitations of these machines, frequency and strain are continuously variable. 

With the earliest machines, the general procedure was to record the stress- 
strain hysteresis loop and calculate the elastomer dynamic properties from 
the dimensions of the ellipse (Section 4.0). These early machines include 
those described by Painter j3 . 3(51 (Figure 39), Roelig Jj>.3l] (Figure 40), 
Lazan Q>.32fJ (Figure 41), and Fletcher and Gent 0.33J (Figure 42). Somewhat 
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Fig. 37 Yerzley Oscillograph (After Reference [5.8]) 


99 



100 







Fig. 39 Forced-Vibration Nonresonant Test Machine 
(After Reference [5.30]) 
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Fig. 40 Roeling Dynamic Test Machine, Forced-Vibration 
Nonresonant Type (After Reference [5.31]) 
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Fig. 42 Forced-Vibration Nonresonant Test Machine (After Reference [5.33]) 




more recently, a test machine was developed by Payne 5.34— (Figure 43) 
which has seen extensive use and is often referred to as the RABRM Sinusoi- 
dal-Strain Machine. This machine utilizes a rather elegant device for 
measuring the phase relationship between stress and strain. From the 
amplitudes of stress and strain (force and deflection) and their phase 
relationship, the dynamic properties of the elastomer specimen can be calcu- 
lated directly. A rather specialized machine, referred to as a dynamic 
viscoelastometer , was designed by Takayanagi Tj.35 ; (Figure 44) and is used 
for testing rubber-band-shaped elastomer specimens in tension. 


Several test machines which have been developed within the past 15 years are 
grouped under the general classification of transmitted force machines (or 
methods). The transmitted force measurement method, as it is commonly used, 
has two main features: 1) the force and displacement measurements are made 

directly across the test specimen, thus eliminating all fixture compliances 
(which were significant for the resonant beam test rig described previously) 
and, 2) an electronic analog circuit that permits the calculation of speci- 
men stiffness and damping at frequencies other than the system resonance. 

(It appears that adoption of the first feature alone to the resonant beam 
test rig would have been a significant improvement.) The development of the 
necessary electronics for the transmitted force measurement method was a 
lengthy and difficult process f 5 * 13, 5.4 i. This was primarily due to the 
need to handle very precisely and very small phase angles between measured 
quantities. This approach, once chosen, also affects the choice of instru- 
mentation and sensors. For example, the use of accelerometers as sensors 
had to be rejected due to both the inherently low phase shift between force 
and deflection (for most elastomer applications) and the difficulty of 
matching two accelerometers and their associated signal conditioning equip- 
ment for such small angle measurements j5.il • The consequences of this 
problem is demonstrated by the following: 


"Consider an elastomer having a stiffness of 1.75 x 10 N/m 
and a damping coefficient of 1.75 x 10^ N-sec/m, that is 
excited at 5 Hz. If the phase error in the measurement is 
0.1°, the computed value for the damping coefficient has an 
error of approximately 5%. Considering that it is not unusual 
to have a phase measurement error of at least 0.5°, this can 
be a serious problem." 15.36 | 


For reduction and improved accuracy in the electronic computation process, a 
minicomputer has been successfully employed with the transmitted force measure- 
ment method (5.5^). 

A typical example of a transmitted force machine is the closed-loop electro- 
hydraulic test machine introduced in 1966 by MTS System Corporation. This 
machine was a logical pairing of the transmitted force measurement method 
with a mechanical test rig that used a hydraulic actuator and servovalve. 

This rig replaced the small electromagnetic shaker and the mechanical multi- 
plier of the resonant beam test rig. A system schematic is shown in Figure 
45. Unless the displacement sensor is mounted across the test specimen, the 
rigidity of the mechanical load bearing frame to avoid measurement errors has 
to be very high relative to the test specimen stiffness. 
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Fig. 43 Payne Dynamic Tester (After Reference [5.11]) 
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Fig. 44 Takayanagl Dynamic Viscoelastometer (After Reference [5.35]) 
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Fig. 45 System Schematic of Closed-Loop Electrohydraulic 
Test Machine (After Reference [5.12]) 


108 











A similar machine, an electrodynamic tester, provides a more pure sinusoidal 
excitation than the electrohydraulie tester and was developed by Fort et al 
[5.37]], A sketch of the fixturing for this test machine is presented in 
Figure 46, 

5.4 Other Test Methods 


There are several additional test methods which do not fall into any of the 
classifications previously discussed. A number of these methods may be 
classified as wave propagation methods [5.15] . Wave propagation methods use 
the velocity and attenuation of waves travelling through an elastomer speci- 
men to determine its dynamic properties. These methods are generally restricted 
to relatively high frequencies (above 500 Hz - which is relatively high by 
rotordynamics standards) and to materials with moderate levels of damping. 

With low damping materials, standing waves occur; with high damping materials, 
precise measurements are difficult. Several devices have been designed which 
use longitudinal waves. For example, an apparatus used by Witte et al (j5.38j 
is shown in Figure 47, while one described by Morris et al [j>.39] is shown in 
Figure 48. Transverse vibrations have also been used in a device described 
by Ballon and Silverman [~5.40j, shown in Figure 49. Some methods have used 
brief vibration pulses rather than continuous wave trains. Such methods have 
been found to be adequate for measuring bulk modulus, but not shear modulus. 

Several test methods have been developed which utilize noncyclic displace- 
ments. The most common of these is the rebound resilience test which in- 
volves the measurement of the rebound of a mass falling on an elastomer 
specimen. Another is the rotary power loss machine (Figure 50), described by 
Bulgin and Hubbard [5.41^, in which an elastomer specimen is subjected to 
periodic but isolated deformations while rotating against a rigid wheel. 

There are also other classes of tests which are not strictly dynamic property 
tests but which relate, directly or indirectly, to dynamic properties. These 
include creep and stress relaxation tests, fatigue tests, vibration trans- 
missibility tests, and tests on cushioning materials (shock tests). These 
test classes are described in detail and several examples are presented by 
Payne and Scott [5.15]]. 

5.5 Base Excitation Resonant Mass Method 

The Base Excitation Resonant Mass (BERM) test method has been shown to be an 
easy to use, effective means of measuring the dynamic properties of elastomer 
specimens jj>.42j. The method provides for excellent control of the test con- 
ditions, while at the same time permitting maximum flexibility in the choice 
of these test conditions. Thus, with a sufficiently large shaker and a pro- 
perly designed test fixture, tests can be conducted for a variety of elastomer 
geometries and a broad range of frequencies, temperatures and strains. With 
relatively simple fixturing, tests can be run for frequencies in the range of 
100 to 1000 Hz. By using a specially designed test specimen, the testing of 
frequencies up to 2000 Hz should be achievable, and the use of support fix- 
tures for large resonant masses should permit tests to be run for frequencies 
as low as 20 Hz (5.42j . Tests have been run for elastomer specimens subjected 
to pure compression and pure shear, as well as for ring cartridge specimens 
and 0-rings [~5,43]]. Tests have also been run for peak-to-peak dynamic 
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Fig. 46 Electrodynamic Tester (After Reference [5.8]) 
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Fig. 47 Dynamic Test Apparatus Using Logitudinal Wave Propagation 
(After Reference [5.38]) 
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Fig. 48 Dynamic Test Apparatus Using Longitudinal 
Wave Propagation (After Reference [5.39]) 
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Fig, 49 Dynamic Test Apparatus Using Transverse Wave Propagation 
(After Reference [5.40]) 
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Fig. 50 Rotary Power Loss Machine (After Reference [5.41]) 
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strain in the range from 0*0005 to 0*08, and temperatures in the range of 
from 5-216°C (41-420°F) . A description of a typical BERM test rig is pre- 
sented, including the associated instrumentation and test procedures, in 
sufficient detail to be used as a guide in the performance of similar tests* 
Much of the design data presented in Chapter 6 is a result of BERM tests 
performed at MTI, 

The basic features of this test rig were developed by Chiang and Tessarzik 
[5.44J and others j5,28, 5*6, and 5.42j« The test rig was designed to impose 
precisely measured uniaxial vibration amplitudes upon a selected elastomer 
sample at desired test frequencies and under selected force preloads. The 
rig is capable of testing elastomer samples of virtually any size and shape 
over a wide frequency range, (The test sample cavity used is a cylinder 
approximately 12,7 cm (5 in) high by 12,7 cm (5 in*) in diameter,) Selection 
of the resonant mass to match elastomer sample properties permits testing at 
very high amplitudes at resonance with correspondingly lower amplitudes at 
off-resonance conditions. Thus, maximum test amplitudes are limited by 
elastomer characteristics and shaker input power rather than by the test 
rig itself. 

The vibration input to the test rig was obtained from a commercially available 
electromagnetic shaker system* capable of delivering 66,700 N (15,000 lb) 
force in the sinusoidal mode of vibration. The test fixture contained ambient 
temperature control facilities in the form of heating elements and cooling 
coils embedded in the sample base holder and flow paths for hot and cold air. 
Uncertainties due to temperature variations were minimized by accurately 
controlling the ambient temperature and by removing much of the heat generated 
in the specimen. Two different holding plate configurations provided shear 
and compression loading of the samples, A slight modification of the test 
rig permitted testing of ring cartridges, and a similar fixture was used for 
testing 0-rings, 

A schematic of the basic test rig can be seen in Figure 35 and a layout draw- 
ing of the same rig appears in Figure 51, Figure 52 shows the test rig 
mounted on the exciter table. 

The test rig was designed to meet the following functional requirements’, 

1, Very low parasitic damping, 

2, Inertia mass loading of the elastomer test sample 
with weights ranging from approximately 0,8 to 
227 kg (1,8 to 500 lb), 

3, Means for accurately controlling force preloading 
of the elastomer test sample, 

4, Temperature control of the ambient air within the 
test specimen cavity* 


* Ling Electronics, Model 335A Shaker with PP-35/70 VC Power Amplifier 
and SCO-100 Servo Control Center* 
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Fig. 52 Elastomer Test Rig Mounted on Vibration Table 
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The rig employs a unidirectional vibrating mass-spring system. In its 
simplest form, this concept might consist of an elastomer spring to which 
various amounts of inertia mass are bonded, and the system could be com- 
pletely free of parasitic damping or stiffness contributions from test rig 
elements. However, the practical requirements of balancing rather large 
inertia masses - up to 227 kg (500 lbs) - on top of the elastomer spring 
and the desired mechanism to provide elastomer sample preloading at all test 
frequencies introduce additional elements which have finite stiffnesses and 
small amounts of parasitic damping. Two radial guide bearings for large 
inertia masses and two air cylinders are required for preloading the elas- 
tomer samples. The stiffness and damping properties of these elements have 
an influence on the test rig response, but these properties can be indepen- 
dently measured and accounted for to minimize these effects on reported 
data. 

The resonant mass bonded to the top of the elastomer specimens consists of 
a number of elements. For high-frequency tests in which only a minimum of 
mass is required, the resonant mass can be reduced to the top plate of the 
specimen holder and the preload piston. The top plate is made of aluminum 
and the preload piston is made of titanium for minimum weight. Increased 
mass for medium- or low-frequency testing is provided by rigidly attaching 
a long rod to the top of the preload piston and by adding steel weights. 

The weights are centered by the rod and are axially restrained by spacers 
of various lengths and a lock nut near the upper end of the rod. The rod 
receives radial support from two n f rictionless" guide bearings; Figure 53 
shows one of the dismounted bearings. Each bearing consists of a hub and 
12 steel spokes. The hub is fitted over the end of the rod and is axially 
clamped. This kind of bearing arrangement provides good radial stiff- 
ness for reasonably high spoke tension but provides only limited freedom 
for axial motion if overstressing of the spokes is to be avoided. Conse- 
quently, changes in preload on the test specimen must be accompanied by 
lowering or raising of the outer frame, to which the guide bearings are 
attached. 

Preload of the elastomer test specimens is controlled by a pair of air 
cylinders. The lower cylinder is located directly above the test specimen. 
Its piston is attached to the specimen’s upper holding-plate and exerts a 
downward force on the specimen. The upper cylinder is located at the top 
of the test rig and its piston exerts an upward force on the resonant mass. 
When the resonant mass is large, the desired preload is usually achieved 
by pressurizing the lower cylinder and exerting a downward force in addition 
to the weight of the mass. The upper and lower preload pistons are sealed 
in their cylinders by two rolling diaphragm seals (Bellofram 3C-600-37-FPT 
and 3-119-119-CBJ) , The air cylinder inlet hole is 0.762 mm (0.030 in.) in 
diameter, making each cylinder essentially a closed cavity under vibration 
conditions. A constant preload force on the test specimen is maintained 
through regulation of air pressure in the cylinders. 

The ambient temperature within the elastomer test specimen cavity is con- 
trolled through the combined effect of both electric heating cables and 
water-cooled tubes that are installed in the base plate below the test speci- 
mens; Figure 54 is a photograph of the installed heating-cooling elements. 

The grooves in which the heating cable and the tubing are embedded have a 
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Fig* 53 Radial Guide Bearing 
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Fig. 54 Elastomer Test Sample Base with Installed Cooling Coils 
and Heating Cables 
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reverse taper to insure secure and rattle-free operation in a vibration 
environment. The electric heater is automatically controlled by a tempera- 
ture controller which obtains the necessary feedback from a thermocouple 
embedded in the metal of the base plate. Additional temperature control is 
obtained through the use of hot and cold air injected into the test rig 
cavity. 

For certain test conditions, not all available features of the test rig are 
needed. In particular, for test frequencies above 100 Hz, neither the mass 
support air spring nor the guide spokes (guide bearings) are required. The 
simplified configuration consists of a shake- table-mounted base plate to 
which the elastomer samples are bonded and a second plate bonded to the free 
side of the elastomer samples. A piston with Bellofram seals and various 
amounts of deadweight attached to the second plate vibrate freely on the 
elastomer spring when the base plate is subjected to sinusoidal vibrations 
from the shake table. 

The piston cavity is externally pressurized with air and serves to preload 
the elastomer sample in a static mode. An overall view of the simplified 
test rig preload cylinder is given in Figure 55, where a small mass is 
attached to the free end of the preload piston seen protruding from the 
cylinder. A certain amount of preload is needed to insure that the com- 
pression specimen is always in contact with the excited mass. Preload 
adjustment also permits the effect of preload on the dynamic characteristics 
of elastomer to be assessed. Preload is optional for shear, cartridge or 
0-ring configurations. For tests made without static preload required, 
removal of the seals and piston further simplifies the rig. 

Four test configurations have been tested at MTI using the BERM setup: 

• Shear-Type Samples 

• Compression-Type Samples 

• Cartridge-Type Samples 

• 0-rings. 

The shear test sample was made from four individual elastomer strips, each 
measuring 25.4 mm high, 48.8 mm long and 3.2 mm thick (1.0 x 1.92 x 0.125 
inch) . Each individual strip was glued with its large face to the side of 
a square aluminum block and with the opposite face glued to one additional 
aluminum holding block; Ren* plastic TWA-1077 adhesive has been used with 
excellent success. An assembled test element is shown in Figure 56. In 
the assembly of the test fixture, the four holding blocks were bolted to 
the shake table while the elastomer-supported center block was attached to, 
and therefore became part of, the resonant mass. 

The compression test sample consisted of 10 individual elastomer cylinders, 
each measuring 12.7 mm in diameter and 6.5 mm high (0.5 by 0.25 inch). These 
cylinders were arranged in a circular pattern and glued with their faces to 


* Ren Plastics, Lansing, Michigan 


121 




22 


Fig, 55 View of Shake-Table-Mounted BERM Test Rig with Preload Cylinder and Small Mass 








two aluminum plates; see Figure 57. The bottom plate was anchored to the 
shake table, and the top plate was connected to the resonant mass. 

Individual cartridge test specimens were made up of an annular elastomer 
body with rectangular cross section held between the cylindrical bore of a 
rectangular shell structure and a round shaft of smaller diameter. The 
elastomer ring was molded to the shaft and bonded to the surrounding shell, 
which was of split design to facilitate initial assembly and, later, to 
facilitate static preloading prior to testing. Figure 58 shows one cartridge 
specimen before and after assembly of the outer shell around the bonded 
elastomer ring on the shaft. 

A manufacturing drawing of the cartridge elastomer specimens is shown in 
Figure 59. The manufacturing process of the cartridge elements evolved 
from discussions with the manufacturer. Since a split shell was desirable 
so that positive radial preload could be applied for testing, molding of the 
elastomer to both shaft and outer shell was not considered feasible because 
tensile stresses, set up in the elastomer due to material shrinkage, would 
not be circumferentially symmetrical. In the selected manufacturing process, 
the elastomer was first molded on the shaft and then ground on the outer 
diameter to the dimension of the circular hole in the assembled shell with 
a small interference tolerance. The outer shell was then clamped around 
the elastomer and the assembly heat cured. Bonding of the elastomer to 
the shell was facilitated by a special coating applied to the bore of the 
shell prior to assembly; this coating is clearly visible in Figure 58. 

Prior to assembly into a complete test sample, all cartridge elements 
received a uniform preload of 5.3% of their uncompressed radial heights. 
Preloading was accomplished by replacing the shims, measuring 1.625 mm 
between the halves of the outer shell, with reduced- thickness shims, only 
1.120 mm thick. The assembly of the individual elastomer cartridges into 
the complete test sample is depicted in Figure 60. A 4-cartridge assembly 
was the largest that could effectively be packaged within the existing size 
constraints of the test rig. 

A similar test fixture was used for determining the dynamic properties of 
elastomer 0-rings. This test fixture is shown in Figure 61. The ambient 
temperature was maintained by means of an insulated thermal cover and an 
insulated base to enclose the test specimen. Both high-temperature heating 
air and controllable cooling air were provided to accurately control the 
temperature. 

The uncomplicated measurement requirements of the base excitation resonant 
mass testing method are one of its primary attributes. These measurement 
requirements are: 

1. Displacement measurement of elastomer support plate 
attached to the vibration table, relative to ground 

2. Displacement measurement of elastomer support plate 
attached to the resonant mass, relative to ground 
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Fig. 61 Assembly Drawing: O-Ring Test Fixture 
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Phase angle measurement between displacement measure- 
ments (1) and (2) above 

4. Displacement measurements between two elastomer support 
plates to insure the presence of required amplitude of 
elastomer motion 

5. Vibration frequency 

6. Temperature of elastomer support plates and ambient 
air temperature in test specimen cavity. 

For convenience, acceleration was measured at points (1) and (2) above 
instead of displacement, but, since the motions were sinusoidal, the dis- 
placements are directly obtained from the accelerations. The very low level 
of motion of the shaker body (from which the shaker table is isolated by air 
springs) was also measured by accelerometers. 

The input acceleration, point (1) , was measured by an accelerometer located 
on the support base for the elastomer test samples. Output acceleration on 
the resonant mass, point (2), was determined from three accelerometers, mounted 
directly to the top plate attached to the elastomer specimens. The use of 
three accelerometers, spaced 90° apart with their output signals displayed on 
one oscilloscope screen, permitted immediate detection of nonaxial motions of 
the resonant mass. When nonaxial motions of the resonant mass occur, they 
manifest themselves either as amplitude or phase angle variations among the 
three signals or as a combination of both. (Testing experience revealed that 
noticeable nonaxial motions occurred only in a few instances at certain test 
frequencies. Test data was generally discarded when nonaxial motions were 
observed.) For the 0-ring tests, two input and two output accelerometers 
were used. The displacement measurement, point (4), between the two elasto- 
mer support plates (relative amplitude across the elastomer) was accomplished 
with a noncontacting capacitance- type sensor. 

Chromel-alumel-type thermocouples were used to measure the temperature in the 
metal directly adjacent to one or two elastomer specimens in each assembly at 
a distance of approximately 1.6 mm from the elastomer surfaces. In the shear 
sample, the thermocouples were embedded at the midface location. In the com- 
pression test samples, the thermocouples were directly above and below the 
center point of the flat surfaces of the cylindrical specimens. In the 
cartridge specimens, one thermocouple was located in the axis of the shaft 
at the midpoint of the elastomer, and two additional thermocouples were placed 
diametrically opposite in the metal of the outer shell, approximately 1.6 mm 
from the outer diameter of the elastomer ring and also at the axial midpoint 
location of the elastomer ring. In the 0-ring fixture, the thermocouples were 
mounted in the inner and outer housings, immediately adjacent to each of the 
0-rings . 

Displacement and acceleration data were displayed on oscilloscopes in real- 
time form and were, together with temperatures from individual thermocouples, 
displayed in digital form for monitoring purposes. Amplitude signals from 
the four accelerometers and the one displacement sensor were sequentially 
switched by a minicomputer-controlled analog scanner into a two- channel 
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tracking filter, which provided visual readout of vibration frequency and 
two filtered amplitude signals. One input accelerometer signal was fed at 
all times into one of the two tracking filter channels and from there into 
a phase meter, where it served as a reference signal for the measurement of 
the phase angle of the output acceleration signals and the displacement 
signal. The dc values proportional to phase angle and amplitude from the 
phase meter and tracking filter were then covered intc digital form in 
digital voltmeters and read automatically by a PDP 11/34 minicomputer. The 
minicomputer acquired the data, performed some specified calculations, and 
printed the data and results of the calculations (including values for 
stiffness and damping) on a computer terminal used by the test rig operator. 

A schematic of the data acquisition system used at MTI* is presented in 
Figure 62. 

One significant feature of the BERM method is the performance of the 
scheduled tests at near-resonance conditions. Since phase angle between 
base excitation and resonant mass response is an accurate indicator of the 
amount of damping in the region of resonance, measurements should preferably 
be made in the phase angle range of 15-165° [5.28]. This generally requires 
that the test frequency be approximately 0.9 to Tt5 times the critical fre- 
quency range of 100-1000 Hz with sufficient test points, therefore, necessi- 
tates changes in the size of the resonant mass. Typically, between five and 
seven different mass arrangements ranging from 0.5 to 15 kg were found 
sufficient to cover the desired test range. The smallest mass consisted of 
the elastomer test sample top plate only (with no additional weight attached) 
Increasingly larger mass combinations were obtained by attaching weights of 
up to 15 kg directly to the top plate. In the case of the compression sample 
which was tested with preload applied through the air piston cylinder arrange 
ment, weights up to 10 kg were directly attached to the preload piston. For 
test frequencies in the range of 20 to 100 Hz, larger resonant masses were 
used with the support tripod (two to three different masses to cover the 
frequency range) . 

The following step-by-step process defines the procedure followed by the 
operator in executing the BERM tests for elastomer dynamic properties: 

1. Select and install a resonant mass. 

2. Statically preload the elastomer sample (for compression 
tests) . 

3. Enter information describing test conditions and resonant 
mass. 

4. Enclose the elastomer test sample and allow the tempera- 
ture control system time to adjust ambient temperature 
to the desired value. 


* Mechanical Technology Incorporated, Latham, N.Y. 
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Conduct frequency scans with low vibration levels applied 
to the base of the elastomer holding fixture until the 
approximate resonant frequency of the system is found. 

It may be noted here that for a base-excited, single- 
degree of freedom spring- damper-mass system, resonance 
occurs at an angle smaller than 90°, The deviation 
from 90° is essentially determined by the amount of 
damping in the system. 

6. Scrag the test specimen (run briefly at high amplitude) 
to eliminate the effect of structure (Section 4) and 
allow the specimen to reset for a few minutes before 
continuing. 

7. Maintain predetermined strain in the elastomer test 
sample by adjustment of the shaker power input level, 
and adjust the vibration frequency to obtain the nearest 
specified test frequency. Provided none of the accelera- 
tion and displacement signals show signs of abnormalities 
(distortions, or indications of nonaxial motion of the 
resonant mass), instruct the computer to acquire data. 

8. Instruct the computer to acquire data in the form of 
amplitude and phase for each sensor and temperature for 
each thermocouple. (The computer provides to the opera- 
tor an immediate calculation of stiffness, damping and 
power dissipation per unit volume of elastomer along with 
the raw data from the sensors and thermocouples.) 

9. Review these data and calculated results and indicate to 
the computer either that the data point is acceptable or 
not acceptable (normally it is acceptable) . 

10. If acceptable, store the data on a disk file; if not 
acceptable, discard the data point. 

11. Change the vibration frequency to the next nearest test 
frequency (with the phase angle between 15 and 165°) . 

12. Repeat steps (5) through (10) for each specified value 
of strain. 

13. Repeat tests, comprising steps (1) through (11), with 
each of the remaining masses in turn, each mass given a 
dynamic system with a different resonant frequency, per- 
mitting data to be taken at other test frequencies. 

14. Adjust the ambient temperature in the test sample cavity 
to the next specified value. 

15. Repeat steps (1) through (13) until data at all desired 
temperatures are obtained. 
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As mentioned, the test data is reduced by an on-line minicomputer. The 
primary objective of data reduction is to interpret these signals in terms 
of stiffness and damping. The structure of this data and the method of 
reduction can best be defined in relation to an analysis of a simple spring- 
mass-damper system subjected to base excitation, as shown in Figure 63. 


The equation of motion for the mass is 

rax 4* k.. (x - x ) + b(x - x ) = 0 (5-1) 

1 m o m o 

where x^ is displacement of mass 
x q is displacement of base. 

It is reasonable to assume that the motions of the base and mass are harmonic 
so that 
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and the equation of motion may be rewritten 
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Now, if <f> is the phase angle between x q and x^, then 
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Fig. 63 Model o! Base Excitation Resonant Mass System 
for Data Reduction Analysis 
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Thus, for known a, the transmissibility , or ratio of output amplitude |x I 
to input amplitude jx^ j , and (j>, the phase angle between input and output 
amplitudes, stiffness and damping (k^ and k2) can be evaluated directly; 
see equations 5-10 and 5-11. 


These expressions based on the simplest one-degree-of-f reedom model with no 
correction for tare damping or stiffness. This is in contrast to the ex- 
pressions used in Reference 5. 2 8 which make corrections for effects of 
damping in preload and load support pistons. The simplicity of the present 
treatment results from the simplicity and accuracy of the test procedure 
performed for the majority of the BERM tests, in which no external restraint 
of any kind was applied to the supporting mass. Only for the compression 
specimens is the preload piston employed. These expressions are, therefore, 
fully accurate for shear, cartridge, and 0-ring specimens and are subject to 
some uncertainty only in the case of the compression specimens. 


Statistical methods can be employed to develop empirical relations for test 
specimen or material dynamic properties as functions of the test parameters. 
Such empirical relations for several elastomer materials have been generated 
from BERM test results and are presented in Section 6.0. 

5 . 6 Summary 

Many test methods have been used for determining the dynamic properties of 
elastomers. These methods generally belong to one of three classifications: 
(1) forced vibration resonant test method, ( 2 ) free vibration resonant test 
method, and (3) forced vibration nonresonant test method. In general, the 
earliest methods which fall into the first classification were easy to use 
but suffered from the problem of severely limited test parameters. This 
limitation also applies to the methods which fall into the second classifi- 
cation. 


The more recent forced resonant methods, particularly the Base Excitation 
Resonant Mass (BERM) method, have used advanced instrumentation, including 
on-line computers, to eliminate these disadvantages. The result is a test 
method which provides significant flexibility in the choice of test para- 
meters, while requiring very straightforward operation and data reduction. 
The only major limitation of the BERM method is that it is not suitable for 
testing at very low frequencies (less than 20 Hz) due to the large resonant 


136 



masses required. However, for most rotordynamic applications, these 
frequencies are of minimal interest. (Fifty to 1000 Hz is generally of 
most interest and covers the operational characteristics of most high 
speed rotating equipment.) 

Forced nonresonant methods have also matured with the advent of advanced 
instrumentation. These methods, however, are more sensitive to measurement 
error (particularly phase measurement error) than the BERM method and are 
restricted to low frequencies (less than 100 Hz) and low dynamic strain 
which makes them generally limited for rotordynamics applications. 

Several other test methods have been used which do not fit into any of the 
three classifications. These include wave propagation methods and noncyclic 
displacement methods. Other classes of tests which are not strictly dynamic 
property tests, but relate, directly or indirectly, to dynamic properties 
include creep and stress relaxation tests, fatigue tests, vibration trans- 
missibility tests, and shock tests. 
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6.0 ELASTOMER MATERIAL PROPERTIES 


Access to reliable elastomer material properties is essential for the design 
of elastomer dampers. A collection of such properties for numerous elasto- 
mer materials is presented in this chapter. These material properties 
include both "general physical" properties and "dynamic" properties. The 
category of general physical properties includes properties such as density, 
environmentally related properties such as fluid compatibility, typical 
composition, static mechanical properties such as static shear modulus, and 
thermal, electrical and magnetic properties. Dynamic properties refer to 
the dynamic shear storage and loss moduli of the material and their variation 
with changes in operating conditions (e.g. frequency, strain, temperature, 
etc. ) . 

The general physical properties data presented in this chapter have been 
collected from a number of references; the dynamic property data presented 
have been determined empirically from results of Base Excitation Resonant 
Mass (BERM) tests with several elastomer materials. (The BERM test method 
was discussed in Section 5.0.) 

6.1 General Physical Properties 

Values for many general physical properties are presented for a number of 
commonly used elastomer materials. This collection of data is intended to 
be a summary for reader convenience. A more complete collection of data 
would constitute a handbook in itself, as is illustrated by Reference 6.1, 
which is the source for much of the data presented in this section. The 
reader is encouraged to refer to Reference 6.1 for additional information 
and for materials not included herein. Also, while most of the data pre- 
sented here is classified by generic elastomer type, Reference 6.1 also 
cross references many material properties to commercially available types 
(trade names) . 

Elastomer materials, including natural rubber, can be grouped into at least 
20 general classes. A brief description of each of the generic classifica- 
tion of elastomers, which includes such information as chemical structure 
and compounding procedure is presented in Reference 6.1. A summary of the 
general properties for 18 generic classes of elastomers is presented in 
Table 5 (6.2]*. 

Specific properties of several generic elastomer types are presented in 
Ref erence 6 . 3, some of which are reproduced here in Tables 6 through 12. 

Table 6 6.2 ! provides a broad classification of desirable and undesirable 

characteristics of the seven elastomers listed. Some general properties, 
including resistance to a number of environmental conditions, is presented 
in Table 7 . The hardness rating in Degrees British Standard (°BS) is 
identical to the International Rubber Hardness (IRHD) and is essentially 
equivalent to the commonly used rating of durometer (or Shore ® A hardness) . 

A direct relationship between elastomer hardness and static shear modulus has 


* Numbers in brackets indicate references found in Section 6.5. 

® Registered Trademark, Shore Instruments and Manufacturing Co., Inc., 
Jamaica, N.Y, 
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THE ELASTOMERS - DESIGNATIONS AND GENERAL PROPERTIES* 



■K 
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Reference 6.2 



Table 6 


SUMMARY OF ADVANTAGES AND DISADVANTAGES OF 
VARIOUS RUBBER TYPES* 


■■■ 

Useful Properties 

Properties that Limit 
Usefulness 

Natural 

High resilience, tensile strength 
and elongation; wide hardness 
range; long flex life 

Moderate to poor heat, 
oil, weather and ozone 
resistance 

Butadiene- 

styrene 

Better abrasion and general 
aging resistance than natural 
rubber 

Moderate to poor oil, 
tear, weather and ozone 
resistance 

Butyl 

High weather, ozone, acid and 
alkali resistance; very low 
permeability to gases 

Low tensile strength 
and resilience; high 
inflammability 

Neoprene 

High heat, aging, weather 
and flame resistance and flex 
life; moderate oil and chemi- 
cal resistance 

Only moderate in 
several important 
properties 

Nitrile 

Moderate to good heat and 
aging resistance; resistance 
to mineral oils, petroleum 
solvents, dilute acids and 
alkalis 

Relatively poor cold 
resistance 

Thiokol 

Maximum oil and solvent 
resistance 

High compression set; 
low tensile strength 
and elongation and 
heat resistance 

Silicone 

Maximum heat and cold 
resistance 

Low strength proper- 
ties; fairly high 
compression set; high 
cost 


* Reference 6.3 
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Table 7 


GENERAL PROPERTIES OF VULCANIZED RUBBERS* 



| 


Butadiene- 



I 

Nitrile 

ThiofeoJ ■ 




Natural 

styrene (a) 

Butyl 

Neoprene 

(b) 

(c) 

Silicone 

Hardness range, °BS 
Tensile strength (lb/in*) : 

* 

30-100 

35-100 

35-85 

40-90 

40-100 

40-90 

40-85 

‘pure gum’ 


>3000 

VL 

2000 

3000 

VL 

VL 

VL 

technical (d) 


800-4500 

800-4000 

800-2500 

800-3000 

800-2500 

600-1500 

400-1000 

Elasticity (e) 


H 

H 

FL 

M 

M 

L 

FL 

tan 5: at 20°C 


0-08 

0-094M4 

0*16 

0*09 

0-10-0*18 

— 

0*09 

(f) at 70°C 


0-06 

0-08-0*10 

0*08 

0*08 

0-09-0-11 

— 

0*08 

Resistance to compression set 


FH 

FH 

M 

FH 

M-FH 

L 

M 

Tear strength 


FH 

M 

M 

H 

M 

L 

L 

Abrasion resistance 


H 

H 

M 

H 

M 

L 

L 

Resilience (rebound): cold 


H 

M 

L 

FH 

M 

L-M 

H 

hot 


H 

M 

FH 

FH 

M 

L-M 

H 

Heat resistance 


FL 

FL 

M 

FH 

M-FH 

L 

VH 

Cold resistance 


H 

H 

H 

FH 

L-M 

FH 

VH 

T lt , °C (g) . . . 


-55 

-50 

-50 

-45 

-20 to 

-30 to 

-60 







-30 

—AS 


Working temperature range, d| 

C 

—55 to 

—45 to 

—50 to 

—20 to 

-20 to 

— 50 to 

— 90 to 



+70 

+ 100 

+ 125 

! +120(k) 

+ 120 

+ 95 

+250 (h) 

Ageing resistance, general 


M 

M-FH 

H 

H 

H 

H 

H 

Resistance to sunlight 


L 

L 

H 

H 

M 

M 

H 

Resistance to ozone and corona 

M 

L-M 

M-H 

FH-H 

L-M 

H 

H 

Resistance to flame 
Resistance to liquids: 

• 

L 

L 

L 

H 

L-M 

L. 

H 

fuels (aliphatic) 

. 

L 

L 

L 

M-H 

H 

H 

L 

fuels (aromatic) . 

• 

L 

L 

L 

L 

M 

H 

L 

Resistance to liquids: 







U 


mineral lubricating oils 


L 

L 

L 

M 

M-H 


L 

animal and vegetable oils 


Lr-M 

L-M 

H 

M-H 

H 

■ 

L-M 

water 


H. 

M-H 

FH-H 

MG) 

M 


M 

Electrical resistivity 


H 

H 

H 

M 

M 

■ | 

H 

Electrical power loss 


L 

L 

L 

M 

FH 

1 

L 

Electrical breakdown strength 


H 

H 

H 

M 

M 

■ t 

H 

Adhesion to metals 


H 

H 

M 

H 

H 

■ 

M 

Adhesion to fabrics 


H 

FH 

FH 

H 

FH 

mm 

M 


No Us. The level of each property is denoted as follows; VH: very high FH: fairly high FL: fairly low 

H: high M: medium L: low 

VL: very low 

(a) Usually known under the generic name GR-S or SBR. 

(b) Copolymers of butadiene and acrylonitrile in various proportions; with increasing amount of the latter the resistance to 
fuels and oils increases but cold resistance deteriorates, i.e. r„ becomes higher. 

(c) ‘Thiokol* is a trade name applied to a variety of synthetic rubbers; most arc thioplasts (Le. polymers of organic sulphides), 
and the information tabulated refers to these. 

(d) ‘Pure gum’ indicates a vulcanisate without fillers; ‘technical’ indicates the approximate range among the more usual technical 
compounds containing fillers. 

(e) In the popular sense, i.e. ability to stretch and recover. 

(f) For rubbers containing 45 parts (per 100 of polymer) carbon black (CK3), except the silicone rubber which contains a 
special filler; test frequency 16§ cps (Ecker, Kant. u. Gummi, 1956, 9, WT2). 

(g) The temperature at which elastic modulus is 10 times its value at 20 9 C. 

(h) No one type covers this whole temperature range; individual types vary from — 90/+150°C to *— 45/+250°G, 

(j) Must be specially compounded to resist water absorption. 

(k) Lower limit depends on type of neoprene, some types being more resistant than others to hardening by ‘crystallisation.’ 

Although the following types of synthetic rubber are not likely to be used for general engineering purposes, high cost being often 
one adverse factor, it may be worth noting their special advantages: 

Fluormated rubbers — excellent resistance to heat (operating range 0— 250°C), fuels, oils and corrosive chemicals. 
Polyacrylate rubbers — good heat and oil resistance; hence useful where oil-resistance of silicone rubbers is inadequate. 
Polyurethane (or isocyanate) rubbers — very high tensile strength and resistance to abrasion and tearing. 


* Reference 6.3 
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Table 8 


SPECIFIC HEATS, THERMAL CONDUCTIVITIES, AND SPECIFIC 
GRAVITIES OF RUBBERS AND COMPOUNDING MATERIALS* 



Specific 
heat (c) 
cal/g °C 

Conduc- 
tivity 
cgs units 

(b) 

Specific 

gravity 

Rubbers, unvulcanised or soft ‘pure gum’ 


xio- 5 


vulcanisates without filler; 




Natural rubber .... 

0-45-0*50 

32 

0*92-0*98 

Rutadiene-styrene (GR-S, SBR) (a) . 

0-44-0-48 

60 

0*94-1*00 

Butadiene-acrylonitrile (‘nitrile’) (a). 

0-47 

60 

0*94-1*10 

Neoprene ..... 

0*52 

45 

1*25-1*30 

Butyl 

0*44-0*46 

22 

0*92-0*98 

Barium sulphate (barytes) 

0*11 

80-250 

4*3 

China clay ..... 

0*21 

60-145 

2*6 

Carbcn black ..... 

0*20 

35-67 

1*8 

Iron oxide (ferric oxide) 

0*16 

185 

5*0 

Lithopone ..... 

0*115 

90 

3*6-4*2 

Magnesium carbonate 

0*30 

105 

2*2 

Whiting ...... 

0*20 

70-85 

2*7 

Zinc oxide ..... 

0*12 

140-165 

5*6 

Ebonite ...... 

0*33 (d) 

38-44 

1*2 (d) 


A 'ales, (a) Properties depend on ratio of butadiene to other constituent; values given 
are therefore approximate. 

(b) ~cal/cm 2 sec °G/cm; to convert to BTU/sqft hr °F/inch multiply by 2900. 

(c) Specific heat for most rubbers increases with temperature by about 
0*003-0*005 per °C; values quoted are for 20°-25°C. 

(d) natural rubber ebonite without fillers. 


Table 9 

LINEAR COEFFICIENT OF THERMAL EXPANSION* 



Coefficient 
per °C 

Rubbers, unvulcanised or soft ‘pure gum* vulcanisates, i.e. 
without filler; 

X10" B 

Natural rubber ........ 

22 

Butadiene-styrene (GR-S, SBR) ..... 

22 

Neoprene ......... 

20 

Butyl .......... 

19 

Silicone rubber ........ 

40 

Ebonite ......... 

7 

Mineral fillers (carbon black, whiting, zinc oxide) 

c. 0*5 
(i.e. negli- 
gible com- 
pared with 
rubbers) 

Aluminium ......... 

2*5 

Brass .......... 

2*0 

Steel .......... 

1*2 


Reference 6,3 
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Table 10 


VALUES OF THE 'CHARACTERISTIC TEMPERATURE' T and 

s 

GLASS-TRANSITION TEMPERATURE T FOR VARIOUS RUBBERS* 

8 



T a 


T s 



°K 

°G 

°G 


Natural rubber: 





unvulcanised . 

248 

-25 

— 71 


vulcanised (soft) 

251 to 260 (a) 

-22 to -13 

—70 to - 

-60 

,, (ebonite) 

360 

-1-87 

— 


,, (plasticised) (b) 

240 

-33 

— 


Butadiene-styrene rubbers (e) .* 





75:25 .... 

261 to 268 

-12 to -5 

—63 to - 

-54 

60:40 

283 

+ 10 

-1 


55:45 .... 

293 

-j-20 

-1 

U) 

50:50 .... 

296 

+23 

— | 

30:70 . 

328 

+55 

-J 

1 

75:25 (vulcanised) . 

75:25 (vulcanised, plasti- 

276 

+3 



cised) (d) . 

257 

-16 



Butadiene-acrylonitrile ('nitrile') 





rubbers, vulcanised (c) ; 
75:25 .... 

278 

+ 5 

—50 to 

— 40 

60:40 .... 

293 

+ 20 

—30 to 

—20 

50:50 .... 

310 

+37 

— 15 to +5 

Neoprene GN , vulcanised 

272 

-1 

-46 


Butyl rubber , unvulcanised . 

— 

— 

— 75 to 

—70 

Silicone rubber 

— 

— 

-123 



Notes, (a) Depends on the amount of combined sulphur; T t approximately equals 
the value for unvulcanised rubber plus 2’5°G for every 1 % sulphur. 

(b) 20 parts di-octyl sebacate per 100 rubber (Fletcher and Gent, Brit. J . 
Appl. Phys ., 1957, 8, 194). 

(c) The numbers following are the ratios of butadiene to styrene (or acrylo- 
nitrile) in the polymer; figures for butadiene-acrylonitrile rubbers have 
been calculated indirectly and are therefore approximate. 

(d) 40 parts ether plasticiser per 100 rubber (Fletcher and Gent, loc. cit.). 

(e) Published data from various sources differ too much to give reliable 
averages. 
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MECHANICAL PROPERTIES OF NATURAL RUBBER VULCANIZATES OF VARIOUS HARDNESS 
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Table 12 


SUMMARY OF EFFECTS OF VULCANIZATION AND COMPOUNDING ON 
DYNAMIC PROPERTIES* 


Property 

Region 

Vulcanisation 

Plasticizer 

Filler 

In-phase modulus 
(G* or E') 

rubbery 

Increase; 

’plateau 1 

lengthened 

Decrease 

Large increase; 
’plateau* leng- 
thened 


glassy 

Little effect 

Slight decrease 

Increase 

Out-of-phase 
modulus 
(G" or E") 

rubbery 

Increase 

Increase 

Increase 


glassy 

Little effect 

Increase 

Increase 

TAN <5 = 1/Q 

approx. ** 

rubbery 

Increase 

Decrease 

Increase 


glassy 

Decrease 

Decrease 

Increase 


transition 

i 

Peak lowered 
and broadened 

Peak lowered 
and broadened 
(at high plas- 
ticiser content 
peak may be 
raised again) 

Peak lowered 
and broadened 


* Reference 6*3 

** Q is the peak (maximum) transmissibility , a quantity sometimes used by 
engineers in place of tan 5 (see Section 2) 
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been demonstrated as shown in Figure 64 (from 6.4 ). Values of specific 
heat, thermal conductivity, specific gravity and linear coefficient of 
thermal expansion for several types of elastomers and compounding materials 
are presented in Tables 8 and 9. The characteristic and glass transition 
temperatures (as discussed in Section 4.0) for some typical mixes of several 
types of elastomers are presented in Table 10; static and low-frequency 
dynamic mechanical properties for natural rubber of various hardnesses are 
presented in Table 11; and the general effects of vulcanization and compound- 
ing are listed in Table 12. Specific lists of general properties for a 
large number of commercial compounds may be found in Reference 6.1 along 
with a listing of the relative resistance of twenty-two classes of elastomers 
to a number of classes of chemical agents. The resistance of a number of 
generic elastomer types to specific oils and solvents and to general environ- 
mental conditions are presented in Tables 13 and 14, respectively (from 6.2^). 

Due to an increased rate of oxidation, the use of elastomer materials at 
elevated temperatures has the effect of reducing the service life. This 
oxidation manifests itself as a decrease in tensile strength, elongation and 
resilience, and an increase in hardness. Eventually, cracks develop in the 
elastomer element and lead to failure. There is a threshold temperature 
below which the acceleration of aging is negligible. The high-temperature 
service limits for various elastomers are presented^ in the bar graph in 
Figure 65, 6.2 . For a reasonable service life, the elastomer temperature 

generally should be maintained well below the high temperature limit indicated 
on Figure 65. 

It has been shown 6.3 that the static compressive stress, f, is given by 

f = -G ftl-e) - (l-e)“ 2 ls (6-1) 

L_ 

where G is the shear modulus, e is the static compressive strain, and s is a 
shape function. For small e, equation (6-1) reduces to 

f - 3Gss (6-2) 

and thus, the term 3Gs becomes the effective compression modulus. Since 
elastomers are nearly incompressible (Poisson’s ratio is approximately 0.5), 
Young’s modulus may be approximated as three times the shear modulus. 

Therefore, the effective compression modulus determined above differs from 
Young’s modulus only by the factor of the shape function, s. The value of 
the shape function differs from unity because of the restriction of deforma- 
tion perpendicular to the load at the bonded surfaces; see Figure 66. Static 
shape functions for several typical geometries, given in Reference 6.3 and 
reproduced in Table 15, has been extended to include frequency-dependent 
dynamic shape functions 6.5 . This dynamic shape function (or factor) 
discussed in detail in Section 4.0, is also covered in Section 6.2. 

Compatibility with adhesives is another important property of elastomer 
materials as related to damper design. This compatibility is particularly 
important when using elastomer elements that are subjected to shear loading 
since the adhesive must carry the entire shear load. In general, elastomer 
elements may be bonded to a structure (usually metal) either by including 
the structure in the vulcanization process or by introducing an adhesive 
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Table 13 
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Aromatic Hydrocarbons 
Aniline, benzol , 


Carbon tetrachloride 
Lacquer solvents 
He thy I ethyl ketone 



ENVIRONMENTAL RESISTANCE OF ELASTOMERS* 
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200 300 400 500 600 


TEMPERATURE (F) 


NR = NATURAL RUBBER 

SBR = STYRENE-BUTADIENE 

CR = NEOPRENE OR CHLOROPRENE 

NBR = BUNA N OR ACRYLONITRILE-BUTADIENE 

MR = BUTYL OR ISOBUTYLENE ISOPRENE 

T = POLYSULFIDE 

CSM = CHLOROSULFONATED POLYETHYLENE 
VMQ = SILICONE FMQ = FLUOROSILICONE 




for Various Elastomers (After [6.2]) 


153 


EXPRESSIONS FOR SHAPE FUNCTION 
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side = ai 2 

much greater than a^ 1 + B(a2 - a^) z /h 

a^ almost as large as a 2 1-33 + D(a2 - ai) 2 /h 2 
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(usually epoxy) between the already vulcanized elastomer and the structure. 

In the former case* bonding agents (such as isocyanates or "chemloks" [6.4j) 
are often used for pretreating the bonding surfaces. In the latter case, 
good success has been achieved for several types of elastomers (particularly 
those for which dynamic properties are given below) with Renwald epoxy 
TWA-1077*. This adhesive has been effective even in an oily environment. 

In some cases (particularly when extreme environmental conditions are likely), 
it may be advisable to perform some tests on the compatibility and durability 
of a particular adhesive, or bonding method, with specific elastomer and 
structural materials, and a particular set of environmental conditions. 

When bonding an elastomer to a structure, it is advisable to clean the 
surfaces with emery paper and a compatible solvent (one which will not 
degrade the elastomer) , to apply the adhesive (usually to both surfaces) , to 
assemble the components, and to place and maintain pressure on the bond until 
the adhesive is fully cured. 

6*2 Dynamic Properties 

Empirical dynamic property data have been determined for several elastomer 
materials by the use of the BERM test method; (see Section 5.0). These 
materials were polybutadiene, flourocarbon (Viton) , nitrile (Buna-N) , EPDM 
(Ethylene-Propylene-Diene) and chloroprene (neoprene) . As dynamic properties 
of elastomers are sensitive to the percent and type of filler used in the 
material composition, batch-to-batch variation may occur in dynamic properties. 
For this reason, all elastomer specimens for each material tested were drawn 
from a single batch of elastomer stock. In each case, additional material 
was taken from the same batch for measuring some of the relevent general 
physical properties [6.6 and 6.7j presented in Table 16. Note that, of these 
materials, Viton has the highest transition temperature (closest to ambient). 
Consequently, the loss coefficient for Viton is the highest but is also the 
most sensitive to changes in frequency or temperature; see Section 4.0. 

Shear specimen BERM tests were conducted for each of these materials to develop 
relationships for the shear storage and loss moduli. These tests were conducted 
for a variety of operating conditions (i.e. frequency, temperature and dynamic 
strain). Using statistical regression analyses, the measured shear moduli have 
been quantitatively related to the values of the operating condition parameters* 
This relationship is described in this section. A set of compression sample 
tests (Section 5.0) has also been conducted for polybutadiene* Using the 
shear moduli values from the shear tests, values for the dynamic shape factors 
(Section 4.0) have been determined; these shape factors have been related to 
the operating condition parameters as discussed for the shear moduli* 

A series of BERM tests was also conducted forv 0-rings composed of fluorocarbon 
(Viton - CFM) and nitrile (Buna-N - NBR) . Each test for each elastomer 
involved material from a single batch, although not the same batch from which 
the corresponding elastomer material for the shear tests was drawn. 0-ring 
geometry is too complex to easily relate engineering parameters (i.e. stiffness 
and damping) to the shear moduli. Consequently , the test results presented are 
in terms of the support parameters as functions of the operating conditions. 


* Ren Plastics, Lansing, Michigan 
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SOME PHYSICAL PROPERTIES OF FIVE ELASTOMER MATERIALS USED FOR 
DYNAMIC PROPERTY TESTS AT MECHANICAL TECHNOLOGY INCORPORATED, LATHAM, NEW YORK 
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The results of shear and compression specimen tests for butadiene (Polybuta- 
diene - BR) are presented in detail; see [6.8] . As the form of the results 
for the other materials is similar, that data will be presented in somewhat 
less detail. For butadiene (Polybutadiene - BR) , typical plots of stiffness, 
damping and loss coefficient as functions of frequency (for a single tempera- 
ture but a variety of dynamic strain values) are presented in Figures 67 and 
68. Similarly, typical plots of stiffness, damping and loss coefficient as 
functions of dynamic strain (for a single temperature but for a variety of 
frequency values) are presented in Figures 69 and 70. All four plots are 
for the compression specimen tests at 66°C. However, the results for the 
shear specimen tests and the remainder of the compression specimen tests are 
qualitatively very similar. With frequency as an independent variable, a 
linear least squares fit has been applied to the common logarithms of the 
data. With strain as an independent variable, a second-order least squares 
fit has been applied to the common logarithms of the data. 

The effects of dynamic strain are summarized in Figures 71 and 72. The 
curves, with strain as an independent variable, show certain distinct 
characteristics for butadiene (Polybutadiene - BR) : 

• Stiffness decreases with increasing strain 

• Loss coefficient increases with increasing strain 

• Stiffness decreases consistently as the ambient tempera- 
ture is increased* 

• Loss coefficient decreases with increasing temperature* 

• Low strain loss coefficient falls by a small amount with 
increasing temperature* 

• Reduction in loss coefficient is more pronounced at high 
strain* 

• A distinct clustering of data points is evident about 
the fitted curves. 

• The plots for the compression test specimen are even 
more closely clustered about the fitted line than for 
the shear specimen 

• The plots as a function of strain provide, with reasonable 
scatter, a uniform trend for data at all frequencies. 

The plots as a function of frequency do not provide a uniform trend for the 
data at all strains; the scatter is very broad, indicating that, for this 
material and in this frequency range, the effects of strain are much more 
pronounced than the effects of frequency. In the case of stiffness, the 
data tend to show increasing stiffness with increasing frequency; but, as 
previously stated, the scatter is considerable. In the case of loss coeffi- 
cient, the data suggests little dependence on frequency. 


* Also see Figures 73 and 74. 
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STIFFNESS AND DAMPING (N/m) 



Fig. 67 Stiffness and Damping as a Function of Frequency for the 
66° C Butadiene (Poiybutadiene) Compression Tests 
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Fig. 68 Loss Coefficient as a Function of Frequency for the 66°C Butadiene (Polybutadiene) Compression Tests 



STIFFNESS 



Fig. 69 Stiffness and Damping as a Function of Strain for the 
66° C Butadiene (Polybutadiene) Compression Tests 
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Fig. 70 Loss Coefficient as a Function of Strain for the 

66°C Butadiene (Polybutadiene) Compression Tests 



STIFFNESS RATIO, K/K € (K € = 0.001) 



DYNAMIC STRAIN (DA) 


Fig. 71 Ratio of Stiffness at Finite Strain to Stiffness at 0.001 Strain as a Function of 

Dynamic Strain. Shear and Compression Specimens, Butadiene (Polybutadiene) 
at 32° C, 66° C, and 80° C. 
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Fig. 72 Loss Coefficient vs Dynamic Strain, Butadiene (Polybutadiene) 
Shear and Compression Specimens, 32° C, 66° C, and 80° C 
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Fig. 73 Effect of Ambient Temperature on Dynamic Properties of Shear Specimen 
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AMBIENT TEMPERATURE 





Based on the two sets of plots, dependence on strain is obviously more pro- 
nounced and more consistent for the range of frequencies investigated than 
the dependence on frequency for the range of strains investigated. 

Figure 75 presents plots of stiffness and damping as functions of frequency 
for a single low value of strain for' the butadiene (Polybutadiene - BR) shear 
test specimen at 66°C. These plots show a more consistent variation with fre- 
quency to which an acceptable straight line fit on log paper may be applied. 

When compared to the polybutadiene shear specimen results, the test data for 
the other four materials (fluorocarbon (Viton - CFM) , nitrile (Buna-N - NBR) , 
chloroprene (Neoprene - CR) , and EPDM) show more sensitivity to frequency and 
temperature. For example, using loss coefficient as an indicator, the effects 
of frequency at three different temperatures are compared for the five 
materials in Figures 76 through 78; similarly, the effect of temperature at 
a fixed frequency (500 Hz) is illustrated in Figure 79; and each of the five 
materials is considered individually in Figures 80 through 84. The following 
observations have been made from the results of these tests: 

• Fluorocarbon (Viton - CFM) is extremely sensitive to 
frequency and temperature (loss coefficient rises with 
increasing frequency and falls with increasing tempera- 
ture) . 

• Unlike the other four materials, loss coefficient for 
fluorocarbon (Viton - CFM) falls at high strain (may be 
related to self-heating and temperature sensitivity) . 

• Unlike the other four materials, the loss coefficient for 
EPDM falls with increasing frequency. 

• In general, butadiene (Polybutadiene - BR) is less sensi- 
tive to frequency and temperature, and more sensitive to 
strain than the other four materials. 

• At low temperature, the loss coefficient for fluorocarbon 
(Viton - CFM) (^0.7) is substantially higher and nitrile 
(Buna-N - NBR) (^0.15) somewhat lower than for the other 
materials . 

• At high temperature, 80°C, the loss coefficients for all 
the materials ranges between 0.1 and 0.2 for the entire 
frequency range. 

• Other than butadiene (Polybutadiene - BR) , nitrile 
(Buna-N - NBR) and EPDM are the least affected by 
temperature. 

Treating the distortion in the shear specimens as pure shear (discussed in 
Section 4.0), the shear storage and loss moduli are directly proportional to 
the measured stiffness and damping values, respectively, and may be calcu- 
lated from 
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Fig. 75 Stiffness and Damping vs Frequency. Low Strain (€ -0.0008), 
Shear Specimen, Butadiene (Polybutadiene) at 66° C 
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at 32° C 
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Fig. 77 Shear Specimen Data at 66° C 
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Fig. 78 Shear Specimen Data at 80° C 





Fig. 79 Temperature Variation of Elastomer Loss Coefficient (at 500 Hz) for Shear Tests 
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Fig. 80 Butadiene (Polybutadiene) Shear Specimen Data 
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Fig. 81 Fluorocarbon (Viton-70) Shear Specimen Data 
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g ' = k ; ( l } 


(6-3) 


G" = K" (f) (6-4) 

S A 

where 

2 

G 1 = shear storage modulus - N/m 

2 

G M = shear loss modulus - N/m 

= stiffness of shear specimen (real part) - N/m 

K" = damping of shear specimen (imaginary part) - N/m 
s 

t = thickness of shear elements - m 

2 

A * total bonded area (one side only) of shear elements - m 

Thus, the shear moduli (like the stiffness and damping) are functions of 
frequency, temperature and dynamic strain. The relationships between the 
shear moduli, the stiffness and the damping values of the compression specimen 
are much more complex. 


K T 

c 

3G'A 

t 

1+6' 

(— ) 2 
v 4t ; 

(6-5) 

K" 

c 

3G ,! A 

t 

1 + 6" 

(_^) 2 

'•4t ; 

(6-6) 

where 8 T and £" are shape 

factors 

(Section 4.0). 

These shape factors are 


also functions of frequency, temperature and dynamic strain. Thus, the 
logical procedure is to develop empirical relations for the shear moduli 
based on the results of the shear specimen tests and then to develop empirical 
relations for the shape factors, based on the results of the compression 
specimen tests together with the shear moduli relations. 

Section 4.0 also showed that the effects of temperature and frequency on the 
shear moduli are inversely related and may be represented by a single inde- 
pendent variable. Reduced shear moduli are used, as given by 

T 

c; = G’-f (6-7) 

T 

G^ = G"-| (6-8) 

where G^ and GJI are the reduced shear storage and loss moduli, T is the 
elastomer temperature and T c is the characteristic temperature (°K) for the 
particular elastomer materials. 

Assuming that dynamic strain is the only other important independent variable 
and that the effects of strain and reduced frequency can be separated, the 
reduced shear moduli are given by 
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G; = f f (cy>) g' (e) 


(6-9) 


= f" (cyo) g" (e) 


(6-10) 


where f and g represent unknown functions, (a^to) is the reduced frequency, 
and e is the dynamic strain. The coefficient a«j. is determined from: 


n a T 


C- (T - T ) 
1 c 


10 T C 0 + (T - T ) 
z c 


( 6 - 11 ) 


where the constants C ^ and C 2 depend on the elastomer material composition. 

In cases where high filler loadings are not present, C-^ and are essentially 

invariant with elastomer material and equal to - 8.86 and 101 . 6 , respectively 

[ 6 . 6 ], 

The logarithms of the shear moduli have been shown to vary approximately 
linearly with the logarithm of the frequency. Extending this relationship 
to the reduced frequency, the logarithm of the functions f f and f" are 
linearly related to the logarithm of the reduced frequency. Thus, the 
functions themselves are of the form 


f f (a T w) = a 1 (a T aj) 2 


( 6 - 12 ) 


f f, (a T w) * b^Ccyo) 2 


(6-13) 


where a and b are constants which are determined empirically for any particular 
elastomer material. 

The logarithms of the shear moduli also have been shown to possess a quadratic 
relationship to the logarithm of the dynamic strain. Thus, the logarithms of 
the functions g f and g ,f take the form 


log in g'(e) - a, + a, 10 g, n e + a (log n e)' 


log 10 g"(e) = b 3 + b 4 lo S 10 e + b 5^ log 10 e ^ 


(6-14) 

(6-15) 


to give 


g'(e) = a^e 


a 4 +(a 5 log 10 e) 


(6-16) 


b +(B lOg e) 

;"(e) = ble 


(6-17) 


where a 3 and b^ are the antilogarithms of a 3 and b^, respectively, i.e. 


a b 

10 , bi = 10 
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Equations (6-12) through (6-15) may be substituted into equations (6-9) and 
(6-10) to give the reduced shear moduli as 


G f 

r 


G" 

r 


( \ 2 4 
= a 6 (a T oj) e e 

^2 ^4 

= b 6 (a T aj) e e 


(a 5 log 10 e) . 

5 

(b 5 log 1() e) 


a 6 a l * a 3 
b 6 = b i* b 3 


(6-18) 

(6-19) 


The relations for the original shear moduli can be determined from equations 
(6-7) and (6-8) as 


G T = a.(— ) (a_aj) c e 


a 2 a 4 ( a 5 log 10 e) 


6 V T ' 
c 


G" » bg(— ) (cyj) £ £ 


^2 ^4 ^ 5 ^°^ 10 £ ^ 


( 6 - 20 ) 

( 6 - 21 ) 


Equations (6-20) and (6-21) can be made linear with respect to the unknown 
coefficients by taking the logarithm of both sides of ^ these equations to give 


log 10 (G V) = a 7 + a 2 (1 ° g 10° t T + + 


a 4 log 10 e + M lo Sin e) 5 a i = lo g-m a 


( 6 - 22 ) 


3 10 


10 6 


log 10 (G"— ) = by + b 2 (log 10 a T + log^w) + 


V° g lO E + b 5 (log 10 e) ; b 7 = log 10 b 6 


(6-23) 


Thus, once a quantity of experimental data has been generated for a particular 
elastomer material, the unknown coefficients of equations (6-22) and (6-23) 
can be found using a standard statistical regression analysis This has 

been done using the results of the shear specimen tests for butadiene 
(Polybutadiene - BR) , fluorocarbon (Viton - CFM) , nitrile (Buna-N - NBR) , 

EPDM and chloroprene (Neoprene - CR) . The values of these coefficients 
determined for these four materials are presented (with 90% confidence inter- 
vals) in Table 17, along with the estimate for the standard deviation of the 
data. That is, there is a 90% likelihood that the true value of each of the 
coefficients falls within the corresponding confidence interval. Using these 
coefficients, the shear moduli log terms can be estimated for any general set 
of operating conditions (for one of these four materials) using equations 
(6-22) and (6-23). 

A 90% confidence prediction interval can be assigned to each of these moduli 
as described in Appendix A to this Section. As an alternative to the above 
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REGRESSION COEFFICIENTS FOR DYNAMIC SHEAR STORAGE AND 
LOSS MODULI FROM DYNAMIC PROPERTY TESTS 
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equations, the charts in Figures 85 through 94 may be used for predicting 
the values of the shear storage and loss moduli for the five materials 
tested. For example, using the illustration in Appendix A of butadiene 
(Polybutadiene - BR) at 60°C and 400 Hz, with a dynamic strain of 0.01, 
gives a T w equal to 1.9 x 10“2 # Then, from Figures 85 and 86, the reduced 
moduli are 5.5 x 10^ and 9.5 x 10^ N/m^ for and G'-r, respectively. Multi- 
plying by T/T c (1.242) gives 6.8 x 10^ and 1.2 x 10& for G ! and G ff , 
respectively. These values compare favorably with the corresponding nominal 
values obtained directly from the equations in Appendix A. 


For the butadiene (Polybutadiene - BR) samples, the filler loading was found 
to be so high that the generally used values for and C 2 in equation (6-11) 
were not adequate and an interactive procedure was used to determine these 
values, as well as the value for T c [6.61. The filler loading was much 
lower for the remainder of the materials tested and the generally used values 
for and C 2 were found to be adequate. The values for T c for these materials 
were estimated from the values for T in Table 16. The appropriate values for 
C]_, C 2 and T c are given in Table 18 ror each of the elastomer materials tested. 


Given the relations for the shear modulus, as specified above, equations (6-5) 
and (6-6) can be solved for the shape factors 8 f and 8" to give 


6' = 

fK 't 
c 

n 

(— ) 

L3G f A 

± 

~ i 

3" = 

"k "t 
c 

— > 
1 

(— ) 

_3G'A 

5 


(6-24) 

(6-25) 


Then, using the results of compression specimen tests (such as those for 
butadiene (Polybutadiene - BR)), relations for 8 1 and 8" as functions of 
frequency, temperature and strain can be found using statistical regression 
analyses in the manner described above for the shear moduli. The relationships 
are again logarithmic and are given by 


log 10 B ’ = d l + d 2 log 10 u + d 3 log 10 + d 4 log 10 e (6_26) 

log 10 e " = e l + e 2 log l0 u + e 3 log 10 + e 4 log 10 e (6-27) 

The values of the coefficients in equations (6-26) and (6-27) for butadiene 
(Polybutadiene - BR) are given in Table 19, along with the values for 
standard deviation. Confidence prediction intervals can be calculated for 8 T 
and 8" as described in Appendix A. These equations may also be used to con- 
struct design charts, such as those in Figures 95 through 98. For example, 
using the conditions in Appendix A for butadiene (Polybutadiene - BR) 

(T = 60 C, 0 ) = 400 Hz, e = 0.01), Figures 95 and 97 gives values of 8 ! and 8” 
as 3.2 and 3.9, respectively. These values compare very favorably to the 
corresponding nominal values from Appendix A. 
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Fig. 86 Butadiene (Polybutadiene) Design Curves 
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Fig. 91 EPDM Design 
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Fig. 94 



Table 18 


VALUES FOR C^, C 2 AND CHARACTERISTIC TEMPERATURE (T ) 
FOR THE DYNAMIC TEST MATERIALS (j6.6 6 .l) 


Elastomer _ 

Material 

"i 

^2 

c (°K) 

butadiene 

(Polybutadiene - BR) 

-7.48* 

90.7* 

268.1* 

fluorocarbon 

(Viton - CFM) 

“8.86 

101.6 

306** 

nitrile 

(Buna-N - NBR) 

-8.86 

101.6 

293** 

chloroprene _ 

(Neoprene - CR) 

-8.86 

101.6 

272** 

EPDM 

-8.86 

101.6 

267** 


Determined by iterative approximation 
Estimated from transition temperature (Tg + 50) 
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Table 19 


REGRESSION COEFFICIENTS FOR DYNAMIC STORAGE AND LOSS SHAPE 
FACTOR FROM POLYBUTADIENE COMPRESSION TESTS 


Subscript 

1 

2 

3 

4 

Standard 

Deviation 

5 


For Storage Shape Factor (d) 
0.615 + 0.175 
-0.142 + 0.049 
2.13 +0.32 

-0.0837 + 0.0148 

0.0799 


For Loss Shape Factor (e) 
0.884 ± 0.145 
-0.180 + 0.041 
1.29 + 0.27 

- 0.101 + 0.012 

0.0662 
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Fig. 95 Shape Factor vs Strain 
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Fig. 96 Shape Factor vs Frequency 
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Fig. 98 Shape Factor vs Frequency 




The determination of the dynamic characteristics of 0-rings requires special 
consideration, 0-rings geometry and resulting loading conditions generally 
represent a much more complex dynamic support structure than do the simpler 
shear or compression buttons. Accordingly, 0-ring dynamic property testing 
requires that considerations of 0-ring material, assembly pre-distortions 
(stretch on squeeze), preload, strain, temperature and frequency all of which 
may greatly influence the dynamic property data. However, as 0-rings are often 
used as low cost dampers and consequently the results of these test are very 
meaningful. Such data has been provided in detail per Reference 6.10. The results 
are presented directly in terms of stiffness, damping and loss coefficient. 

The test plan for 0-rings was designed to provide maximum information for a 
minimum of test condition changes . With the large number of test parameters 
(7+ frequency), it was apparent that, unless an effort was made to hold the 
number of changes to a minimum, a very large test program would be required. 

To meet these constraints, a "parameter perturbation" approach was under- 
taken. Under this approach, a reference combination of all test parameters 
was defined and data for this reference combination was generated as a func- 
tion of frequency; the reference condition is defined in Table 20. One 
single parameter was then varied about its reference value and data generated 
as a function of frequency for each variation. This parameter was then set 
back to its reference value and a second parameter varied about its reference 
value. 

The variation was repeated for each parameter, in turn, according to the 
schedule of Table 21, for a total of nineteen tests. 

This parameter perturbation test method was advantageous in that it returns 
a large amount of information from a small number of tests. Its disadvantage 
was that it does not reveal the interactions which result from varying more 
than one parameter at a time. Thus, the results of this series of tests 
should be regarded as baseline data from which major trends and influential 
parameters can be identified; the less significant parameters can be identi- 
fied and eliminated from further tests in which interaction effects between 
parameters are investigated. 

The test frequency range generally lies between 70 and 1,000 Hz. A plot of the 
test data for the nominal test condition is presented in Figure 99. Actual 
test data points are shown; and a line which gives the minimum RMS deviation 
between test data and the fitted line is drawn through the full set of data 
points. All test results are for a pair of 0-rings, as tested, since this is 
the most likely configuration to be encountered in practice. 

Significant characteristics of most of the data are a strong increase in 
stiffness with frequency (typically by a factor of 2.5) and a high value of 
loss coefficient (in the range of 0.7 to 1.1). Of further interest is the 
fact that the 1,000-Hz stiffness is over 9 times the static stiffness of a 
pair of 0-rings (static stiffness measurement is discussed in the following 
paragraphs) . These characteristics are indicative of an elastomer in its 
transition region ;6.11j and are consistent with others measured for 
fluorocarbon (Viton - CFM) 6.12K 
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Table 20 


REFERENCE OR NOMINAL TEST CONDITION 


Material. . . . 

Viton-70 

Stretch 

5 percent 

Temperature . . 

25°C 

X-Section Diameter. . 

1/8 M Nominal (0.353 cm** 

Amplitude . . . 

7.62 x 10 ^ m 

Groove Width* .... 

135 percent 

Squeeze .... 

15 percent 

0-Ring 0D 

2 1/2" Nominal (6.35 cm) 


Based on actual X-Section diameter 
Measured Average 
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Table 21 


PERTURBATION PARAMETER VALUES 


Material 

Temperature 

Amplitude 

Squeeze 

Stretch 

X-Section Diameter 
Groove Width 


nitrile 
(Buna-N - NBR) 

38°C 66°C 

2.54 x 10 m 

5% 10% 

0 % 10 % 


fluorocarbon 
(Viton-90 - CFM) 

149°C 216°C 

1.27 x 10~ 4 m 

20% 30% 


1/16" nominal (1.778 x 10_? m) 

3/16” nominal (5,334 x 10 m) 

115/^ 150% of actual Cross-Section diameter 
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Fig. 99 O-Ring Test Results: Nominal Case 





In addition to the dynamic stiffness tests described above, the effect of 
a static load upon 0-ring deflection was investigated. Even under a static 
load the visco-elastic nature of the elastomer leads to time-dependent 
deflection; see Figure 100 where the deflection resulting from a sudden 
load application is plotted as a function of time. It is apparent that the 
deflection continuously increases with time and doubles in the elapsed time 
between 3 seconds and 300 seconds. Even after 7 minutes, the further 
deflection which occurs in the period between the 7 minutes and 8 minutes 
amounts to over 1% of the deflection at that point in time. 

In order to obtain some meaningful comparison between the elastomers, the 
deflection occurring 6 minutes after load application has been plotted 
against static load, and the results are presented in Figure 101. These 
curves illustrate some softening of the 0-rings as the static load is 
increased; for example, the 6 minute deflection of a nitrile (Buna-N - NBR) 
ring under 200 Newtons is five times the 50 Newtons for the same time. It 
is seen that nitrile (Buna-N - NBR) and fluorocarbon (Viton 70 - CFM) , which 
are both 70 durometer materials, have very similar static stiffness. The 
curves further show that the 90 durometer fluorocarbon (Viton - CFM) has 
close to 4 times the stiffness of the 70 durometer materials. 

Table 22 provides a summary of effective static stiffness values (load/ 
deflection) for the three materials under loads of 50, 100, and 200 Newtons. 
As previously discussed, the static stiffnesses are much lower than the 
dynamic stiffnesses. Table 23 provides the ratios of dynamic stiffness at 
100 Hz and 1,000 Hz to the effective static stiffness under 200 Newtons for 
the three different materials, and it may be seen that these ratios range 
from over 4 to almost 11. Ratios such as these can be used (when better 
data is not available) to estimate the dynamic stiffness of an elastomer 
specimen based on its static stiffness. A good estimate of the static shear 
modulus of an elastomer material can be made from hardness (durometer or 
Shore A®) as shown in Figure 64 [6.4). The static stiffness of the element 
can then be calculated from the geometry. 

Inspection of the test data plots for individual test conditions revealed 
that certain parameters had a significant effect upon 0-ring dynamic charac- 
teristics and that other parameters had a less pronounced effect. To bring 
these trends more sharply into focus, plots have been prepared in which 
power law curves for each of the seven parameter perturbation ranges are 
presented on a single plot. These "trend summary" plots are presented in 
Figures 102 through 108 and clearly reveal the important effects of the 
governing parameters. The effect of all seven parameters are discussed in 
the following section. 

6.3 Effects of Specific Parameters 

6.3.1 Material 


The selection of material and durometer value have a pronounced effect on the 
dynamic performance achieved from an 0-ring flexible support, as shown in 
Figure 102. Fluorocarbon (Viton 70 - CFM) and nitrile (Buna-N - NBR) have 
a similar average dynamic stiffness of 6 x 10^ N/m (39,000 Ib/inch) but that 
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Fig. 100 Typical Trace of Deflection vs Time Under Static Load for O-Ring Tests 


DEFLECTION, mm AT 6 MINUTES 



LOAD — NEWTONS 


Fig. 101 “Static” Deflection as a Function of Applied Load for O-Ring Tests 
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Table 22 


STATIC STIFFNESS SUMMARY FOR O-RING TESTS 


Load, 

fluorocarbon 

Effective Stiffness, N/m 

fluorocarbon 

Newtons 

(Viton-70 - CFM) 

nitrile 
(Buna-N - NBR) 

(Viton-90 - CFM) 
/* 

50 

1.52 x 10 6 

1.52 x 10^ 

7.69 x 10 

100 

1.29 x 10 6 

1.38 x 10 6 

5.71 x 10 6 

200 

1.14 x 10 6 

1.21 x 10 6 

4.55 x 10 6 


Table 23 

RATIO OF DYNAMIC STIFFNESS TO EFFECTIVE STATIC STIFFNESS 
UNDER 200 NEWTONS FOR 0-RING TESTS 


Material 



Stiffness Ratio 



100 Hz 

1000 Hz 

fluorocarbon (Viton 70 - 

CFM) 

4.35 

10.89 

nitrile (Buna-N - NBR) 


4.91 

6.98 

fluorocarbon (Viton-90 - 

CFM) 

4.73 

7.69 
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STIFFNESS (N/m) LOSS COEFFICIENT, 


1.5 




FREQUENCY (Hz) 


Fig. 102 Trend Summary: The Effect of Material on Stiffness 
and Loss Coefficient for O-Ring Dampers 
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STIFFNESS (N/m) LOSS COEFFICIENT 




FREQUENCY (Hz) 


Fig. 103 Trend Summary: The Effect of Temperature on Stiffness 
and Loss Coefficient for O-Ring Dampers 
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STIFFNESS (N/m) LOSS COEFFICIENT 




50 100 2 4 6 8 1000 

FREQUENCY (HZ) 


Fig. 104 Trend Summary: The Effect of Amplitude on Stiffness 
and Loss Coefficient for O-Ring Dampers 
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50 100 2 4 6 8 1000 

FREQUENCY (Hz) 


Fig. 105 Trend Summary: The Effect of Squeeze on Stiffness 
and Loss Coefficient for O-Ring Dampers 
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LOSS COEFFICIENT, 




Fig, 106 Trend Summary: The Effect of Stretch on Stiffness 
and Loss Coefficient for Fluorocarbon (Viton-70) 
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STIFFNESS (N/m) LOSS COEFFICIENT, 77 




Fig. 107 Trend Summary: The Effect of Cross-Section Diameter 
on Stiffness and Coefficient for Fluorocarbon (Viton-70) 
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STIFFNESS (N/m) LOSS COEFFICIENT, 




Fig. 108 Trend Summary: The Effect of Groove Width on Stiffness 
and Loss Coefficient for Nitrile (Buna-N) 
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fluorocarbon (Viton 70 - CFM) has a stronger frequency dependence over the 
range covered, and a loss coefficient almost three times that of nitrile 
(Buna-N - NBR) . Provided temperature and amplitude can be controlled, the 
fluorocarbon (Viton - CFM) appears to be a better damping material. If 
higher stiffness is required, an increase factor of 3 to 4 can be achieved 
by a shift to 90 durometer fluorocarbon (Viton - CFM); but the loss coeffi- 
cient for fluorocarbon (Viton 90 - CFM) is about half that of fluorocarbon 
(Viton 70 - CFM), and the frequency dependence of the fluorocarbon (Viton 90 - 
CFM) properties is less pronounced than those of fluorocarbon (Viton 70 - CFM) . 

6.3.2 Temperature 

The effects of temperature on stiffness and loss coefficient for fluorocarbon 
(Viton 70 - CFM) are very strong, particularly for small increases in tempera- 
ture above 25°C as shown in Figure 103. Both stiffness and loss coefficient 
fall sharply with increasing temperature; stiffness by a factor of two or 
more between 25 and 66°C; and loss coefficient by a factor of four. Above 
66°C the effects of temperature are less pronounced, and no change in stiff- 
ness was observed above 149°C. 

The trends in stiffness and loss coefficient as a function of temperature 
are emphasized by a cross plot of these quantities versus temperature, at 
200 and 800 Hz (Figure 109). Figure 110 compares the loss coefficient 
results with those from the shear test. These results indicate consistency. 

6.3.3 Amplitude 

The effect of amplitude on stiffness and loss coefficient is strong, as shown 
in Figure 104. In the range from 100 Hz to 1000 Hz, increasing amplitude 
consistently decreases stiffness and damping and decreases their frequency 
dependence. At 500 Hz, even a modest increase in amplitude from 7.62 x 10“^ mm 
(0.3 mil) to 25.4 x lCF^ mm (1 mil) cuts the stiffness by over 30%. An increase 
to 127 x 10*“ 3 mm (5 mils) causes an additional reduction of over 50%. The loss 
coefficient is less affected; an increase from 7.67 x 10”3 mm to 127 x 10”6 mm 
reduces the loss coefficient from 1 to 0.43. This effect is also consistent 
with the results of the fluorocarbon (Viton - CFM) shear specimen tests, 

6.3.4 Squeeze 

The radial precompression, or squeeze, imposed on the 0-ring affects its 
dynamic behavior, particularly its stiffness; see Figure 105. The Hertzian 
nature of the contact, in which increasing squeeze increases the area of 
contact between the ring and metal, thereby enables the force to be developed 
for a given relative radial displacement. For squeeze values in the region 
of 15% (Figure 111 defines squeeze) , the stiffness is least sensitive to 
squeeze (which probably makes 15% a good design value of squeeze) . As 
squeeze is reduced to 5%, the stiffness starts to fall off sharply, because, 
under even light gravity loads, contact between ring and metal is not maintained 
over 360° when 5% squeeze is coupled with the reference value of 5% stretch. 

The interaction of squeeze and stretch must be recognized. When the ring is 
stretched, its effective cross-sectional diameter is reduced and the effective 
squeeze is less than the squeeze as defined in Figure 111. 
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Fig. 109 Stiffness and Loss Coefficient vs Temperature for Fluorocarbon (Viton-70) 




LOSS COEFFICIENT 


1.2 



FREQUENCY (Hz) 

Fig. 110 Comparison of Fluorocarbon (Viton-70) Shear Specimen and O-Ring Test Data 
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When squeeze is increased to 30% , stiffness increases more rapidly than at 
15%. The reason could be contact between groove walls and ring or simply, 
the increased rectangularity of the squeezed 0-ring cross section and an 
increased influence of elastomer bulk modulus on deflection characteristics. 
The trends in stiffness as a function of squeeze are emphasized in Figure 
112 which is a cross plot of stiffness versus squeeze at 200 Hz. It shows 
the turning point at around 15% squeeze. 

6.3.5 Stretch 


The influence of stretch on stiffness and loss coefficient is shown in 
Figure 106. It is apparent that no clear trend can be defined since 5% 
stretch causes a higher stiffness and lower loss coefficient than either 
0 or 10%. Also, the spread for all stretch values is only 20% in stiffness 
and 10% in loss coefficient. Stretch is, therefore, shown to be of minimal 
importance. 

6.3.6 Cross Sectional Diameters 


As shown in Figure 107, cross sectional diameter is another parameter which 
develops no clear trend within the available test results. The reference 
value of 1/8 for cross section diameter develops the highest stiffness, 
and either a smaller or larger cross section causes a reduced stiffness. 

For loss coefficient, the spread of values is very small and some inter- 
section of the lines occurs. 

Simple analysis would suggest that cross sectional diameter should have 
little effect provided squeeze remains the same. Both stressed and strained 
dimensions of the cross section are proportional to cross section diameter, 
and stiffness should be independent of this diameter. Secondary effects 
such as stretch-squeeze interaction and the dependence of strain on diameter 
(for constant amplitude) may cause some secondary trends, which could be 
clearly identified by more extensive testing. 

6.3.7 Groove Width 


As would be expected, groove width has a negligible effect on stiffness and 
loss coefficient. This statement would have to be modified only if inter- 
action of an 0-ring with the groove walls occurred. Since the smallest 
groove width investigated was 1.15 times the cross sectional diameter, this 
interaction did not occur and, as shown in Figure 108, very similar stiff- 
ness and loss coefficient values were obtained for all three groove widths 
considered. 

The most notable feature of Figure 108 is that the stiffness and loss 
coefficient do not correspond to the reference values for other tests. This 
is because nitrile (Buna-N - NBR) 0-rings were inadvertantly used for these 
tests instead of fluorocarbon (Viton 70 - CFM) rings. The conclusion of a 
negligible effect of groove width is expected to be the same for fluorocarbon 
(Viton 70 - CFM). 

The recommended value for groove width (from 0-ring design guides) is 1.35d, 
shown to be an entirely satisfactory value to use. If space is at a premium, 
1,15 is also an acceptable value from dynamic performance considerations. 
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Fig. 112 Stiffness vs Squeeze for Fluorocarbon (Viton-70) 


6 . 4 Summary 


Accurate elastomer material properties are an essential element in the design 
of elastomer dampers. These property values are generally obtained from 
empirical data collected in a reference, such as this chapter, or determined 
from specific tests. The designer should be certain to use reliable elastomer 
property data. Reliable data will insure that the damper design is valid. 
Thus, it may be advisable to conduct some amount of independent experimental 
verification of the more important material properties. 
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APPENDIX 6A 


CALCULATION OF PREDICTION CONFIDENCE INTERVALS 
FOR DYNAMIC PROPERTIES OF ELASTOMERS 


When calculating the dynamic properties of elastomer materials from the 
equations in Section 6.0, it is useful to be able to estimate the uncertainty 
inherent in these calculations (due to the experimental error in the original 
test data). Thus, a 90% confidence prediction interval can be assigned to 
each of the modulus and shape factor values by attaching the term: 

+ 1.645 sA (6A-1) 


to the common logarithm of each of the reduced modulus and shape factor pre- 
dictions. In equation (6A-1) , s is the appropriate standard deviation from 
Table 17 or 18, H is a four element column vector given by: 


for reduced modulus: 


1 

lo8 10 a T + 1O8 10 U 
log 10 £ 
( lo gl0 £ ) 2 


(6A-2) 


for shape factor : 


A = 




1 

lo 8 io“ 

log 10 

- lo8 10 e - 


(6A-3) 


Q is a square matrix of rank four which is given for each of the materials 
for the modulus calculations in Table 6A-1 and for the shape factor calcula- 
tions in Table 6A-2. 


T T 

Note that the term jl Q where A is the transpose of A, is the same for 
either the storage or loss modulus of a given material. 


The center and extremes of the confidence interval for the actual moduli or 
shape factors can then be determined by taking the antilogarithm base 10 of 
each of the corresponding values. For the moduli, the reduced moduli are 
then converted to actual moduli by multiplying by the ratio of the elastomer 
temperature to the characteristic temperature. Alternately, the predicted 
values for the moduli or shape factors can be found directly from the equa- 
tions in exponential form, and the extremes of the confidence intervals found 
by multiplying by the antilogarithm of equation (6A-1) . 
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Table 6A-1 


Elements of Matrix from Dynamic Property Tests 


Elastomer Material (J Matrix 


butadiene 

0.478 

0.00260 

0.460 

0.103 

(Polybutadiene - BR) 

0.00260 

0.00788 

0.00445 

0.000805 


0.460 

0.00445 

0.461 

0.106 


0.103 

0.000805 

0.106 

0.0252 

fluorocarbon 

0.225 

-0.00314 

0.162 

0.0273 

(Viton - CFM) 

-0.00314 

0.00336 

0.00163 

0.000138 


0.162 

0.00163 

0.129 

0.0229 


0.0273 

0.000138 

0.0229 

0.00436 

nitrile 

0.365 

-0.00290 

0.292 

0.0550 

(Buna-N - NBR) 

-0.00290 

0.00405 

0.000988 

0.0000368 


0.292 

0.000988 

0.250 

0.0490 


0.0550 

0.0000368 

0.0490 

0.0100 

chloroprene 

0.650 

0.0500 

0.559 

0.114 

(Neoprene - CR) 

0.0500 

0.0575 

0.0243 

0.00292 


0.559 

0.0243 

0.507 

0.107 


0.114 

0.00292 

0.107 

0.0235 

EPDM 

0.505 

0.00516 

0.438 

0.0886 


0.00516 

0.00822 

0.00496 

0.000947 


0.438 

0.00496 

0.396 

0.0826 


0.0886 

0.000947 

0.0826 

0.0177 


222 









Table 6A-2 


Elements of (^Matrix from butadiene (Polybutadiene - BR) 

Compression Tests 


1.78 

-0.480 

•0.480 

0.140 

1.60 

0.280 

0.0223 

0.00326 


-1.60 0.0223 

0.280 0.00326 
6.07 -0.0406 

-0.0406 0.0128 
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For example, for butadiene (Poly butadiene - BR) at a temperature of 60 degrees 
centigrade, a frequency of 400 HZ, and a dynamic strain of 0.01, the logarithms 
of the reduced moduli become, from equations (6-22) and (6-23) 

log 1Q (0.805 G* ) = 6.76 (6A-4) 


log 10 (0.805 G") = 5.97 
The vector Z, for modulus calculations, becomes 


1.0 


Z 


0.280 


- 2.0 


4.0 


(6A-5) 


(6A-6) 


Equation (6A-1) then gives 

+ 0.044 for storage modulus (6A-7) 

+ 0.090 for loss modulus (6A-8) 

Thus, the extremes of the 90% confidence intervals for G T and G" are 

6.46 x 10 6 N/m 2 < G' <7.91 x 10 6 N/m 2 (6A-9) 

0.95 x 10 6 N/m 2 < G" <1.44 x 10 6 N/m 2 (6A-10) 

while the nominal values for G T and G n are 

G' - 7.15 x 10 6 N/m 2 (6A-11) 

G" = 1.17 x 10 6 N/m 2 (6A-12) 

Similarly, for the shape factor calculations, for the same specifications 

log 10 B’ = 0.500 (6A-13) 


From equation (6A-3) 


log 10 e" = 0.595 


1.0 



£ = 


3.40 

0.094 


- 2.0 


(6A-14) 


(6A-15) 
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Equation (6A-1) then gives 


+ 0.131 for storage shape factor (6A-16) 

+ 0.109 for loss shape factor (6A-17) 

Thus, the extremes of the 90% confidence intervals for 6' and 3" are 

2.34 < 3' 4.28 (6A-18) 

3.07 < 3" 5.06 (6A-19) 

while the nominal values for 3' and 3" are 

3' = 3.16 (6A-20) 

3" - 3.94 (6A-21) 


The nominal values for G’ , G", 3' and 3" can also be easily extracted from 
the design charts in Section 6.0. 
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7.0 PRACTICAL DESIGN CONSIDERATIONS AND PROCEDURES 


This section provides a detailed discussion of the steps required for the 
proper design of elastomer dampers and the practical considerations implicit 
in such a design. A brief discussion related to applications for control of 
unidirectional vibrations is presented; however, the major emphasis is 
placed on rotordynamic applications, the subject of this handbook. The 
specific procedure required for the design of elastomer dampers for rotating 
machinery is described in detail, including both machinery dynamics and 
hardware-related considerations. Several examples of actual designs, in- 
cluding details of the design procedure, designed hardware, and test results, 
are presented in Section 8.0. 

The basic steps involved in elastomer damper design are to: 

e Analyze the system dynamics and establish the best or optimal 
levels of support damping and stiffness to be provided by the 
elastomer elements under the given design requirements 

9 Choose the elastomer material and geometry to best approximate 
this stiffness and damping, based on reasonable ranges of the 
operating parameters (e.g. frequency, temperature and strain) 
and any appropriate environmental considerations (e.g., fluid 
immersion, ozone) 

• Predict the system response with the actual expected ranges 
of elastomer stiffness and damping 

« Design the actual damper hardware based on manufacturing, 
assembly, and other practical considerations. 

To assist the engineer, each step is considered in detail in this section. 

7.1 Determine Best or Optimal Damping and Stiffness 

For unidirectional, as well as rotordynamic applications, it is first necessary 
to develop an analytical model which approximately represents the dynamic 
characteristics of the vibratory system. For simplicity, linear models are 
often used and are generally adequate. For unidirectional applications, 
elastomer dampers are often used in the form of foot pads on which relatively 
rigid machines are mounted. In such cases, a single-degree-of-freedom, damped, 
forced vibration model, as described in Section 2.0, may be sufficient. If 
the most likely forcing frequencies are known (from previous analysis or 
machinery experience) , it is a straightforward procedure to determine the 
stiffness and damping levels which minimize the transmissibility of the damper 
[7.1]*. 

More complex unidirectional systems, such as those employing constrained-layer 
damping [7.2J, may require fairly complicated analytical models. In such cases, 
numerical tools using discretized models (e.g., finite difference or finite 
element representations) and digital computers are generally used l 7 : 3 J- 
Although these predictive tools are more expensive than those approximations 
for the single-degree-of-freedom model, the objective of the analysis is the 

* Numbers in brackets indicate References found in Section 7.6 
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same: to minimize the system vibration (using some predefined critical 

response criteria) through variation of the damping and stiffness of 
dampers. Essentially, this optimization procedure involves analyzing the 
system model with a number of realistic combinations of damping and stiffness 
for dampers. The optimum combination of stiffness and damping can then be 
selected. Awareness of the effect on the system response of small changes in 
these parameters is essential, since a certain amount of uncertainty exists 
in the specification of elastomer properties (as there is in the design of 
any damper). Further, elastomer materials tend to be moderately sensitive 
to changes in operating conditions and to be somewhat nonlinear. If there 
are alternatives available in the choice of damper location, it would be 
advisable to locate the damper (s) as near to the antinode (s) of vibration 
as is practical, or at any location which would insure that motion of the 
damper would maximize the energy dissipated. 

The first step in the analysis of rot or dynamic systems also involves the 
development of an analytical model. The amount of detail required in such 
a model is dependent on the complexity of the rotordynamic system. For 
systems with rigid rotors, for example, it is often sufficient to lump the 
mass and inertia of the rotor and to consider only the one or two degrees-of- 
freedom provided by the rotor supports. However, more detailed rotor models 
are necessary for flexible rotors. Rigid rotors are those for which no 
significant shaft deformation is exhibited for running speeds up to the 
maximum operating speed of the rotor. Conversely, flexible rotors are those 
which exhibit appreciable shaft deformation. Numerous computer software 
packages are available for analyzing rotordynamic systems. The actual system 
modelling procedure is dependent on the particular type of analysis used. 

The most common of these analytical tools employs transfer matrix [7.4, 7.5, 
7.6] or finite element ^7.7j representations. These representations are 
generally composed of linear elements , although nonlinear analyses may also 
be conducted. Several examples involving the use of such linear analytical 
tools are presented in Section 8.0. 

Selection of the most appropriate analytical tools require that the principal 
motivation for the addition of damping to the system be ascertained. For 
example, the damping may be necessary to aid in the control of synchronous 
vibrations due to rotor mass unbalance or to stabilize otherwise unstable 
nonsynchronous, or even transient, vibrations. Some form of acceptable 
vibration or stability criteria must then be established as a guide in the 
design process. 

Part of the modelling procedure necessitates the determination of damper 
location (s). Rotor dampers are generally located in the bearing supports and, 
accordingly, may also be used to support the superstructure of the rotordynamic 
system. Several alternate locations can be available for the dampers, with 
the final decision being based on the analytical results and other considera- 
tions (temperature, maintainability, etc.). 

To begin this analytical procedure, realistic ranges of stiffness and damping 
values (and possibly damper locations) are selected. The rotordynamic model 
is then analyzed for support parameters spanning these ranges and a pre- 
established critical response (based on the criteria mentioned above) is 
plotted as a function of these support parameters. The optimum values of 
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stiffness and damping can then be identified. It is also useful to note the 
effect of small changes in support stiffness or damping on the critical 
response. A commonly used critical response criterion is system logarithmic 
decrement 3 which gives an indication of overall damping effectiveness of 
each critical speed, as discussed in Section 2.0. Thus, the log decrement 
is indicative of the sensitivity of the rotor response to unbalance, as well 
as the system stability. Section 8.0 discusses several examples in which 
log decrement is used as a design criteria. 

7.2 Specify Elastomer Material and Geometry 

Proper specification of the elastomer material and geometry is dependent on 
a number of factors. Material selection is affected by both dynamic and 
environmental considerations. The desired dynamic properties of the damper, 
particularly loss coefficient, have a major impact on the choice of elastomer 
material. The damper operating parameters (frequency, temperature and strain) 
are also important. 

Elastomer dampers are often subjected to severe environmental conditions, such 
as exposure to fluids or gases, which may be incompatible with the elastomer 
material. Even in a relatively clean environment, exposure to ozone and ultra- 
violet light (at normal ambient levels) over an extended period of time can 
result in aging of the elastomer material. Such aging is generally manifested 
by an increase in the stiffness of the elastomer. The aging effect may be so 
significant for some elastomer materials that a limitation must be placed on 
the shelf life of the elastomer specimens (generally not less than several 
years) . 

Data presented in Section 6.0 concerns the resistance of many elastomer 
materials to degradation from prolonged exposure to a variety of liquids and 
gases. Using such data, in cases when environmental considerations are 
expected to be important, elastomer material selection can be narrowed to 
those best suited to a particular set of environmental conditions. Similarly, 
the effects of frequency, temperature and dynamic strain on the dynamic pro- 
perties of an elastomer material can be evaluated using empirical relation- 
ships (Section 6.0). 

Once the appropriate elastomer material has been selected, a suitable mode 
of bonding, using an adhesive or bonding agent, must be chosen. A brief 
discussion on bonding of elastomer materials is also included in Section 
6.0. If elevated elastomer temperatures are likely to be a problem, 
cooling, using either forced air convection or liquid in cooling channels, 
should be considered. 

There are many possible geometric configurations for the elastomer elements 
in a damper. Four typical examples are: 1) Shear cartridge, 2) Compression 

button cartridge, 3) Continuous compression cartridge and 4) 0-rings. The 
particular choice of geometric configuration is determined somewhat by the 
application. However, the shear cartridge and compression button cartridge 
configuration generally are recommended since they inherently lend themselves 
to a more straightforward evaluation of dynamic properties and to less com- 
plicated damper hardware design, assembly and maintenance. 
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A shear cartridge configuration can be in the form of a continuous cartridge 
(Figure 113) or a segmented cartridge (Figure 114) . In this segmented shear 
cartridge design, the elastomer elements are subjected to pure shear loading 
from rotating radial forces. Using a slightly different design (Figure 115), 
the elastomer elements are subjected to a combination of shear and compressive 
loading. The pure shear design simplifies the analysis required to predict 
dynamic properties. Any elastomer element which is subject to an oscillatory 
compression (and tension) load must be sufficiently preloaded in compression 
to preclude the possibility of subjecting the element to an absolute tensile 
load; elastomer materials and bonding interfaces tend to be less reliable In 
tension than in compression (or shear). 


For the pure shear designs (either continuous or segmented) , the calculation 
of the stiffness and damping values for the damper is straightforward. Pre- 
sumably, the shear storage and loss moduli for the chosen elastomer material 
are available from test results or from data as presented in Section 6.0. 

The stiffness and damping of the damper are then given by 


K’ = G ' 
s 


A 

t 


(7-1) 


K" = G" 
s 


A 

t 


(7-2) 


where 

K = K T + iK ft 
s s s 

and 

= damper shear stiffness (real part) - N/m 

K" = damper shear damping (imaginary part of stiffness) 
s 

2 

G T = shear storage modulus - N/m 

2 

G" = shear loss modulus - N/m 

A = total bonded area (one side only) of shear elements 
t = thickness of shear elements - m 

For example, if a segmented pure shear cartridge is constructed with three 
cylindrical buttons, 10 mm in diameter and 3 mm thick, the total area, A, 
would be 2.36 x 10~^ m^, and the ratio of stiffness component (K T or K") to 
shear modulus (G ? or G") would be 0.0785 m. The analysis of the combined 
shear and compression configuration. Figure 115, is somewhat more complicated 
and is very similar to the analysis for the compression button cartridge 
discussed in the following paragraphs. 

A typical compression button cartridge configuration is illustrated in 
Figure 116. The application of a rotating force results in the elastomer 
elements being subjected to a combination of compression and shear loading. 
Also, as discussed, preloading the elastomer buttons in compression is 


- N/m 


2 

- m 
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Fig. 113 Example of Continuous Shear Cartridge Elastomer Damper Configuration 
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necessary to exceed the largest anticipated dynamic tensile load. The 
principal attraction of this configuration is the simplicity it lends to 
the damper hardware design; see two examples in Section 8.0. 

When calculating the effective stiffness and damping of this elastomer damper 
conf iguration, it is convenient to separate the effects of shear and compres- 
sion loading and to superimpose the separate results. The stiffness and 
damping of an individual button in shear are calculated directly from equations 
(7-1) and (7-2). However, the relationships between the shear moduli and the 
compression stiffness and damping of an elastomer button are considerably more 
complex; this complexity is due to the resistance to out-of-plane deformation 
provided by the bonded surfaces (Section 4.0). These relationships are of the 
form 

K C = (1 + e,<s2) (7-3) 

K c " (1 + 3 " 6 } (7 ~ 4) 

where the moduli of elasticity are taken to be three times the shear moduli 
(Poisson’s ratio of about 0.5), the 3 f and 3" are dynamic shape factors, and 
5 is a dimensionless geometric factor (D/4t for a short cylinder) . Values 
for the 3 ? s (as functions of the operating conditions) may be obtained from 
empirical data such as that presented in Section 6.0 or directly from experi- 
mental results. 

Once the stiffness and damping for an individual button have been determined 
separately for compression and shear loading, the overall stiffness and 
damping of the elastomer damper may be calculated {j.8j . The expressions K 
and K g are used to represent the complex stiffness of a single button in C 
compression and shear, respectively, where 

K = K f + iK" (7-5) 

s s s ' 

K = K 1 + iK" (7-6) 

c c c 

K = the elastomer button stiffness in shear 
s 

K^ = the elastomer button stiffness in compression 

Consider the support geometry given in Figure 117a which shows three elastomer 
buttons, spaced at 120° intervals and separating a bearing housing from the 
surrounding pedestal. When the housing is given a downward displacement of 
some distance X, the individual buttons are deformed in some combination of 
compression and shear. This is illustrated in Figure 117b for one of the 
buttons, where: X is the total downward deflection, X c is that component of 

the total deflection taken up by elastomer compression, and X s is that com- 
ponent accommodated by shearing. If the individual buttons are labeled 1, 2, 
and 3 starting from the top, the compression and shear deflection of the 
individual buttons can be described in terms of the imposed bearing housing 
deflection X. The deflection of the first button is simply as follows: 


2 34 
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Support Geometry and Button Deformation for Stiffness Calculations of 
Compression Button Cartridge Elastomer Damper Configuration (after [7.8]) 



(7—7) 


X 


lc 


= X 


X ls - 0 (7-8) 

where X lc is the compressive, or normal, deflection and X ls is the shear 
deflection. For buttons 2 and 3, the corresponding deflections are 

X 2c ' X 3c ‘ Xc ° s60 ° * I < 7 “ 9 > 

X 2s = X 3s = Xcos30 ° 2 X (7-10) 

The forces produced by these displacements are the products of the individual 
button spring constants (K c for compression and K g for shear) with the normal 
or shear displacements. For the first button, the forces are 

F lc - K c X (7-11) 





F = 0 
Is 

(7-12) 

where F^_ c 
buttons 2 

and 

and 

F^ s are the 
3, the same 

compressive and shear forces, respectively, 
procedure is followed: 

For 




F„ = F_ = K f 

(7-13) 




F 2s " F 3s ‘ Vi X 

(7-14) 


The total force in the direction of displacement X is the sum of the force 
through button 1 and the components of force in the X direction from buttons 
2 and 3. 


F T ■ F lc + <F 2c + F 3c> “ s 60 ° + (F 2s + r 33> cos 30 ° 


K X + f 
c 2 


K 

c 



(7-15) 


The above equation may be simplified to give 

F^ * 1.5 (K + K )X (7-16) 

T c s 

If the elastomer mount consists of cartridges with a number of individual 
axial elastomer buttons or elements, the total force is, simply, the sum of 
the individual element forces, 

F T = 1.5 N b (K c + K g ) X (7-17) 

where Nb is the number of buttons or elements per cartridge at each circum- 
ferential location. The radial stiffness of the assembly, K^, is the 
resultant force divided by the imposed displacement. 

h m r ' 1 - 5 h b « c + V (7 - 18) 
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If the circumferential arrangement of elastomer buttons is different from 
that shown in Figure 117, the above analysis must be modified to account 
for the actual configuration. While the final equation will be somewhat 
different, the general procedure will be essentially the same. 

Some elastomer damper designs utilize continuous compression cartridges 
(Figure 118) or 0-rings (Figure 119) as the elastomer elements. 0-rings 
provide certain conveniences as they are often used as seals (and may thus 
serve a dual purpose) and may be obtained quickly and cheaply in a large 
variety of sizes and materials. However, the mode of dynamic deformation 
of continuous compression cartridges and 0-rings is very complicated, and 
attempts at representing such deformation analytically have not been very 
successful. Consequently, instead of using analytical procedures for 
predicting stiffness and damping, it has been generally necessary to generate 
large quantities of experimental data (as presented in Section 6.0) directly 
in terms of stiffness and damping. Because 0-rings do not have a constant 
cross-sectional dimension (relative to radial position), they are more 
sensitive to changes in preload and more nonlinear than for the other 
geometric configurations. In fact, the effect of modest preload (10% or less) 
on the properties of pure shear and compression specimens has been shown to 
be relatively unimportant 7.9~], while that for 0-rings (in the form of 
squeeze) has been shown to^be very significant 7.10 i. Consequently, the 
use of such elastomer elements in a rotor damper' design requires the avail- 
ability or generation of test data for the specific geometry. 

When considering elastomer damper geometry, it is important to keep in mind 
the limitations within which the analyses are valid. The elastomer element 
must be short enough that beam effects become negligible when considering 
shear and that column effects become negligible when considerating compres- 
sion (Section 4.0). 

7.3 Reanalyze System with Specified Elastomer Damper Configuration 

Once the elastomer material and geometric configuration have been specified, 
the corresponding support stiffness and damping values can be determined; 
these values can then be compared with the optimal stiffness and damping 
values. Invariably, there will be some difference, although it may not be 
large. Consequently, reanalysis of the system dynamics will generally be 
necessary to consider these actual support parameter values. This is true 
for both unidirectional and rotordynamic systems. When analyzing rotordynamic 
systems with elastomer dampers, consideration of the effects of angular 
stiffness and damping in the dampers may be important, particularly when 
axial rows of elastomer shear or compression elements are used. (See Section 
8.0 for an example.) 

7.4 Design of Damper Hardware 

Many factors must be considered in the design of actual elastomer damper 
hardware. These factors include both the general design considerations 
(simplicity, ease of assembly, durability and cost) and considerations peculiar 
to elastomer dampers, such as the choice of geometric configuration for the 
elastomer elements, and whether those elements should be permanent or replace- 
able. Elastomer dampers for control of unidirectional vibrations usually 
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involve simple foot mounting of machines and thus require rather straightforward 
designs* Occasionally, constrained layer damping (e.g, surface coatings) 
may be appropriate [j7.2\ 

The required damper designs are generally more complex for rotordynamic sys- 
tems; space limitations and maintenance requirements tend to be stricter for 
rotordynamic dampers. For most applications, a segmented elastomer cartridge 
configuration (with either shear or compression buttons) is recommended. The 
segmented cartridge is well suited for the design of compact dampers with 
replaceable elastomer elements. While the use of replaceable elastomer 
elements results in a small increase in the number of components in the damper, 
this factor is more than offset by the resulting simplification of the assembly 
and alignment procedures and the improved maintainability of the damper. The 
testing of alternative elastomer element materials and dimensions is also 
possible with replaceable elements, permitting in-service optimization of the 
damper. The elastomer elements may also be easily replaced after a specified 
length of running or calendar time during scheduled maintenance operations, 
thus minimizing the effects of elastomer aging, work hardening, and fatigue. 

0-rings may also be utilized as replaceable elastomer elements. However, 

0-rings are not as convenient to use since they generally require splicing 
when being replaced. Also, splicing of 0-rings results in some variation in 
circumference which can, in turn, produce some variation in the dynamic pro- 
perties of the damper. 

Several examples of actual elastomer damper designs with replaceable com- 
pression or shear button elastomer elements are presented in Section 8.0. 

Test results and operating experience with these dampers are also discussed. 

7 . 5 Summary 

Several steps are involved in the proper design of elastomer dampers for 
rotating machinery. A recommended set of steps is summarized in the flowchart 
in Figure 120, corresponding to the steps listed below: 

1. Review elastomer materials available for the required 
parameters and operating environment (temperatures, 
lubricant compatibility, etc.). 

2. Construct a rotordynamics model of the shaft system 
in its bearings. 

3. Compute the damped natural frequencies and logarithmic 
decrements of the shaft system for a range of support 
coefficients within the elastomer material limits. 

4. Choose optimum support properties in terms of elastomer 
stiffness and loss coefficient. 

5. Confirm this choice of parameters by performing an 
unbalance response analysis to obtain bearing loads, 
steady-state orbits, etc. 
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Fig. 120 Flow Chart Showing Steps of Elastomer Damper Design 
(Numbers Refer to Number Steps In Text) 
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6 e Select candidate materials and configurations for the 
damper; perform tests to confirm properties and suit- 
ability for the expected operating environment (if such 
data are not available) - 

7. Lay out possible configurations for the most promising 
elastomer; select dimensions to give stiffness and 
damping properties as close as possible to optimum values 
for each material; and evaluate the performance of the 
chosen configurations (using, for example, unbalance 
response analysis) . 

8. If the results of step 7 are satisfactory, choose one 
or more configurations and design the damper in detail; 
if the results of step 7 are not satisfactory, repeat 
steps 6 and/or 7 until acceptable combinations are found. 

9. Build and test the dampers to experimentally verify their 
predicted performance, and make final selections. 

Several specific examples of this design procedure and of the resulting 
damper hardware and test results are presented in Section 8.0. 
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8.0 EXAMPLES OF ELASTOMER DAMPER DESIGNS 


The step-by-step procedure for the design of elastomer dampers for rotating 
machinery was discussed in detail in Section 7.0. Three such dampers, which 
were designed using this procedure, are described in this section. These 
dampers were built, tested, and operated satisfactorily. The three dampers 
described were used on a gas turbine simulator [8.ll*, a power transmission 
shaft [8.2J, and an actual power turbine shaft fronT a gas turbine engine 
[8.3], 

8.1 Gas Turbine Simulator f8.lj 

An elastomer damper was applied to a gas turbine simulation rig adapted 
from an earlier test rig as described in Reference 8.4. The rotor was 
designed to represent the power turbine of an advanced gas turbine engine 
and consisted of a 2.22 cm (7/8 in.) x 58.4 cm (23 in.) shaft with a 26 cm 
(10.25 in.) diameter disc mounted at one end. A small air turbine at the 
opposite end was used to drive the rotor to a maximum speed of 30,000 rpm. 
Figure 121 shows this rotor during assembly. The rotor was supported by a 
pair of angular contact ball bearings at either end, originally mounted in 
pedestals with metallic flexures and damped with squeeze-film dampers. 

These pedestals were put aside, and a new set was designed to substitute 
elastomeric mounts for the original flexure-squeeze-film arrangement. 

The new mounts were designed by creating a mathematical model of the rotor. 
This model, shown schematically in Figure 122, consisted of 31 stations 
defining 30 discrete sections. Each section was defined by length, outer 
diameters, inner diameter, elastic modulus, density and shear modulus. 
Additionally, discs described by concentrated mass and inertia values were 
located at appropriate stations. The mathematical model of this rotor is 
shown in Table 24. Elastomeric pedestals were represented by radial stiff- 
ness and equivalent viscous damping coefficients. The ball bearings were 
represented by a radial elastic stiffness and a small viscous damping 
coefficient. The bearings and pedestals were assumed to have had zero 
angular stiffness and damping. 

The mathematical analysis of the elastomer damper test rig was composed of 
three phases. The first phase consisted of a damped critical speed analysis 
8.5j which allowed the optimum support properties to be determined, along 
with the selection of elastomers to be tested. The second phase consisted 
of performing a response analysis to compute expected unbalance response of 
the test rig ! 8.6 j, and the third phase involved the actual design of the 
elastomer supports. 

A damped natural frequency analysis was performed to assess the values 
required for elastomer stiffness and damping. The pedestals were assumed to 
have stiffnesses of 1.75 x 10^, 5.25 x 10^ and 1.75 x 10^ N/m (10^, 3 x 10^ 
and 105 lb/in.) at both ends. Damped natural frequencies were computed for 
pedestal damping coefficients of 1750, 3500 and 8750 Ns/m (10, 20, and 50 
lb-sec/in.). Table 25 shows the results of these calculations. Damped 


* Numbers in brackets indicate references found in Section 8.5. 
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TEST ROTOR MATHEMATICAL MODEL 
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natural frequencies were calculated by assuming a rotor speed and calculating 
the natural frequencies corresponding to that speed. Damped natural frequen- 
cies are complex quantities and of the form o + ito, where a denotes system 
stability (negative o implies stability) and w denotes the frequency as the 
solution form, is exp (a + iw) T . 

The stiffness range of 1.75 x 10 6 to 1.75 x 10 7 N/m (10^ to 10 5 lb/in.) was 
established as the practical elastomer stiffness range for this size and type 
of rig. 

For a particular elastomer material, the loss coefficient (ratio of damping, 
coB, to stiffness, K) is relatively fixed. The equivalent viscous damping, 

B, generally falls off rapidly with increasing frequency. Consequently, the 
decision was made to target the elastomer damper design for the higher stiff- 
ness of 1.75 x 10 7 N/m, which maximizes the available damping. 

If the log decrement is plotted against pedestal damping for an elastomer 
stiffness of 1.75 x 10 7 N/m (100,000 lb/in.), the curves in Figure 123 are 
the result. From this figure, it is apparent that an optimum value of 
damping exists, approximately 1.75 x 10^ N-sec/m (100 lb-sec/in.), and that 
the corresponding log decrement is 0.965. Figure 124 shows how damping, B, 
varies with the loss coefficient, n. For a loss coefficient of 0.8, which 
is typical for fluorocarbon (Viton 70 - CFM) at 70°F, the equivalent viscous 
damping at 20,000 rpm is 6655 N-sec/m (38 lb-sec/in.), and, for a loss 
coefficient of 0.15, which is typical for butadiene (Polybutadiene - BR) at 
70°F, the damping is 1261 N-sec/m (7.2 lb-sec/in.)*. 

For this series of tests, two elastomer materials, fluorocarbon (Viton - CFM) 
and butadiene (Polybutadiene - BR) , were considered since they represent the 
practical extremes of elastomer capability in terms of internal damping. Both 
elastomers incur a penalty with respect to the damping, as compared to the 
optimum log decrement. Table 26 presents this penalty by comparing the 
expected log decrement of an elastomer damper with the optimum value of log 
decrement. The table shows a larger penalty with butadiene (Polybutadiene - 
BR) than with fluorocarbon (Viton -CFM). 


A series of unbalance response calculations were performed to estimate the 
unbalance sensitivity of the elastomer-mounted rotor. Unbalance may be 
added in the form of threaded masses screwed into tapped holes at several 
locations along the rotor (Figure 120) . Mass addition locations are at the 
large disc, two discrete locations along the shaft and at the drive turbine. 
Figures 125 through 128 illustrate the results of unbalance response calcula- 
tions using an unbalance of 2.54 gram centimeters (1 gram- inch) in either 
the disc or the shaft (Stations 29 or 11, respectively, in Figure 122). These 
figures show the response of the disc (Station 31 of Figure 125) for a mass 
added to the disc or the response of the center of the shaft (Station 13) for 
a mass added to the shaft. 


* The use of equivalent viscous damping to represent hysteretic damping for 
a single frequency is discussed in Section 2.0. 
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ELASTOMER STIFFNESS = 1.75 * 10 7 N/m (100,000 Ib/in.) 
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Fig. 123 Elastomeric Damper Test Rig, Log Decrement as a Function of Elastomer Damping 
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Fig. 124 Damping Coefficient as a Function of Loss Coefficient 


Table 26 


COMPARISON OF MODE-2 LOG DECREMENT FOR BUTADIENE (POLYBUTADIENE - 
AND FLUOROCARBON (VITON 70 - CFM) WITH OPTIMUM DAMPING 


(1) Optimum 


j N-sec 
1 m 


lb-sec 1 
i in. / 


6 


17513/100 .965 


(2) butadiene 

(Polybutadiene - BR) 1261/7.2 


.075 


(3) fluorocarbon 

(Viton 70 - CFM) 


6655/38 


.31 
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Figures 125 and 126 show the calculated response plots of amplitude against 
rotor speed for the test rig on butadiene (Polybutadiene - BR) mounts (loss 
coefficient n - ,15). Figures 127 and 128 show the response to these same 
unbalances, but with fluorocarbon (Viton 70 - CFM) supports (r\ - .8). These 
sets of curves indicate that the second critical should be the most trouble- 
some and that the rig should be much more sensitive to unbalance when mounted 
on the butadiene (Polybutadiene - BR) supports. Table 27 gives the computed 
values of the rig critical speeds and the expected sensitivity to unbalance 
in units of microns (peak-to-peak) per gram of unbalance at the balance hole 
radius . 


The elastomer supports were designed by mounting each bearing housing on 
three elastomer cartridges spaced 120° apart, as shown in Figure 129a. The 
elastomers 1 elements were composed of one or more buttons, glued to upper and 
lower platens to form easily replaceable cartridges as shown in Figure 129b. 
The number, diameter, and thickness of the buttons in each cartridge were 
chosen to satisfy the requirement of a 1.75 x 10? N/m (100,000 lb/in.) 
overall stiffness at 20,000 rpm under specified ambient conditions. The 
elastomers were preloaded about 10% by a pair of preload screws on each 
cartridge so the elastomer buttons would always remain in compression. 

For the elastomer mounting configuration shown in Figure 129, the overall 
stiffness was calculated by using the relationship defined in equation (7-18) . 


In considering the design for individual buttons with the required support 
stiffness of 1.75 x 10^ N/m (100,000 lb/in.), equation (7-18) gives 


N (K + K ) - 11.7 x 10 
B e s 


N/m (66,666 lb/in.) 


( 8 - 1 ) 


The empirical correlations for the compressive and shear stiffness of 
butadiene (Polybutadiene - BR) buttons presented in Section 6.0 were not 
available at the time that this damper was designed. Consequently, less 
complete correlations for the shear storage modulus and shape factor 
developed in Reference 8.7 were used and are given below: 


K 

c 


ttD 


fl + 12 . 33u 0-29 



( 8 - 2 ) 


K = G' 
s 



butadiene (Polybutadiene - BR) at 32°C; 
G' = 3.686 x 10 6 a) 0 * 2037 Pa (N/m 2 ) 


(8-3) 


(8-4) 


butadiene (Polybutadiene - BR) at 50°C; 
G’ = 1.902 x 10 6 m°- 2727 Pa(N/m 2 ) 


(8-5) 
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Table 27 
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Fig. 129 Elastomer Damper Construction 




where D is the diameter and h the height of a cylindrical elastomer button 
(meters) ; where G ? is the shear storage modulus (N/m^) and m is the frequency 
(rad/sec) . While the shear modulus is frequency-, temperature-, and dynamic- 
strain-dependent, the temperature and strain levels encountered in this test 
rig were similar to those used in Reference 8.7. Thus, the correlations from 
Reference 8.7, which were not complete in terms of temperature and strain, 
were shown to be adequate for this application. 

Equations (7-18), (8-2), and (8-3) were combined to find the radial stiffness 
for butadiene (Polybutadiene - BR) buttons as 

2 

= 4 . 5N g G ' ll.33 + 12.33~°' 29 (|^) 2 (8-6) 

Similar test data was not available for fluorocarbon (Viton 70 - CFM) ; 
however, the radial stiffness correlation for butadiene (Polybutadiene - BR) 
given in Equation (8-6) was assumed to apply to fluorocarbon (Viton 70 - CFM) . 
Only a value for the shear modulus at the frequency of interest had to be 
found. The static shear modulus for 70-durometer elastomers (from Reference 
8.8) is 1.861 MPa (270 lb/in. 2); (see also Figure 64). As the squeeze film 
is supported on fluorocarbon (Viton - CFM) 0-rings, a dynamic multiplier of 
6.084 is necessary to correct the static shear modulus to the operating 
condition of 2,000 rad/sec at 32°C, Reference 8.9. An additional factor of 
0.614 is also required to correct for an operating temperature of 50°C at that 
frequency. Table 28 lists the assumed values of shear modulus for both 
elastomer types at 2,000 rad/sec. 

Equation (8-6) was applied for the appropriate elastomer buttons design. 
Figures 130 through 133 plot radial stiffness against button diameter for 
the case of a single button per cartridge (Ng - 1) . Each figure consists of 
three curves corresponding to elastomer thicknesses of 2.38, 3.18 and 4.76 mm 
(3/32, 1/8 and 3/16 inch). Figure 130 corresponds to fluorocarbon (Viton 70 - 
CFM) at 32°C, Figure 131 to fluorocarbon (Viton 70 - CFM) at 50°C; Figure 132 
to butadiene (Polybutadiene - BR) at 32°C; and Figure 133 to butadiene 
(Polybutadiene - BR) at 50°C. The required button diameter for any given 
stiffness is readily found from these curves. Table 29 gives the necessary 
stiffness per button to achieve the 100,000 lb/in. that was desired for the 
overall stiffness value. 

Table 30 gives the necessary button diameters, based on Table 29 and Figures 
130 though 133, to achieve the 1.75 x lO' 7 N/m (100,000 lb/in.) stiffness. 

Preliminary layout work indicated that a 3-button-per-cartridge arrangement 
using 1/8- inch thick elastomer stock would be dimensionally convenient. This 
combination was selected, and, for simplicity, the decision was made to use 
a single diameter for fluorocarbon (Viton - CFM) and butadiene (Polybutadiene- 
BR) buttons. A diameter of 15 mm (0.59 in.) was selected (see Table 30) to 
insure that a minimum stiffness of 1.75 x 10? N/m (100,000 Ib/in.) existed 
at all times. 

Figure 134 shows the elastomer cartridges being assembled. On the right are 
two platens and three individual buttons, while the assembled cartridge can 
be seen on the left. Figure 135 shows the two pedestals and bearings: the 
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Table 28 


VALUES OF SHEAR MODULUS G' AT 2000 RADIANS /SECOND 


Elastomer 


Temperature 

(°C) 


Shear Modulus 



fluorocarbon (Viton 70 - CFM) 
fluorocarbon (Viton 70 - CFM) 
butadiene (Polybutadiene - BR) 
butadiene (Polybutadiene - BR) 


32 

11.33/1,643 

50 

6.96/1,009 

32 

17.31/2,511 

50 

14.00/2,031 


Table 29 

STIFFNESS VALUES TO BE USED IN FINDING BUTTON DIAMETERS 


Number of Elastomer 
Buttons Per Cartridge 

1 

2 

3 

4 


Stiffness Value to be Used in 

Figures 140 through 143 (— /— • ) 
m in. ■ 

1.75 x 10 7 /100,000 

8.76 x 10 6 /50,000 
5.84 x 10 6 /33,333 
4.38 x 10 6 /25,000 
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STIFFNESS (Ib/in.) 





STIFFNESS (Ib/in.) 



Fig. 132 Plot of Support Stiffness versus Button Diameter for Butadiene (Polybutadiene) at 32° C 
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STIFFNESS (N/m) 



STIFFNESS (Ib/in.) 









Table 30 


BUTTON DIAMETERS TO ACHIEVE A 100,000 LB/IN. RADIAL STIFFNESS 


Button 

Number of 

flourocarbon 

Thickness 

Buttons Per 

(Viton) 

(mm/ in • ) 

Cartridge 

at 32°C 


1 

14. 7/. 58 

2. 38/^| 

2 

12. 2/. 48 

3 

10. 7/. 42 


4 

9. 9/. 39 


1 

18.0/. 71 

3.i8/i 

2 

3 

14. 7/. 58 
13.0/. 51 


4 

11. 7/. 46 


1 

24. 4/. 96 

4 ' 76/ lf 

2 

3 

10. 6/. 77 
16. 8/. 66 


4 

15. 2/. 60 


flourocarbon butadiene butadiene 

(Viton) (Polybutadiene) (Polybutadiene) 
at 50°C at 32°C at 50°C 


17.0/. 67 

13. 2/. 52 

14.0/. 55 

14.0/. 55 

10. 7/. 42 

11. 4/. 45 

12. 4/. 49 

9. 4/. 37 

9. 9/. 39 

11. 4/. 45 

8. 6/. 34 

9.1/. 36 

21. 6/. 85 

16.0/. 63 

17.0/. 67 

17.0/. 67 

13.0/. 51 

13. 7/. 54 

15.0/. 59 

11. 2/. 44 

12. 2/. 48 

13. 7/. 54 

10. 2/. 40 

10. 9/. 43 

28.2/1.11 

21. 3/. 84 

22. 6/. 89 

22. 9/. 90 

16. 5/. 65 

18.0/. 71 

20.1/. 79 

14. 5/. 57 

15. 5/. 61 

18.0/. 71 

13.0/. 51 

14. 2/. 56 
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disc end on the left and the turbine end on the right. Additionally, the 
figure shows two preload screws with their calibrated washers (the washers 
were trimmed on assembly to get the correct preload) as well as a partially 
assembled elastomer cartridge. Figure 136 shows the disc end bearing 
housing assembled within its pedestal. Note that solid steel blocks have 
been inserted instead of the elastomer cartridges. These blocks allowed 
comparisons in unbalance response to be made between hard-mounted and 
elastomer-mounted rotors. 

Figure 137 shows a side view of the test rig as finally mounted and instru- 
mented in the test cell. The disc is inside a vacuum box on the left, and 
the drive turbine is on the right. Figure 138 shows the rig from the 
front, with the vacuum hoses in the foreground. Figure 139 shows 
a detailed view of the turbine-end bearing housing and pedestal. Two elas- 
tomer cartridges are visible in the 12 o’clock and 4 o’clock positions. The 
two noncontacting probes were used to measure the motion of the housing 
relative to the pedestal, and the two visible flexible tubes are the bearing 
oil feed and drain lines. 

A series of tests were conducted with this test rig to evaluate the effec- 
tiveness of the elastomer dampers. The test procedure was as follows: 

1. Install butadiene (Polybutadiene - BP) damper cartridges 
at both locations (turbine end and disk end) . 

2. Balance the rotor. 

3. Create an unbalanced condition in the disk in four discrete 
steps. 

4. Create an unbalanced condition in the shaft in four separate 
steps after the unbalance has been removed from the disc. 

5. Remove shaft unbalance and install fluorocarbon (Viton 70 - 

CFM) damper cartridges in place of the butadiene (Polybutadiene - 
BR) cartridges . 

6. Rerun with disc unbalances, 

7. Repeat shaft unbalance runs. 

8. Replace fluorocarbon (Viton - CFM) dampers with solid steel 
blocks and run with disc unbalance, 

9. Replace turbine end steel block with fluorocarbon (Viton - CFM) 
damper and run with a limited number of unbalances. 

10. Record data from selected runs on magnetic tape. 

11. Plot displacement amplitude against rotor speed for test 
runs. 

12. Photograph disc and shaft orbits for selected runs. 
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Fig* 137 Side View of Assembled Test Rig 









The plan for reducing elastomer damper test rig data involved calculation 
of sensitivity to unbalance masses, calculation of modal damping values, 
and comparing these values with predictions. The sensitivity, S u , to un- 
balance masses is defined as: 


S 

u 


Modal Response 
Unbalance Magnitude 


Mils P-P 
gram 


(8-7) 


The sensitivity for a particular mode depends on the locations of both the 
unbalance mass and the measurement probe. This sensitivity is valuable in 
that it is an indication of the difficulty to be expected in balancing a 
particular mode (at that location). 

Damping information was extracted from the test runs by making use of the 
width of the resonance curve at the half-power points and also by using 
the rate of change of the phase angle through the resonant speed. With 
the half-power point method, the resonant speed, f n (for a particular 
critical speed), is located by finding the peak of the resonance curve for 
that mode and noting the frequency at which it occurs. The half-power 
points are the two locations on the resonance curve where the amplitude is 
equal to the peak amplitude divided by the square root of two. Then Af is 
defined as the difference in frequency between the two half-power points. 
With this information, the system log decrement, 6, may be found with the 
equation : 

6 — tt (8-8) 

n 

An alternate method takes advantage of the fact that the logarithmic decre- 
ment is inversely proportional to the rate of change of the phase angle 
through the resonant speed. Equation (8-9) expresses the relationship 
between the log decrement, the resonant frequency, and the rate of change 
of phase through the resonant speed. 


6 


360 



(8-9) 


Using either equation (8-8) or (8-9), the quality or amplification factor, 
Q, may be calculated from 


Q * f (8-10) 

The initial or baseline balancing was done with butadiene (Polybutadiene - 
BR) dampers using an influence coefficient balancing program. Because the 
first critical speed could not be initially negotiated, influence coeffi- 
cients were calculated at 7,000 rpm for the first three correction mass runs. 
Correction masses were only placed in the disc plane. 
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Fig. 139 Close-up of Turbine in Bearing Housing Showing Elastomer Cartridges at 12 o ? clock and 
4 o f clock Positions 



Influence coefficients were then calculated at 11,500 rpm and one addi- 
tional correction mass was installed. The next correction mass run cal- 
culated a mass of 0.058 grams which was too small to install* The vibra- 
tion level was low throughout the speed range (less than 20 microns (.001 
inch or 1 mil) peak to peak) . Test runs were made with unbalance masses in 
either the disc (plane 1 of Figure 140) or in the shaft (plane 3 of 
Figure 140) . 

Table 31 shows the magnitudes of the unbalance masses used. In the disc, 
four separate and distinct masses were used alternately in the same holes, 
but the small hole size on the shaft required simultaneous use of from one 
to four masses in adjacent holes to obtain the necessary unbalance. The 
corresponding shaft unbalance mass shown in Table 31 is the equivalent 
mass (if only one hole were used? i.e., the vector sum of the unbalance 
masses) . 


Table 31 

SUMMARY OF MASSES USED TO UNBALANCE 
THE TEST RIG DISC AND SHAFT 


Unbalance in Disc 

Unbalance in Shaft 

(grams) 

(grams) 

0.2 

0.183 

0.3 

0.353 

0.4 

0.583 

0.5 

0.621 


Figure 141 shows the installation of one of the shaft unbalance masses. 

Table 32 shows the location of the first two critical speeds of the test rig 
with the two different damper materials. The critical speeds were higher 
with the fluorocarbon (Viton - CFM) dampers than with the butadiene (Poly- 
butadiene - BR) dampers. This indicates the fluorocarbon (Viton - CFM) 
dampers were stiffer than the butadiene (Polybutadiene - BR) . The butadiene 
(Polybutadiene - BR) shear moduli used in the design (Equation 8-4) were 
based on tests of a previous material batch although while changes in carbon 
black had been made between batches by the manufacturer. Moduli for the 
pertinent batch, given in Section 6.0, were only 60% of the values in 
Equation (8-4) . 


Table 32 

LOCATIONS OF CRITICAL SPEEDS FOR DIFFERENT ELASTOMER SUPPORTS 


Elastomer 


First Critical 
Speed 


10,200 

12,050 


Second Critical 
Speed 


butadiene (Polybutadiene - BR) 
fluorocarbon (Viton 70 - CFM) 
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20,500 

24,000 




2 






Figures 141 through 145 show the rig's response to the unbalances intro- 
duced in the different planes with the two different damper materials. 

Note that the locations of the critical speeds drop in frequency with 
increasing unbalance, indicating a small strain-softening effect in the 
elastomer mounts as described in Section 4.0. At highest strain, the drop 
in resonant frequency is approximately 8%, indicating at least 16% loss 
in stiffness (and probably more). Peak amplitudes at the damper were 
approximately 10 to 20 microns (0.4 to 0.8 mils), which is a dynamic strain 
of 0.003 to 0.007 in a 3.2-mm button. Previous test results [s.l indicate 
that such a strain could result in a 9 to 32 percent loss in stiffness in a 
compression specimen. Thus, the test results for vibration control testing 
of an elastomer damper were consistent with the earlier component test 
results for strain effects. 

The influence of unbalance mass on peak amplitude is shown in Figure 146. 

The slopes, which are equal to the unbalance sensitivities, are approxi- 
mately the same for fluorocarbon (Viton - CFM) and for butadiene (Polybuta- 
diene - BR) in this application; this is a discrepancy from the predicted 
behavior with these materials discussed previously. There is some reduction 
in sensitivity with increasing amplitude; this reduction is consistent with 
the increase in log decrement with increasing strain observed in earlier 
tests. Table 33 compares the observed and predicted sensitivity to unbalance. 

In all cases the rotor is less sensitive than predicted. This lesser sensitivity 
indicates that the system damping achieved by the elastomer dampers exceeded 
expectations of the original system dynamic model. Two reasons are hypo- 
thesized. First, the computer model assumed purely radial stiffness and 
damping, with no moment restraint. In fact, the elastomer cartridges and 
double-row ball bearings probably provided some finite amount of moment 
stiffness and damping in addition to the radial values. Second, the elasto- 
mers ran at average temperatures of 27°C and 16°C at the disc and turbine 
end, respectively (as opposed to the design. temperature of 50°C) , contri- 
buting to both greater stiffness and damping than was originally expected. 

Table 34 shows the calculated values of the log decrement for butadiene 
(Polybutadiene - BR) and fluorocarbon (Viton - CFM) , using both the half- 
power point and the rate of change of phase angle methods, and the dis- 
crepancies between the two methods used. The agreement was generally good. 

The average log decrements and Q values are shown in Table 35, The predicted 

values for the logarithmic decrement 6 in mode two are: 5 - - 0.057 and 

6 . =0.298. P0ly 

Viton 

The repeatability of the test rig during acceleration and deceleration is 
shown in Figure 147. This figure was obtained with the damper material 
(butadiene (Polybutadiene - BR)) at approximately equilibrium temperature 
with a 0.4-gram unbalance in the disc. This plot is a good representation 
of the rig T s repeatability. The only significant deviations observed 
occurred when there were large temperature differences in the elastomers 
between the initial run up-to-speed and the coast-down, which were avoided 
during testing by running the rotor to speed at least once before a test 
sequence began. 
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Fig, 142 Effect of Disc Unbalance With Butadiene (Polybutadiene) Dampers 
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Fig. 143 Effect of Shaft Unbalance With Butadiene (Polybutadiene) Dampers 
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Fig. 144 Effect of Disc Unbalance With Fluorocarbon (Viton-70) Dampers 
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Fig. 145 Effect of Shaft Unbalance With Ruorocarbon (Viton-70) Dampers 
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Fig. 146 Sensitivity of Test Rig to Unbalance 
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Fig. 147 Repeatability of Elastomer-Mounted Rotor Running To and From Speed 
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* Signal level too low 

** Approximate values, could not reach second half-power points 
***Diff iculty with taped data 



AVERAGE VALUES FOR LOG DECREMENT AND AMPLIFICATION FACTOR 
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After the elastomer damper data were taken, solid steel blocks were sub- 
stituted in place of the elastomers. See Figure 136 for the pedestal with 
steel blocks in place of the elastomer cartridges. With the steel blocks 
installed, the rig had only one resonance in the speed range, at approxi- 
mately 19,000 rpm, and greatly increased response amplitude. 

Figures 148 and 149 show the response at the disc and shaft, respectively, 
when an unbalance mass was added to the disc. The data were not taped or 
digitally acquired so the phase change is unknown, but the response appeared 
to be nonlinear. The log decrement values are included in Tables 34 and 35, 
but confidence in these numbers is limited because of the support non- 
linearity. It was noticed during these runs that the shaft was not as free 
to rotate with the steel blocks installed as it was with the elastomer 
dampers. The increased resistance was the result of a small amount of 
binding in the bearing caused by misalignment of the pedestals and shaft. 

The relatively soft elastomer mounts were able to tolerate this misalignment, 
but the stiff steel blocks could not. 

Figure 150 shows the rig response with fluorocarbon (Viton - CFM) dampers in 
the turbine end and steel blocks at the disc end. This arrangement simulated 
the use of elastomeric supports at only the cold end of a gas turbine power 
shaft. The response amplitude is between the values observed with both 
elastomers' steel blocks at both ends. Thus, damping benefits can be gained 
with elastomers at the cold end only. 

The amplitude data for the response plots were handled in two distinct ways. 

In some cases the probe output was not taped, and the plots were made "on-line”, 
directly from the probe output after conditioning by a tracking filter to 
extract the synchronous component. In other cases, the signals were tape- 
recorded, and the tape-recorded signals were analyzed by feeding them into a 
spectrum analyzer and using the analyzer to create a "peakhold" average spec- 
trum for each run, which was plotted using a digital plotter. The tape- 
recording method is valid when the synchronous component of the signal is 
larger than any other signal components at all speeds so the analyzer is 
effectively acting as a tracking filter. This was the case for the elastomer 
rig, and two sample instantaneous spectra are shown in Figure 151. The upper 
spectrum shows the frequency content of the disc probe and the lower spectrum, 
the content of a probe at the turbine. For both cases, the rotor speed is 
6,000 rpm, well away from the first critical speed where the synchronous 
component is greatly amplified. The disc probe shows a large signal at run- 
ning speed and very little elsewhere. The turbine probe shows very low 
vibration and, indeed, was not considered in the data analysis because of 
its low amplitude for both modes. However, the output of this probe is 
interesting because it shows very small peaks at multiples of running speed. 
These peaks are the result of forced vibration due to the impact of the 
compressed air jets on the machined "buckets" of this Impulse air turbine. 
Traces of these spikes can be seen in the disc spectrum as well but, at all 
frequencies, were exceeded by the amplitude of the synchronous component. 

Figure 152 shows the orbit of the disc at three different speeds: below, at, 

and above the first critical speed of the rig when a 0.5-gram unbalance mass 
was installed in the disc (Plane 1 of Figure 140) . Figure 153 shows the 
orbit of the shaft center at three speeds which are below, at, and above the 
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Fig. 150 Effect of Using Viton Damper Cartridges at the Turbine End Only 
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Fig. 152 Disc Orbit Through First Critical With Unbalance in Disc 
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Fig. 153 Shaft Orbit Through Second Critical With Unbalance in Shaft 




second critical speed. The unbalance for this case was 0.621 grams in the 
shaft (Plane 3 of Figure 140). These photographs are presented to illus- 
trate the fact that the orbits are essentially circular; therefore, the 
elastomer support can be considered isotropic. The maximum amplitudes of 
the orbits for both figures are given in Table 36. 


Table 3 6 

SPEEDS AND MAXIMUM AMPLITUDE OF SHAFT ORBITS 


Probe 

Location 

Shaft Speed 
(rpm) 

Maximum Amplitude 
(microns/mils p-p) 

Disc 

9,000 

56/2.2 

Disc 

10,950 

94/3.7 

Disc 

12,000 

33/1.3 

Shaft 

21,000 

51/2.0 

Shaft 

23,000 

79/3.1 

Shaft 

26,500 

46/1.8 


The small irregularities in the orbits are due to a phenomenon called "elec- 
trical runout", in which small irregularities in surface finish or variations 
in electrical conductivity of the rotor are seen by the eddy-cur rent- type 
displacement pickup as if they were vibration signals. Since they are fixed 
to the rotor surface, these small "ripples" in the orbit make phase changes 
easily visible. This is seen in Figure 153a where a small bump appears at 
the top of the orbit. In Figure 153b, the bump appears at some midway point 
(apparently smaller since the orbit has enlarged but the "electrical runout" 
signal has not). This bump appears at the bottom in Figure 153a. These 
correspond to speeds below, at, and above the second critical speed, 
respectively. Note also that a total phase change of about 180° has occurred 
while traversing the resonance. 

8.2 Supercritical Power Transmission Shaft [8.2 j 

A test rig was designed and built for demonstrating the use of supercritical 
power transmission shafting. During the operation of this test rig, nonsyn- 
chronous vibrations from a variety of sources were observed. The test rig 
was first run with a minimum amount of external damping and then with a 
larger amount of external damping. In this way, the effect of external damp- 
ing on nonsynchronous vibration and rotordynamic stability was observed. 

The external damping was provided by a self-contained, oil squeeze-film 
damper and an elastomer damper, both designed to permit operation through the 
fourth bending mode. The elastomer damper is of principal interest here, but 
the oil damper was used as the basis of comparison of damping effectiveness. 
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Tests were performed to identify any problems, needed technology, or limita- 
tions inherent in the use of supercritical shafting. One of the principal 
areas of interest involved demonstrating the balancing of supercritical 
shafts. The tests included balancing of a very flexible shaft with both 
damped and undamped supports. 

A sketch of the test rig is presented in Figure 154. A 224-kW (300 hp) 
electric motor is the prime mover for this test rig. The motor drives a 
variable-speed magnetic coupling, whose output speed is continuously variable 
from 50 to 3,600 rpm. A gearbox with a ratio of approximately 5.7 to 1 
produces a drive-shaft speed of up to 20,000 rpm. The drive shaft is formed 
from a section of aluminum tubing, 7.62 cm (3 in.) in diameter and 3.66 m 
(12 ft.) in length with a wall thickness of 3.175 mm (0.125 in.). 

Early testing without a damper did not permit rig operation above 1,200 rpm 
due to the instability of the subsynchronous excitation of the 16-Hz first 
critical speed. Stability problems which occurred during this stage of 
testing clearly established that some form of external damping, coupled with 
balancing, would be required for control of vibrations. Several methods for 
applying external damping were investigated, and an oil-squeeze-f ilm-type 
damper was chosen first based upon its flexibility and predictability in the 
level of damping achieved [8.10 through 8.13j. This damper was subsequently 
replaced with an elastomer damper which proved to be equally effective. 

For this particular application, some special considerations arose beyond 
the normal considerations discussed in earlier sections. These considera- 
tions were that the damper should have the following capabilities: 

• Easily added to the existing test rig. 

• Compact and self-contained to be practical in helicopter 
drive shaft applications. 

• Providing support stiffness as well as damping. 

The approach selected was to flexibly couple an extension shaft with a damper 
to the existing test shaft and perform undamped critical speed analysis of 
the new shaft geometry. Figure 155 shows the predicted mode shapes and 
critical speeds for the first five modes. It should be noted that, for all 
of these modes, some significant amplitude exists at the location of the 
damper as represented by bearing number 2. This is important to achieve the 
necessary dissipation of vibration energy at the damper. 

The initial model of the test rig included a support of 1.75 x 10^N/m (1,000 
lb/in.). A damped natural frequency analysis was conducted using a series of 
damping values from 875 to 10,500 N-s/m (5 to 60 Ib-sec/in.) with this stiff- 
ness. The second critical speed was found to be critically damped. The log 
decrement is plotted as a function of damping for the remainder of the first 
five critical speeds in Fiqure 156. The optimum value of damping appears to 
be about 8,750 N-s/m (50 lb-sec/in.). 

The damped natural frequency analysis was then conducted for a series of 
values of support stiffness from 1.75 x 10^ to 2.625 to 10^ N/m (1,000 to 
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Fig. 155 Mode Shapes of the First Five Critical Speeds, Original 
Shaft with Extension to Carry Squeeze-Film Dampers 
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Fig. 156 Log Decrement as a Function of Damping 
(Stiffness = 1.75 x 10 s N/m [1000 lb/in.]) 


15,000 lb/in.) with a constant damping value of 8,750 N-s/m (50 lb-sec/in*)- 
Again, the second critical speed was critically damped* Log decrement is 
plotted as a function of stiffness for the remainder of the first five 
critical speeds in Figure 157. The only critical speed for which there was 
a substantial detrimental effect from increasing support stiffness was the 
first critical speed. However, this effect was serious enough to require 
that the support stiffness be kept below 8.175 x 10 5 N/m (5,000 lb/in.). 

Based on the rotordynamic analysis, the decision was made to first design a 
squeeze film for a stiffness of 7 x 10^ N/m (4,000 lb/in.) using an 0-ring 
as the spring member. Studies recently performed by Smalley, Darlow and 
Mehta (j8.14j indicate that this is about as low as can be practically 
achieved for a continuously supported 0-ring in this configuration. 

The squeeze-film damper was designed to be a sealed damper with no circula- 
tion of oil. In this way, the damper would require no support hardware (such 
as oil supply pumps) which would prohibit its use in helicopter or other 
space- and weight-limited applications. 

A sketch of the test rig damper is shown in Figure 158. The squeeze-film 
damper was evacuated, filled and sealed, and pressurized by the use of a 
bladder. The damper was predicted to generate less heat than the bearings, 
so the dissipation of heat was not expected to be a serious problem. The 
0-ring retainers are radially adjustable, so that the squeeze film damper 
may be centered manually to compensate for static deflection of the 0-rings. 
The damper was designed to achieve a level of damping in the range of 8,750 
N-s/m (50 lb-sec/in.) for silicone oil with a viscosity of about 80 centi- 
stokes, using the short bearing theory and assuming no cavitation. The 
radial clearance is 0.635 mm (25 mils) and the length and diameter are 7.62 cm 
(3 in.) and 10.16 cm (4 in.), respectively. Although the squeeze film damper 
was designed to allow for the two 0-rings, it was used initially with just 
one 0-ring to obtain a low parallel support stiffness of 7 x 1GP N/m (4,000 
lb/in.) for a continuously supported 0-ring. 

Testing showed that the squeeze- film damper improved the dynamic stability 
of the test rig. However, the addition of external damping alone could not 
provide complete control of all test shaft vibrations. The test shaft also 
had to be balanced to achieve safe operation throughout its speed range. 

Through the combined use of external squeeze- film damping and balancing, the 
rig was safely run to over 12,000 rpm while negotiating the first four flex- 
ural critical speeds of the test shaft. This speed was over 10 times that 
which could be achieved without any external damping. 

Substantial synchronous growth of the response of the fifth critical speed 
occurred at 12,000 rpm. A very large subsynchronous response of the first 
and fourth critical speeds also occurred. In fact, the subsynchronous 
response of the first critical speed, while running at 12,000 rpm, was such 
that the test rig could not be run higher in speed. The subsynchronous 
response of the first critical speed possibly was due to internal damping 
similar to that observed with the undamped test rig. If this is the case, 
then the effect of the added damping was to raise the threshold of this 
instability from just over one time the first critical speed to almost 13 
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Fig. 157 Log Decrement as a Function of Support Stiffness 
(Damping = 8760 N-s/m [50 lb-sec/in.]) 
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in Supercritical Power Transmission Shaft 



times the first critical speed, which was a substantial improvement. The 
sub synchronous response of the fourth critical speed was most significant 
at a rotational speed of about 11,800 rpra, as illustrated by the frequency 
spectrum plot at the top of Figure 159. However, at a rotational speed of 
12,000 rpm, the sub synchronous response of the first critical speed became 
predominant; that of the fourth critical speed had become much less signi- 
ficant, as illustrated by the frequency spectrum plot presented at the 
bottom of Figure 159. In any case, it was demonstrated. that the test rig 
could be run to more than 11,000 rpm with no large synchronous or nonsynchro- 
nous response of the test shaft. Additional details of the observed nonsyn- 
chronous effects noted with the power transmission shaft were provided by 
Darlow and Zorzi [8.15~J. 

The design of the fluorocarbon (Viton 70 - CFM) elastomers utilized in the 
supercritical shaft was based upon extrapolated data obtained from 0-ring 
testing during previous elastomer testing [8.14], and the stability analysis 
performed is shown in Figure 156. The supercritical shaft squeeze-film 
damper (Figure 160) was modified to accept replacement elastomer dampers. 

The design requirements included consideration of shaft disassembly and 
reassembly unbalance changes. The elastomer dampers were required to 
replace the squeeze-film damper without the necessity of shaft/damper dis- 
assembly. This permitted a direct comparison of elastomer and squeeze-film 
performance. The six fluorocarbon (Viton 70 - CFM) buttons shown in Figure 
161 were 0.635 cm (0.25 in.) in diameter and 0.635 cm thick. Six buttons 
provided a stiffness of approximately 7 x 10^ N/m (4,000 lb/in.) with a 
loss coefficient of 0.75 used for design purposes. Figure 162 shows an 
elastomer button as assembled in the squeeze-film damper housing, Figure 160. 

The elastomer consistently permitted higher operational speeds than the 
squeeze-film mount. Figure 163 illustrates the synchronous response measured 
when the power transmission shaft was balanced through the third critical 
speed with an elastomer damper. A change to the squeeze-film mount did not 
permit operation through the third mode. With the squeeze film operable, 
trim balancing of the first and third modes was continued, and the elastomer 
was again installed. Figure 164 illustrates the fact that the elastomer 
once again outperformed the squeeze-film mount, although for this run, the 
elastomer did not reduce vibration through the first mode as well as the 
squeeze film. 

Later balancing with the elastomer permitted operation to 13,000 rpm (Figure 
165), which is 1,000 rpm higher than that achieved with the squeeze-film 
Figure 159. However, as was the case with the squeeze film, the unstable 
subsynchronous excitation of the 16-Hz first critical was again the factor 
which limited higher speed operation of the power transmission shaft 
(Figure 166) . 

These results clearly demonstrate that elastomer dampers can control nonsyn- 
chronous, as well as synchronous, vibration of high-speed rotating machinery. 
Further, the elastomer damper used in this application continually outper- 
formed the squeeze film and permitted higher speed operation of the super- 
critical shaft. 
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Fig. 159 Frequency Spectrum Plots of Test Shaft Vibrations 
With a Damper Running Near 12,000 rpm 
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Fig. 160 Elastomer Damper/Squeeze-Film Damper Schematic 
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Fig. 163 Synchronous Response 
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Fig. 164 Synchronous Response 
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Fig. 165 Synchronous Response on Elastomer Damper 
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Fig. 166 Frequency Spectrum at 13,000 rpm With Elastomer Damper Operational 






8.3 T55 Helicopter Power Turbine Shaft '8.3; 


A turbine engine balancing rig was used as a test facility for this elastomer 
testing. The major hardware components of the test facility consist of a 
drive motor gearbox/belt drive speed increaser, vacuum chamber and computer- 
ized data acquisition and balancing system; see Figures 167, 168, and 169. 

The drive motor is a 45-kW (60 hp) variable-speed motor. Motor speed is 
varied from zero to 1700 rpm by a solid-state silicon-controlled rectifier 
circuit. The speed increaser provides approximately a 12:1 increase in 
speed via two components: a standard belt drive and a high-speed gearbox. 

The belt drive provides a 2:1 step-up in speed from the motor output to the 
input to the gearbox. The gearbox (Sunstrand Model LMV311) provides an 
additional 6:1 increase to the power turbine. The vacuum chamber reduces 
windage and permits operation with the low-power drive system. The automated 
data acquisition and balancing system is configured about a PDP 11/03 mini- 
computer with 32K words of memory and a RXV11 dual floppy disk drive for 
mass storage. 

The rig had successfully operated the T55 power turbine to the speed of 

16.000 rpm on its production hardware supports. Typical response data is as 
shown in Figure 170. Figure 170 is a trace of the synchronous response of 
an acceleration pattern for a 90° pair (vertical and horizontal probes) at 
approximately midshaft location. The two dominant features are the peaks 

at approximately midshaft location. The two dominant features are the peaks 
at approximately 4,000 and 6,000 rpm. It is significant that the horizontal 
peak at 4,000 rpm occurs when the vertical probe response is minimum. The 
vertical peak at 6,000 rpm shows similarity in the lack of response of the 
horizontal probe. This highly elliptical horizontal orbit rapidly changing 
to a dominant vertical ellipse is typical of rotor characteristics encouraged 
by asymmetric support characteristics. The response of the turbine end at 
these two speeds is dominant, whereas the roller bearing end (cold end) 
shows little activity indicating a strong precession of the turbine such as 
a rigid body conical mode shape. A subsequent set of peak occurs in the 

8.000 to 9,000 rpm speed range, with all probe sets along the shaft indicat- 
ing motion. The activity at this speed has not been predicted analytically, 
but is generally subordinate to the two main peaks at 4,000 and 6,000 rpm; 
see Figure 170. 

To assess possible rotor modes and the acceptability of the proposed rotor 
dynamic model, a shaker was attached to the CCAD rig to excite a nonrotating 
T55 power turbine. Swept sine-wave excitation, Figure 171, indicated three 
possible structural modes: 

• 115-Hz bounce mode of turbine 

9 210-Hz bending mode of shaft 

• 600-Hz bending mode of shaft 

Figure 172 indicates the location of these structural frequencies when com- 
pared to a tuned analytical model of the T55 (whirl speed marked by "X" at 
zero spin speed axis). This model, also predicted a first retrogarde 
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Fig. 170 Synchronous Response T55 on CCAD Rig 
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(backward) precision turbine conical mode at 4500 rpm, a first forward mode 
at 6000 rpm and a first bending (2nd critical) at 21*000 rpm. The produc- 
tion engine support rigidity needed to provide this model was reasonable* 
with a roller bearing stiffness of 8.75 x 10? N/m (500,000 lb/in.) and a 
ball bearing (turbine) pedestal stiffness of 1.4 x 10? N/m (80,000 lb/in.). 

Thus, it is reasonable to predict the modes at 4000 to 6000 rpm yet not 
encounter the bending mode with operation to 16,000 rpm. Gyroscopic forces 
affect the speed of the first backward and forward so that a minimum of 
asymmetry is necessary to induce excitation of both the 4500 rpm and 6000 
rpm peak. 

Further evidence that the location of the second mode for the baseline T55 
configuration was above the operating speed of 16,000 rpm was obtained from 
observing the supersynchronous excitation of the small 2/rev component of 
motion throughout a normal acceleration pattern to maximum speed. The 
2/rev peaked while operating at 7100 rpm indicating the presence of a mode. 
This is shown as a circled dot located on the dashed line extending from 
the lower left to the upper middle of Figure 122. The dashed line repre- 
sents a whirl speed of twice the corresponding spin speed (2/rev excitation) . 

Therefore, from observation of rotor synchronous response, static shaker 
excitation data and observation of the 2 /rev vibration, the rotordynamics 
model, Figure 173, was considered as acceptable for analytic examination of 
possible elastomeric-damper designs. 

Previous elastomer material testing of shear specimen using the Ling Shaker* 
was reviewed for candidate elastomer material selection. From Figures 76 
through 78 the fluorocarbon (Viton 70 - CFM) data for 32°C showed that it 
had the largest measured value of loss coefficient for the entire range of 
frequencies. As temperature increased, the selection of material for 
damping became somewhat arbitrary as the loss coefficient for all materials 
nitrile (Buna-N - NBR) , butadiene (Polybutadiene - BR) , EPDM, chloroprene 
(Neoprene - CR) and fluorocarbon (Viton - CFM) ranges between 0.1 and 0.2 
(Figure 76) for the entire frequency range. Therefore, the selection to 
use fluorocarbon (Viton - CFM) as the damping material was based upon the 
following: 

® Superior damping properties than the other materials 
tested for 32°C or 66°C operation (Figures 76 and 77). 

« Dissipation characteristics as good as any of the other 
materials tested for operation at 88°C (Figure 78). 

« Previous use as a damper in successful testing of high- 
speed rotor dampers. 


* Mechanical Technology Incorporated, Latham, N.Y. 
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The range of loss coefficient from data obtained from shear specimen shaker 
testing of fluorocarbon (Viton 70 - CFM) indicated that a range of loss 
coefficients from 0.15 to 0.75 was reasonable for the 0-16,000 rpm opera- 
tion of the CCAD balancing rig. Accordingly, the T55 power turbine rotor- 
dynamics model was modified to determine the "optimal" damper design. 

Figures 174 and 175. All system damping was assumed to be due to the 
elastomer damper; consequently no structural damping was modelled for the 
turbine or gearbox bearing. The two disc pack couplings were included as 
elastic elements which offered a minimum of bending rigidity but offered 
full shear restraints. Figure 173 indicates these couplings at stations 
8 and 19 of the rotor model. 

From the rotor damper optimization studies. Figures 174 and 175, a number 
of facts surfaced regarding the design of the elastomeric damper: 

$ The first mode is a conical precession of the turbine. 

The nodal point is extremely close to the elastomer- 
damper test bearing location. Accordingly, the first 
critical speed offers little damping and is insensitive 
to the elastomer-damper stiffness or damping character- 
istics. 

• The second and third mode shapes offer significant activity 
of the elastomeric damper. Consequently, the critical speed 
location and log decrement of those modes are sensitive to 
damper selection. 

$ A "best" or "optimal" tradeoff between the second and third 
mode is obtained by a fluorocarbon (Viton 70 - CFM) elasto- 
meric damper of between 5.25 x 10^ - 7.0 x 1C)6 N/m 
(30,000-40,000 lb/in.). 

« With this value of stiffness, the T55 rotor mounted on the 
elastomeric damper is expected to traverse the second 
critical, which is not within the operational range of the 
roller-bearing-mounted production T55 configuration, 

Figure 170. 

® The range of log decrement expected for the second critical 
speed is between 0.2 to 1.1 depending upon the loss coeffi- 
cient (temperature of operation) for the elastomer damper. 

The T55 elastomer-damper design offered uniquely challenging problems. This 
damper configuration was to be more compatible with actual gas turbine 
design technology and constraints. Buttons require added radial envelope 
and are not necessarily a desirable geometry for a bearing which may also 
be required to tolerate thrust. Most gas turbines do not have large amounts 
of radial envelope available, particularly if they are multishaft, turbofan, 
or turboshaft engines, such as the T55. Design maneuverability and flexi- 
bility are generally less restricted in the axial direction. The design must 
also accommodate possible thrust loading (although the T55 cold-end mount 
does not take thrust) in combination with radial loads and must be a design 
capable of controlled preload. These requirements were satisfied by a shear 
damper configuration. 
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Fig, 174 T55 Power Turbine — Critical Speeds and Log 
Decrement vs Bearing Stiffness for Fluorocarbon 
(Vi ton-70) Elastomeric Damper ( 77 = 0.15) 
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Fig. 175 T55 Power Turbine — Critical Speeds and Log 
Decrement vs Bearing Stiffness for Fluorocarbon 
(Viton-70) Elastomeric Damper (nrj- 0.75) 


CRITICAL SPEED (rpm) 



Figure 176 illustrates the concept used for the elastomer-damper mount on 
the roller bearing support for the T55. Two fluorocarbon (Viton 70 - CFM) 
shear rings mounted along the entire circumference satisfied all require- 
ments, Axial preload was accommodated by controlled machining of the outer 
flanges which attached to the housing by 12 equally spaced socket head cap 
screws. An overload protector was installed by using an 0-ring with a 
prescribed clearance 0.254 mm (0.010 in.) that would not interfere with 
normal operation of the damper but which would be available as a back-up 
protection if the shear damper failed, (As all testing proved successful, 
the 0-ring was never required to dampen rotor vibrations.) 


To obtain fluorocarbon (Viton 70 - CFM) shear damper with the required 
stiffness, the results of previous shear specimen data were used. Using 
the values for material properties and characteristic temperature for 
fluorocarbon (Viton 70 - CFM) and the unified design curves in Figures 
87 and 88, the analysis showed that an I.D. of 6.98 cm (2-3/4 in.) and an 
O.D. of 9.52 cm (3-3/4 in.) for 3.175 mm (1/8 in.) thick fluorocarbon 
(Viton - CFM) specimen would produce a stiffness in shear of approximately 
6.12 x 10^ - 7.0 x 10^ N/m (35,000 to 40,000 lb/in.) depending upon strain 
amplitude experienced in operation. 


Construction of the elastomeric damper. Figure 177, consisted of machining 
the aluminum bearing supports and shear restraining plates, bonding the 
fluorocarbon (Viton - CFM) using KEN H998* hardner with a KP-138 resin 
adhesive, and inserting four thermocouples and two capacitance probes into 
the assembly. The overload 0-ring was also attached using the REN adhesive, 
Figure 178. As shown, the damper was split radially into two halves for 
ease of assembly and replacement with the production T55 support (shown on 
the left of Figure 179). Prior to installation of the damper on the CCAD 
rig, the elastomer was subjected to static loading for calibration purposes. 
Figure 180 presents the static data for shear loading on the damper. The 
static values of stiffness determined were between 6.12 x 10^ and 6.47 x 
10^ N/m (35,000 and 37,000 lb/ in.). Static loads were applied and held 
for up to 5 minutes. No reduction in load was observed on the pressure 
load gauge, Figure 181. For the purpose of shear loading, a mandrel was 
machined to support the inner block on a n V n -block setup, Figure 182. In 
addition to shear loading, a static compression load was applied to the 
two 3.175 mm (1/8") elastomer pads. An axial compression stiffness of 
1.26 x 107 N/m (72,000 lb/in.) was determined. Figure 183. Again, holding 
load for a period of up to 5 minutes did not produce measurable load 
relaxation. 


With fabrication of the elastomer damper complete, the CCAD balancing rig 
with a balanced T55 power turbine was made available for this test. Figure 
184. Figure 185 illustrates the location and designation of the proximity 
probes which were active throughout this testing. After establishing a 
baseline condition for the T55 production configuration operating in the 
CCAD rig, the elastomer damper was installed, and the power turbine generated 
to 16,000 rpm. Imbalance was introduced at balance planes #1, #2 and #3 at 
various stages in testing ^8.:TJ, as shown in Figure 185. 
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Fig* 178 Cross Section of the T55 Elastomer Damper 
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Fig. 179 Production T55 Roller Bearing Support and the Replacement Elastomer Damper 



500 



(qi) avert 


329 
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Fig, 183 Axial Static Loading of Fluorocarbon (VRon-70) T5S Elastomer Mount 







The baseline run of the T55 power turbine supported on the production roller 
bearing housing showed the characteristics noted in Figure 170* Figures 
186 through 188 show the runout-compensated synchronous response for probes 
#1 through # 6, The data from Figure 186 for probes #1 and #2 indicates 
that the 4,000 and 6,000 rpm peaks are very small. These occur in Figure 

188 for probes #5 and #6 and are larger amplitudes, indicating strong 
precessing of the turbine end. The intermediate peaks occur for all probes 
in the 8,000-9,000 rpm range. For the production T55 engine power turbine 
operating in the CCAD rig, no critical speeds have ever been observed above 
the peaks recorded in the 8,000-9,000 rpm range, regardless of the imbalance 
conditions. 

The T55 roller bearing housing was removed and the elastomer damper housing 
was installed. Two of the four thermocouples monitored elastomer tempera- 
tures, whereas the other two thermocouples monitored outer race bearing 
metal temperatures. The two additional probes (probe #7 and #8 of Figure 
185) were installed with output recorded on the Sanagamo Sabre VI tape 
recorder. The T55 power turbine was then operated to 16,000 rpm. The 
filtered synchronous vibration response obtained is documented in Figures 

189 through 191. Comparative results from the baseline T55 run (Figures 
186 through 188) and the elastomer-damper test results (Figures 189 through 
191) are provided for probes #1 and #5 in Figures 192 and 193. From these 
results a number of facts emerged: 

• A substantial reduction in synchronous response occurs 
for the 6,000 and 8,000 rpm peaks with the elastomer 
operation (Figures 192 and 193). 

• That a reduction in vibration occurs for the entire 
16,000 rpm operational range for probe #5 (Figure 193). 

• A definite peak occurs (Probe #1) in the 13,500-14,000 
rpm range which has provided an increase in rotor vibra- 
tion with the elastomer damper operational (Figure 192) . 

By adding weight imbalance, confirmation of a second bending mode was obtained. 
In particular, the response of probes #7 and #8 of Figure 194 clearly indicates 
the presence of a mode critical speed at 13,500 rpm. 

Thermocouple readings on the elastomer indicated a temperature of approximately 
65° + 5.5°C (150° + 10°F) was typical for operation of this damper. This 
temperature was consistently within 3°C of the bearing metal temperature. 
Further, changes in axial thrust load of the thrust bearing did not alter 
significantly the response characteristic of the rotor. 

A comparison of the analysis performed for the elastomer-damper supported 
system with the measured data provides good correlation for mode shapes. 

Figure 195. Table 37 provides a comparison of critical speeds and logarith- 
mic decrements for the two criticals observed. 
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Fig. 187 T55 Power Turbine Baseline Data — Probes #3 and #4 
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Fig. 189 T55 Elastomeric Damper Baseline Data — Probes #1 and #2 
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Fig. 190 T55 Elastomeric Damper Baseline Data — Probes #3 and #4 


PROBE #5 
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Fig. 191 T55 Elastomeric Damper Baseline Data — Probes #5 and #6 
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Fig. 192 Vibration Plot Comparing Baseline and Elastomer-Mounted Rotor Response at Probe #1 



Fig. 193 Vibration Plot Comparing Baseline ana Elastomer-Mounted Rotor Response at Probe #5 
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Fig. 194 Response of T55 Elastomer-Damped Power Turbine 


T55 SCHEMATIC 



FIRST CRITICAL 



SECOND CRITICAL 



Table 37 


First Critical 
Second Critical 


Analysis 

Speed (rpm) Log Decrement 


6,022 0,001 

13,535 0,201 


Measured 


Speed (rpm) Log Decrement 


6,007 0.090* 

13,666 0.206* 


* Average value for all data obtained 


There is excellent correlation in whirl speed, but the first critical’s log 
decrement is significantly different in value between measured and analysis. 
This is due to the fact that analysis has indicated a node at the damper for 
the first critical. Slight changes in rigidity of the shaft or location of 
roller bearing support would alter appreciably the analytically predicted 
value of log decrement for the first critical. Consideration of structural 
damping which occurs at the turbine-end bearing (which is heavily participat- 
ing in the first critical (Figure 195)) was not included in the analysis. 
Therefore, the analysis does not reflect this form of damping which can 
significantly effect the log decrement of the first critical speed. 

Typical measured response indicated dominant 1/rev excitation of the rotor 
through the testing with the elastomer damper active. Figure 196 presents 
a typical frequency spectrum from probe #5 at 13,000 rpm. As shown, 
relatively minor subsynchronous and super synchronous excitations were 
experienced, and the operation of the T55 power turbine was as predicted 
throughout the entire test phase, 

8 . 4 Summary 

Three different elastomer dampers were designed and tested for three different 
test rigs. The three test rigs were a gas turbine simulator, a supercritical 
power transmission shaft, and a T55 power turbine shaft. The results of these 
tests clearly demonstrated the simplicity, convenience, flexibility and 
effectiveness of elastomer dampers. The detailed descriptions of the design 
procedures and actual hardware for these elastomer dampers can also serve as 
a valuable guide in the design of other elastomer dampers. 
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Fig. 196 Spectrum Plot on Elastomer Probe #5; 3 g at 300° Plane 1 
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